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ABSTRACT 


The cross-shore exchange of material is examined on beaches of varying morphology and 
hydrodynamics. On a dissipative, rip-channeled beach in Monterey, California, field 
measurements of the rip current vertical structure reveal depth-uniform flows inside the 
surf zone, and surface-dominated flows with rotation over depth outside the surf zone. 
Aperiodic, very low frequency pulsations in the near-surface velocity were found to 
control the mean rip current flow. Observations of Lagrangian surface drifters released on 
the inner shelf resulted in a new theory of material transport on a rip-channeled beach, 
where rip current pulses episodically expel material offshore through rip channels, and 
momentum-driven surfzone circulations pull the material back onshore over shoals to 
satisfy a mass balance, with no material being completely removed from the nearshore 
system. Advection and diffusion processes on a steep, reflective beach at Carmel River 
State Beach, California, are also investigated using measurements of the temporal and 
spatial evolution of dye released in the surf zone. A net offshore transport did exist, and 
the mixing and transport processes inside and outside the surf zone were independent. 
Overall, this work provides new insights into the fate of material transported between the 
surf zone and inner shelf on natural beaches. 
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offshore edge of a rip ehannel, just outside the surf zone (blue), 2) in an 
alongshore line outside the surf zone, spanning multiple shoals and rip 
ehannels (magenta), and 3) in a eross-shore line outside the surf zone, in 

line with a rip ehannel (green).105 

Figure 20. Conditions measured during the field experiment of (a) signifieant wave 
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depth. Gray shaded regions represent times of drifter deployments.106 

Figure 21. Drifter positions and speed, where the color of the line represents the 
speed, for drifter deployments exhibiting the locally-contained cross-shore 
exchange drifter pattern on yeardays 122, 129, and 132. Bathymetry 
contours are shown in the background in black, the shoreline is shown as 
the solid black line, and the approximate surfzone boundary is shown as 
the dashed gray line. Colorbar represents the drifter speed. Green circles 
show drifter release locations and red circles show drifter retrieval 
locations. Note there are different cross-shore and alongshore scales 

between plots.107 
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show drifter release locations and red circles show drifter retrieval 
locations. Note there are different cross-shore and alongshore scales 

between plots.108 
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relative drifter position displacements outside the surf zone for a relative 
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shore and alongshore exchange drifter pattern on yeardays 125, 126, 130, 

and 137. Contours are \ogw[P{rx', ry \ t'=2 s)] = -4.5, -4, -3.5, ... 1, 1.5.109 

Figure 24. Absolute diffusivity (one-particle) statistics for the cross-shore (Kxx, red) 
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alongshore exchange drifter pattern on yeardays 125, 126, 130, and 137.110 

Figure 25. Drifter cross-shore velocity magnitude a function of cross- 

shore distance, measured in surfzone-widths Lx, for each drifter 
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deployment, where onshore velocities shown in red and 

offshore velocities shown in blue, and error bars are shown. 

Black line represents the theoretical estimate of the Stokes drift averaged 
over the upper 0.5 m of the water column, Ma(0.5 m), which correspond to 
the depth of the surface drifters. Rows 1 and 2 show drifter deployments 
exhibiting the locally-contained cross-shore exchange drifter pattern on 
yeardays 122, 129, and 132, rows 3 and 4 show drifter deployments 
exhibiting the cross-shore and alongshore exchange drifter pattern on 
yeardays 125, 126, 130, and 137, and row 5 shows drifter deployments 

with no distinguishable exchange pattern on yeardays 135 and 136.Ill 

Figure 26. 3-hour mean vertical profdes, centered around low tide, of cross-shore 
currents measured at each ADCP, sorted by HJh and bin-averaged for 
HJh ranges given in the legend, where ADCPs 2, 3 and 4 were located at 
the surfzone boundary inside the rip channel, ADCPs 1 and 5 were located 
at the surfzone boundary over shoals, and ADCP 6 was located outside the 

surf zone offshore of the rip channel. Note the variable velocity scales.112 

Figure 27. Normalized depth-averaged cross-shore Eulerian currents ^ measured 

at each ADCP and bin-averaged by HJh, where ADCPs 2, 3 and 4 were 
located at the surfzone boundary inside the rip channel, ADCPs 1 and 5 
were located at the surfzone boundary over shoals, and ADCP 6 was 
located outside the surf zone offshore of the rip channel. Solid black line 
represents theoretical Stokes drift that balances ^ according to 

Lentz et al. (2008), and the dashed black line represents theoretical ^ 

through the rip channel that would balance the continuous onshore Stokes 
transport over the shoals according to the conventional rip current theory. ..113 
Figure 28. Drifter observations spatially binned into 10 m by 10 m bins and averaged 
for yearday 132, with the mean cross-shore movement in each bin 
represented as a percentage of time during the drifter deployment that the 
drifters were moving seaward. Red indicates drifters moved seaward 
100% of the time, blue indicates drifters only moved shoreward (moved 
seaward 0% of the time), and lighter colors indicate a tendency for drifters 
to move in both directions. Bathymetry contours shown in black and 

ADCP locations shown as black circles.114 

Figure 29. Eulerian velocity measurements at each ADCP (columns) during the 
drifter deployment on yearday 132, showing (top) the depth-relative mean 
vertical profile averaged over the duration of the drifter deployment, and 
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0.004 Hz) surface cross-shore velocity (black), representing VLE rip 
current pulsations, where ADCPs 3 and 4 were located at the surfzone 
boundary inside the rip channel. Note the variable velocity scales in the 

top plots.115 

Eigure 30. Carmel River State Beach, California, field site shown with a Google 
Earth image and bathymetry contours in meters overlaid in white. The 
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Figure 31. 


Figure 32. 


Figure 33. 
Figure 34. 


Figure 35. 


Figure 36. 


Figure 37. 


inset shows the mean eross-shore beaeh profile, relative to mean sea level. 
Red dots represent the alongshore array of stationary dye sensors, eyan 
squares represent the eross-shore array of ADCPs, and the green square 
represents the loeation of the EMCM array. The local coordinate system is 

shown, where positive is offshore and to the south.116 

Hourly-mean (a) root-mean square wave height, Hrms, (b) mean wave 
period, Tmoi, (c) depth-averaged cross-shore velocity, (d) depth-averaged 
alongshore velocity, and (e) tidal elevation measured at ADCPl in 3 m 
water depth. Red dots and error bars in (c) and (d) represent surfzone 
currents estimated with the dye. Green lines in (c) and (d) represent depth- 
averaged velocities measured by the EMCM array in 0.4 m water depth. 
Vertical dashed magenta lines represent times of dye releases. Velocities 
are in the local coordinate system where positive is offshore and to the 

south.117 

Photos showing examples of (a) the initial dye slug release, (b) the well- 
mixed surf/swash zone shortly after the dye release, and (c-e) a dye sensor 
attached near the sea-bed on a pole (shown by the green arrows) in the surf 
zone, highlighting the intermittency of the dye observations as the sensor 

is alternately submerged and not submerged.118 

Example of measured dye concentration as a function of time (black) and 

the dye envelope fit to the measurements (magenta).119 

Example of cross-shore velocities measured by the EMCM array for times 
when the EMCM array was (top) inside the surf zone and (bottom) outside 
the surf zone. Cross-shore velocity is in the local coordinate system where 
positive is offshore. The height of each EMCM sensor above the sea-bed 

is given in the legend.120 

(a) Water depth measured at ADCPl (dashed line) and at EMCM array 
(solid line), with times when the EMCM array was inside of the surf zone 
shown by the thick black line, (b) half-hour mean cross-shore velocity 
profiles measured by ADCPl (open symbols) and EMCM array (solid 
symbols), and (c) half-hour mean alongshore velocity profiles measured 
by ADCPl (open symbols) and EMCM array (solid symbols), as a 
function of depth-normalized height above the sea-bed. Symbols in the 
velocity profiles correspond to times in the water depth time series. 
Velocities are in the local coordinate system where positive is offshore 

and to the south.121 

Dye concentration as a function of time (left) and as a function of 
alongshore distance from the dye release (right) measured by the 
alongshore array of stationary dye sensors for (top) yearday 162a during 
weak alongshore current conditions, and (bottom) yearday 166c during 
strong alongshore current conditions. Colored circles in left plots represent 

tmoi computed with Equation (14).122 

Estimates based on the dye concentration time series measured by the 
alongshore stationary dye sensor array of (a) the surfzone averaged 
alongshore current V using Equation (13), (b) the alongshore decay of dye 
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mass M{y) using Equation (15), and (c) the surfzone averaged eross-shore 
eurrent U using Equation (12), as a funetion of alongshore distanee j from 
the dye release. In (b) the horizontal blaek dashed-dotted line represents 
where Mg deereases by 1/e, and the in (b) and (e) the vertieal blaek dashed 
lines represent the division between the near-field and far-field regions of 

the flow.123 

Eigure 38. Dye dispersion inside of the surf zone measured with the alongshore array 
of stationary sensors as a funetion of weight-mean eoneentration time tmoi- 
The dashed line represents the line of best-fit for all measurements on 
yearday 166, whieh has an = 0.44 and is signifieant in the 95% 
eonfidenee interval. The alongshore diffusivity Ky is equal to half of the 

slope of the best-fit line.124 

Eigure 39. Example of data eolleeted by swimmers for dye releases (left) 172a and 
(right) 172b. (top) Repeated eross-shore transeets performed by different 
swimmers, and (bottom) the eorresponding mean dye eoneentration 
profiles as a funetion of eross-shore distanee for 5 minute time windows. 
Eaeh eolor represents a different swimmer, and solid lines represent their 
mean eross-shore transeet with ± one standard deviation represented by 
shaded areas. In bottom plots, solid vertieal blaek lines represent 
approximate loeation of the shoreline, dashed vertieal blaek lines represent 
approximate edge of surf zone, and asterisks indieate the first moment of 
the dye eoneentration profile representing the eross-shore loeation of the 
dye plume jix- The loeal eoordinate system is used where positive is 

offshore and to the south.125 

Eigure 40. Dye dispersion outside the surf zone measured with the eross-shore 
swimmer transeets as a funetion of weight-mean eoneentration time tmoi, 
with error bars representing the 95% eonfidenee levels for the dispersion 
estimates for eaeh transeet. Solid lines of eaeh eolor represent the 
dispersion estimates from the swimmer data for that dye release. Dashed 
lines of eaeh eolor represent the lines of best-fit for the swimmer data for 
that dye release. Coeffieients of determination are given for eaeh 
swimmer, with all being signifieant in the 95% eonfidenee interval. The 
eross-shore diffusivity Kx outside the surf zone is equal to half of the slope 

of the best-fit line.126 

Eigure 41. A eoneeptual diagram of the mixing and transport of dye (pink) that was 
observed on the sandy, steep beaeh at CRSB. Dye released as an 
instantaneous point souree at the shoreline (star at t = 0, x = 0, y = 0), 
whieh is assumed to be rapidly vertieally well-mixed, initially dispersed in 
two dimensions in the near-field (t < 5 min, y < 50 m), where adveetion 
and diffusion proeesses were important; then when the dye was 
eompletely mixed aeross the surf zone, the dye was transported due to 
adveetion in the far-field (t > 5 min, y > 50 m). In the near-field, the dye 
spread in the eross-shore (kx) by turbulent diffusion due to breaking waves 
and spread in the alongshore (Ky) due to shear in the alongshore eurrent, 
while also being adveeted downstream with the surfzone averaged 
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alongshore current (V ) and also being advected offshore, resulting in a 
relatively large amount of cross-shore exchange (Uiy)). In the far-field, the 
dye was advected downstream with the surfzone averaged alongshore 
current {V ) and the cross-shore transport leveled off to a constant (U). 
Outside the surfzone (x > Xw), the dye moved as a constant patch, with a 


small amount of cross-shore spreading.127 

Figure 42. Frequency band-passed (top) cross-shore velocities and (bottom) 
excursion lengths computed from the time-integration of cross-shore 
velocity, measured at ADCPl for the given frequency bands.128 
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I. INTRODUCTION 


The transport of material, ineluding suspended sediment, pollutants, biological 
matter, and nutrients, across the surf zone and the continental shelf has a significant 
impact on beach composition, water quality, and ecosystem stability. Pollution in the 
ocean is an environmental hazard, degrading ecosystems and harming wildlife, as well as 
a human health hazard. Land-based pollution, including agricultural run-off and fecal 
matter, enters the ocean in the surf zone by means of run-off, storm drains, and coastal 
tributaries (Culliton 1998; Schiff et al. 2000; Grant et al. 2001), where it can then be 
carried offshore into the inner continental shelf. Nutrients are supplied to continental 
shelf ecosystems, which are some of the most productive on Earth, via cross-shore 
exchange between the surf zone and the inner shelf (Falkowski et al. 1998). Harmful 
algae blooms and offshore operational discharges, such as oil, originate on the inner shelf 
and are transported into the surf zone (Anderson 2009). Cross-shore exchange on natural 
beaches plays an important role in all of these processes; however, it remain s poorly 
understood. 

Onshore and offshore exchange occurs by various processes, depending on beach 
morphology, beach slope, wave conditions, and resulting current patterns. In this work, 
cross-shore exchange is investigated on a dissipative, alongshore-inhomogeneous, rip- 
channeled beach, and on a reflective, alongshore-homogeneous, steep beach. Dissipative 
beaches are defined as having low-sloping profiles and relatively wide surf zones, where 
waves break well offshore and continue to lose energy as they travel as spilling breakers 
onshore. The type of dissipative beach examined in this work is the transverse bar and rip 
state, where the transverse bars are relatively dissipative and strong rip currents occur in 
the deeper embayments. Dissipative, sandy beaches are popular recreation sites. 
Reflective beaches are characterized as having a steep (slope > 1/10), usually linear, 
beach face, and a considerably shallower sloping profile seaward of the step. Steep 
beaches typically have plunging to surging breakers and very narrow surf zones, 
sometimes confined to the zone of run-up on the beach face. Reflective beaches 
commonly include rocky shorelines where biology and organisms commonly dwell. 
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Understanding the transport and retention of hazardous pollutants and biologieal matter 
on these beaehes is important for reasons ranging from human health faetors to sustaining 
eeosystems. 

The eross-shore exehange of material due to the eomplex eireulation patterns on a 
dissipative, rip-ehanneled beaeh is investigated in Chapters II and III. Rip eurrents are 
defined as strong seaward-direeted flows and are generally eoupled to the underlying 
beaeh morphology, eharaeterized by deeper, eross-shore oriented, ineised ehannels in 
shallower, shore-eonneeted shoals or bars. Rip eurrents typieally oceur when waves 
approaeh the eoast shore-normally, and alongshore perturbations in the bottom 
topography induee alongshore variations in the depth-indueed wave breaking, leading to 
alongshore gradients in pressure and momentum driving water from the shoals toward the 
rip ehannels (Bowen 1969). A nearshore eireulation system exists where broader regions 
of wave-driven, onshore-direeted flows over the shoals transition to alongshore-direeted 
feeder eurrents near the shoreline, whieh eonverge at the base of the rip ehannel and flow 
seaward as a rip eurrent. Rip eurrents eommonly oeeur in some variation, on most 
beaehes, and have been observed in nature and studied extensively for many years (e.g., 
Shepard 1936; Shepard et al. 1941; Smith and Largier 1995; Brander and Short 2001; 
MaeMahan et al. 2005). However, it still remains largely unknown how these eurrents aet 
as a transport meehanism, both temporally and spatially, between the surf zone and the 
inner shelf. 

The vertieal strueture of the seaward-flowing rip eurrent ean have a signifieant 
effeet on the volumetrie estimates of exehange on rip-ehanneled beaehes, and is 
examined in Chapter II. Field measurements within the surf zone obtained by MaeMahan 
et al. (2005) indieate the rip eurrent vertieal profile is depth-uniform below the wave 
trough, with a signifieant amount of vertieal shear in the wave erest-trough area due to 
the presenee of waves and onshore Stokes drift eountering the offshore rip eurrent flow. 
These findings support the laboratory investigation performed by Haas and Svendsen 
(2002) of rip eurrent vertieal structure. Inside of the rip channel where waves were 
breaking, the rip currents were nearly depth-uniform, whereas outside of the surf zone 
there were strong depth variations in the flow, with strong seaward velocities in the upper 
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water column and weak seaward or shoreward velocities near the bottom. An extensive 
set of field measurements obtained at Monterey, California, of vertical profiles of rip 
current flows inside and outside the surf zone over a range of wave and tide conditions 
are presented in Chapter II, which supplement the very limited number of existing rip 
current profile measurements, particularly offshore of the surf zone. The variation of the 
magnitude and rotation of the rip current flows over depth are examined. 

Conventional rip current theory suggests that seaward flowing rip currents are the 
compensating return flow that balances the shoreward Stokes drift, which continuously 
transports water into the surf zone. Rip currents have long been believed to extend far 
offshore, and thereby be a means of continuously transporting material out of the surf 
zone (Shepard et al. 1941; Inman and Brush 1973). However, MacMahan et al. (2010) 
deployed a fleet of position-tracking surfzone drifters in different open coast beach rip 
current systems and obtained spatially synoptic estimates of rip current flow patterns. The 
results revealed semi-enclosed, large-scale vortices that were contained within the surf 
zone and largely retained drifters. Smith and Largier (1995) used a sector-scanning 
Doppler sonar to observe jets of water extending seaward from the surf zone (rip 
currents) and determined that the rip currents occurred episodically and recurred 
aperiodically, lasting several minutes each. The estimated flushing time due to the rip 
currents was on the order of 3 hours. 

Rip current surfzone retention and exchange evaluated with a three-dimensional 
(3-D) numerical model demonstrated the importance of including Stokes drift and very 
low frequency (VLF) motions to accurately describe the behavior of rip current exits 
(Renters et al. 2009; Reniers et al. 2010), where excluding Stokes drift or VLF motions 
resulted in an over or under estimation, respectively, of drifter exits (Reniers et al. 2009). 
These studies did not investigate the fate of the material (drifters) that was transported 
outside the surf zone by the rip currents. Talbot and Bate (1987) obtained qualitative 
aerial observations of dye released in rip currents and examined exchange rips that 
discharged dye offshore of the breaker line, which either travelled a short distance past 
the breaker line and was recycled into the surf zone within seconds, or moved up to 2 km 
beyond the breaker line and then moved alongshore and sometimes re-entered the surf 
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zone. MacMahan et al. (2010) found similar results, where drifters that did episodieally 
exit the surf zone were transported about 2 surfzone widths offshore from the shoreline, 
and oeeasionally some re-entered the surf zone; drifters that were retained within the surf 
zone were eoneentrated in the rip eurrent vortiees, and were less likely to be transported 
offshore due to the eoupling of the rip eurrent vortex to the underlying morphology. 
Reniers et al. (2010) used the 3-D numerieal model to produee Lagrangian Coherent 
Struetures (LCSs) within rip eurrent pulse surfaee veloeity fields. LCSs represent 
transport barriers hidden in the flow that eontrol the motion of passively adveeted 
partieles in unsteady flows. Examining the LCSs, the oeeasional drifter exit from the surf 
zone was shown to be due to VLF eddy motions, whieh detaeh from the nearshore rip 
eurrent eireulation eell when the VLFs extend beyond the surf zone. Therefore, VLFs 
were suggested to form the dominant exehange meehanism between the surf zone and the 
inner shelf. Field measurements of VLF pulsations in rip eurrent veloeity and their effeet 
on the magnitude and vertieal strueture of the mean rip eurrent flow are investigated in 
Chapter 11. 

Material transport and eross-shore exehange are investigated on a rip-ehanneled 
beaeh using field observations of position-traeking drifters and dye as traeers in Chapters 
111 and IV, respeetively. Observations of drifters released on the inner shelf of the rip- 
ehanneled beaeh in Monterey, California, are examined in Chapter 111. The fate of 
material that exits the surf zone is evaluated, ineluding the eross-shore extent and the on 
and offshore movement between the surf zone and inner shelf of the drifters, and it was 
found that no drifters eompletely left the nearshore region. A new theory of how material 
is exehanged within the nearshore region on a rip-ehanneled beaeh is presented in 
Chapter 111, where the mass transported offshore by the rip eurrent exits is balaneed by 
the surfzone rip eurrent eireulations over the shoals pulling material baek into the surf 
zone. Observations of dye released inside the surf zone on a rip-ehanneled beaeh in 
Destin, Florida, are examined in Chapter IV. In-situ measurements of dye eoneentration 
in a rip ehannel reveal VLF pulses in the seaward flowing rip eurrent, and aerial images 
of the temporal and spatial evolution of the dye reveal the adveetion and diffusion 
proeesses assoeiated with the rip eurrent eireulation system and that the dye is not 
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completely removed from the nearshore system. The eombination of the measurements 
are used to estimate surfzone flushing times, times for surfzone re-entry of material 
transported offshore, and mass transport quantities for a rip-ehanneled beaeh. 

The transport meehanisms and exehange behavior on a refleotive, alongshore- 
homogeneous, steep beaeh is examined in Chapter IV. There is almost no field data on 
how wave refleetion from the beaeh impaets the eross-shore movement of material. 
Previous dye traeer studies of surfzone mixing and material transport on wide, 
dissipative, alongshore-homogeneous beaehes with milder beaeh slopes (Harris et al. 
1963; Inman et al. 1971; Boehm 2003; Grant et al. 2005; Clarke et al. 2007) may be 
applieable to refleotive, steep beaehes. Dye released in the surf zone dispersed 2-D 
horizontally by adveotion and diffusion until the surf zone was saturated (Harris et al. 
1963; Inman et al. 1971; Clarke et al. 2007; Clark et al. 2010), then was predominantly 
adveoted by the alongshore ourrent and transported alongshore, beooming one¬ 
dimensional (1-D) in the alongshore direotion (Harris et al. 1963; Inman et al. 1971; 
Grant et al. 2005; Clarke et al. 2007; Clark et al. 2010). The alongshore volumetrie flux 
of water was found to be 50 times greater than the eross-shore volumetrie flux (Grant et 
al. 2005). The results suggest that most of the material remains within the surf zone and 
travels alongshore between the shoreline and aetive wave breaking with only a small 
amount of eross-shore exehange oeeurring eontinuously aeross the surfzone boundary. 
By evaluating the surfzone diffusivity, Clark et al. (2010) found that horizontal rotational 
motions with length-seales of the surfzone width or less are important for eross-shore 
traeer dispersion in the surf zone. Field measurements of surfzone mixing and transport, 
and eross-shore exehange aeross the surfzone boundary on a sandy, steep, refleetive 
beaeh at Carmel River State Beaeh, California, are deseribed in Chapter IV for a range of 
wave and alongshore flow eonditions. Fluoreseent Rhodamine dye was released as a slug 
in the surf zone, and the temporal and spatial evolution was measured. The adveetion and 
diffusion proeesses inside the surf zone and aeross the surfzone boundary are deseribed 
and estimates of the amount of eross-shore exehange are given. 

A summary of the findings from the field investigations of eross-shore exehange 
on natural beaehes is given in Chapter V. The mixing and transport proeesses on the 
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dissipative, alongshore-inhomogeneous, rip-channeled beach at Monterey, California, 
will be contrasted to those on the steep, reflective beach at Carmel River State Beach, 
California. Overall, the results from this work provide a better understanding of the cross¬ 
shore exchange between the surf zone and inner shelf on natural sandy beaches, including 
new insights into the fate of material transported offshore of the surf zone on dissipative 
and reflective beaches. 
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II. VERTICAL STRUCTURE OF RIP CURRENTS 


A. INTRODUCTION 

Rip currents are eommonly known to be a meehanism for transporting material, 
sueh as sediments, pollutants, biologieal matter, and even people, out of the surf zone. 
Numerous field studies have been performed foeusing on the temporal and spatial aspeets 
of rip eurrent transport (Shepard and Inman 1950; Sonu 1972; Talbot and Bate 1987; 
Brander and Short 2001; MaeMahan et al. 2010a; amongst others), with many reeent 
studies using Lagrangian drifters to examine the surfaee transport assoeiated with rip 
eurrents (Johnson and Pattiaratehi 2004b; Spy dell et al. 2007; Brown et al. 2009; Austin 
et al. 2010). In partieular, several field experiments were performed at the same field site 
on a rip-ehanneled beaeh at Monterey, California (MaeMahan et al. 2004a,b; MaeMahan 
et al. 2005; Brown et al. 2009; MaeMahan et al. 2010a), and numerous numerieal 
modeling efforts have foeused on simulating results from these experiments (Reniers et 
al. 2007, 2009, 2010), whieh have greatly inereased the understanding of rip eurrent 
dynamies. While the horizontal eireulation flow pattern and surfaee transport of rip 
eurrents is beeoming inereasingly better understood, limited field measurements of the 
vertieal strueture and sub-surfaee transport of the rip eurrent as it moves offshore and 
outside of the surf zone exist, whieh is eritieal for aeeurately deseribing the three- 
dimensional (3-D) material transport behavior of rip eurrent eireulations. 

In this work, field measurements of rip eurrent veloeities throughout the water 
eolumn, both inside and outside of the surf zone, obtained at the same field site at 
Monterey, California, are evaluated to examine the temporal and depth variation of rip 
eurrent flows. The results presented here will be given in eontext to the past work 
examining the role of very low frequeney (VLF) motions as the dominant meehanism 
responsible for rip eurrent transport out of the surf zone, as well as to previous field and 
laboratory measurements of the vertieal strueture of rip eurrents, whieh are discussed 
next. The goal of this work is to evaluate the presenee of VLF pulsations and their effeet 
on the mean rip eurrent flow, and to use the measurements of the rip eurrent vertieal 
strueture inside and outside the surf zone to deseribe how material is transported offshore 
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via rip currents. Ultimately, these findings will contribute to the existing knowledge of 
rip current exchange. 

B, BACKGROUND 

1. Rip Current Dynamics 

Rip currents are defined as strong seaward-direeted flows and are part of a 
dynamic nearshore eirculation system. On open coast beaches, they are generally coupled 
to the underlying beaeh morphology, characterized by deeper, cross-shore oriented, 
incised channels in shallower, shore-eonneeted shoals or bars (MacMahan et al. 2006). 
Morphologically controlled rip eurrents occur when waves approaeh the coast shore- 
normally, and alongshore perturbations in the bottom topography induce alongshore 
variations in the depth-indueed wave breaking, leading to alongshore gradients in 
pressure and momentum driving water from the shoals toward the rip channels (Bowen 
1969; Haller et al. 2002). A cireulation pattern exists where broad regions of wave- 
driven, onshore-directed flows over the shoals transition to alongshore-directed feeder 
eurrents near the shoreline, whieh converge at the base of the rip channel and flow 
seaward as a rip current. Mean rip current veloeities increase with increasing wave height 
(Shepard 1936; Shepard et al. 1941; Shepard and Inman 1950; MaeMahan et al. 2005), 
and are tidally modulated, with deereases in tidal elevation resulting in increases in rip 
current flows (Shepard and Inman 1950; Sonu 1972; Aagard et al. 1997; Brander and 
Short 2001; MacMahan et al. 2005). The interaetion between the seaward-flowing rip 
current and the incident waves has an influence on the offshore transport, where a 
negative feedback exists with the waves reducing the strength and offshore extent of the 
rip eurrent (Yu and Slinn 2003). An understanding of the vertieal strueture of the rip 
eurrent is necessary to understand the near-surface wave effects on the vertically varying 
rip current flow, whieh influences the vertical structure of the mean cross-shore current 
(Uehiyama et al. 2010). 

Spatially synoptie estimates of rip eurrent cireulation and Lagrangian transport by 
MacMahan et al. (2010a) provided field evidenee that rip eurrents are semi-enclosed, 
large-scale vortices that are largely contained within the surf zone. The rip eurrent 
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circulations resulted in a high amount of surfzone retention, with only episodie transport 
offshore (19% of drifters exited the surf zone per hour, on average). These observations 
were corroborated by drifter measurements in a laboratory experiment by Castelle et al. 
(2010) examining rip eurrent eireulations foreed by irregular, shore-normal ineident 
waves propagated over a moveable bed. Renters et al. (2009, 2010) found that VLF 
motions in the rip eurrent flow are the dominant meehanism responsible for the 
oeeastonal surf zone exit. These energetie motions for frequeneies less than 0.004 Hz 
have been observed in rip eurrent flows in the field on many beaehes (Smith and Largier 
1995; Brander and Short 2001; MaeMahan et al. 2004b), and are mostly eonlined to the 
surf zone eorresponding to surfzone eddies (SZEs) (MaeMahan et al. 2004b). Several 
mechanisms have been investigated for the generation of VLF motions inside the surf 
zone, ineluding shear instabilities of the rip eurrent (Haller and Dalrymple 2001; Renters 
et al. 2007; Bruneau et al. 2009) and wave group foreing (Renters et al. 2007; Long and 
Ozkan-Haller 2009). Using a numerieal model. Renters et al. (2007) input various types 
of wave foreing (monoehromatie, biehromatie, random waves) and found that foreing by 
directionally broad random waves best explained the high VLF intensity inside the surf 
zone and rapid decay outside the surf zone observed by MaeMahan et al. (2004b). The 
ineident wave energy appeared as random spatial blobs, and the spatially varying 
breaking waves on the wave group timeseale (0(10 minutes)) were responsible for 
foreing the VLF motions (Renters et al. 2007; MaeMahan et al. 2010a). The VLFs are 
outside the gravity restoring region in frequency-wavenumber spaee, indieating they are 
horizontal vortiees (i.e., SZEs). The VLF energy oeeurs at alongshore wave numbers 
eorresponding to the alongshore spaeing of rip ehannels, suggesting the SZEs are eoupled 
to the underlying bathymetry (MaeMahan et al. 2004b). Relaxation within the 
hydrodynamies allows the entire rip eurrent eells to oscillate in the eross-shore within the 
surf zone on VLF timeseales, but the rip eurrent eireulations are topographically 
controlled in the alongshore by the rip ehannel spaeing, therefore preventing them from 
migrating freely alongshore (MaeMahan et al. 2004b), as they would on a homogeneous- 
bathymetrie beaeh (Johnson and Pattiaratehi 2006; MaeMahan et al. 2010a). 
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The episodic transport of surface material out of the surf zone by a rip current was 
examined by Reniers et al. (2010) by evaluating Lagrangian Coherent Structures (LCSs) 
in the rip current surface velocity field produced by a 3-D numerical model capable of 
resolving wave groups. LCSs represent transport barriers in the flow that control the 
motion of passively advected particles. Model predictions of the fate of virtual surface 
drifters propagated by the Lagrangian velocity field showed that in the rip channels inside 
the surf zone, the virtual drifters clustered together and travelled seaward at the 
convergence of counter-rotating VLF SZE pairs. The presence of the VLF motions 
occasionally caused the rip current to pulse and extend beyond the surf zone. Smith and 
Largier (1995) observed well-defined jets of water extending episodically seaward from 
the surf zone, while Reniers et al. (2010) found that material trapped by the LCSs that 
become detached from the surfzone rip current circulation at the outer surf zone is carried 
offshore in narrow streaks. The field observations and numerical computations 
summarized here indicate the episodic, pulsating nature of rip currents and emphasize the 
importance of VLF motions in forcing rip current transport out of the surf zone. The 
presence and effect of individual VLF pulsations throughout the water column as the rip 
current moves offshore will be investigated in this work. 

2, Rip Current Vertical Structure 

The existing in-situ field measurements of the vertical structure of rip currents are 
limited, and were obtained primarily inside the surf zone (Aagard et al. 1997; Brander 
and Short 2001; MacMahan et al. 2005). Brander and Short (2001) used bi-directional, 
ducted impellor current meters to obtain vertical rip current velocity profiles in a rip neck 
and observed an increase in velocity away from the sea-bed and then a decrease towards 
the water surface, however measurements were made at only three locations in the water 
column. MacMahan et al. (2005) used a portable acoustic Doppler current profiler 
(ADCP) to obtain Eulerian rip current velocity measurements throughout the water 
column in a rip channel in the surf zone. They found that the vertical rip current profile 
was depth-uniform below the wave trough, with a significant amount of vertical shear in 
the wave crest-trough area due to the presence of waves and onshore Stokes drift 


10 



opposing the offshore rip eurrent flow. The rip current velocity profile as it moves 
offshore and outside of the surf zone has not been measured simultaneously in the field. 

A laboratory investigation of 3-D rip current flow was performed by Haas and 
Svendsen (2002), where rip currents were forced by normally-incident, monochromatic 
waves propagating over an alongshore bar with incised rip channels. The rip current flow 
was measured with three acoustic Doppler velocimeters (ADVs) positioned at different 
depths and moved to different horizontal locations during repeated wave runs. A similar 
laboratory experiment was conducted by Dronen et al. (2002) in a narrow wave tank, 
where half of an alongshore bar and rip channel system were represented. The laboratory 
results found that inside the surf zone in the rip channel where waves were breaking, the 
rip currents were nearly depth-uniform, whereas outside the surf zone there were strong 
depth variations in the flow, with strong seaward velocities in the upper water column 
and weak seaward or shoreward velocities near the bottom (Dronen et al. 2002; Haas and 
Svendsen 2002). The steady monochromatic wave forcing resulted in VLF motions 
outside the surf zone (Haas and Svendsen 2002) due to shear instability of the rip current 
jet (Haller and Dalrymple 2001; Reniers et al. 2009). This is a different forcing 
mechanism than surfzone wave-group induced forcing, which results in an aperiodic 
pulse that is transported beyond the surf zone, and the resulting rip current vertical 
structure may differ. Here, field observations of rip current vertical structure under 
directionally spread random wave conditions are examined. 

3, Rip Current Volumetric Transport 

The time variation of the vertical structure of rip current flows as they move 
seaward from the surf zone can have a significant effect on the volumetric estimates of 
exchange on rip-channeled beaches. Talbot and Bate (1987) used semi-quantitative 
measurements of dye releases to roughly estimate the flux associated with rip currents 
that exited the surf zone by assuming vertical uniformity everywhere, and found a half¬ 
residence time of 22 minutes for surf zone water. Smith and Largier (1995) observed 
episodic jets of water that extended seaward from the surf zone (i.e., rip currents) with a 
sector-scanning Doppler sonar. The rip currents lasted several minutes and transported 
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water rich in suspended sediments and bubbles offshore, suggesting they could cause a 
significant amount of exchange between the surf zone and offshore. An estimate of the 
surfzone flushing time was 0(3 hours) based on surface rip current velocities estimated 
by the sonar and assuming the flow spanned the entire water depth. Renters et al. (2009) 
estimated a surfzone flushing time of 0(5 hours) based on the assumption that the 
surfzone retention of surface drifters (80% retention per hour) decayed exponentially 
with time. These different estimates were based on different assumptions of the vertical 
structure of the rip current flow exiting the surf zone. Additionally, the retention/fiushing 
of material such as sediments, algae, larvae, bubbles and organic material are likely 
retained/flushed at the same rate, however their retention depends on vertical advection 
(Fujimura et al. 2014). Therefore, an understanding of the transport and vertical mixing 
associated with rip currents is important. 

C. FIELD EXPERIMENT 

1. Field Site 

Field observations were obtained at Monterey, California, from late April to late 
May 2007, during the Rip Current Experiment (RCEX). The field site was chosen owing 
to the year-round persistent rip currents and existing knowledge from previous 
experiments performed at the site (MacMahan et al. 2004a,b; MacMahan et al. 2005; 
Brown et al. 2009; MacMahan et al. 2010a). The beach was composed of a relatively 
steep foreshore (1/10 slope) with straight and parallel contours, flattening out inshore 
(1/100 slope), with quasi-periodic, 0(125 m), incised rip channels and continuing with 
a 1/20 offshore slope and straight and parallel contours seaward of the breaker zone 
(Eigure 1). Eeeder channels were present near the shoreline and converged with the 
incised rip channels. The local coordinate system used in this work is defined as positive 
offshore and to the south. 

2. Wave, Wind, Tide, and Bathymetry Measurements 

Offshore waves were measured with an upward-looking ADCP deployed in 

approximately 13 m water depth, which measured co-located pressure and velocity data. 

Hourly root-mean-square (rms) wave height, Hrms, in the sea-swell frequency band (0.04 
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</< 0.25 Hz) was computed by transforming the measured pressure speetrum to sea- 
surface elevation using linear wave theory (Guza and Thornton 1980). Hourly mean wave 
period, Tmoi, in the sea-swell band was computed as the first-order moment of the 
pressure energy density speetrum and the hourly mean wave angle, Omoi, in the sea-swell 
band was eomputed from the veloeity auto and eross-speetra (Herbers et al. 1999). 

Wind speed and direetion data was obtained onshore of the field site atop a 10 m 
tower and rotated into the loeal eoordinate system. Wind stress was ealeulated using the 
method by Large and Pond (1981). Tidal elevation data, refereneed to mean sea level 
(MSL), was obtained from NOAA/NOS tidal station 9413450, loeated approximately 3 
km south of the field site. Bathymetry was measured throughout the experiment using 
sonar-equipped personal watereraft (MaeMahan 2001), a walking person, and an all- 
terrain vehicle, all equipped with a global positioning system (GPS). 

During the 19 day long field investigation, the hourly estimated wave eonditions 
in 13 m water depth (Figure 2) ineluded Hrms of 0.2 to 1.2 m, T^oi of 6.3 to 12.8 s, and 
dmoi of -4.2 to 5.9°. Shore-normal ineident waves are represented as 6mo = 0°, with 
positive wave angles approaehing from the north. The eross-shore winds were generally 
onshore, with cross-shore wind stress, of -0.29 to 0.02 N m'^. Spring tides were 
about 2.5 m and oecurred near the end of the experiment. 

3, Rip Current Measurements 

In-situ measurements of the rip eurrent vertieal structure were obtained inside and 
outside the surf zone with two bottom-mounted ADCPs sampled continuously at 1 Hz {dt 
= 1 s). The ADCPs were deployed in a eross-shore line in a rip channel inside the surf 
zone (ADCPin) in approximately 2 m water depth and just offshore outside of the surf 
zone (ADCPout) in approximately 4 m water depth (Figure 1). The ADCPs inside 
(outside) the surf zone were mounted approximately 0.35 (0.5) m above the sea-bed, and 
measured veloeity in a total of 30 (20) measurement eells throughout the water eolumn, 
with a eell size of 0.2 (0.25) m, with the first veloeity measurement taken approximately 
0.65 (1.23) m above the sea-bed. 
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The ADCPs measured eo-loeated pressure and horizontal veloeity data. The time- 
varying water depth, h, at each ADCP was computed as the sum of the measured time- 
varying pressure, p, and the height of the instrument above the sea-bed. The elevation of 
each ADCP velocity measurement cell above the sea-bed, d, where J = 0 m corresponds 
to the sea-bed, was computed as the sum of the height of the instrument above the sea¬ 
bed, the measurement blanking distance from the sensor head, and the measurement cell 
size multiplied by the given cell number. 

Throughout the experiment, the location of ADCPm actually varied between being 
inside and outside the surf zone as a function of the tidal elevation and incident wave 
height. Times when ADCPin was inside the surf zone were determined using the local 
wave-breaking criterion HrmJh > 0.4 (Thornton and Guza 1982), which typically 
occurred during large waves and low tides (Figure 3). 

a. Depth-Relative Normalization of Velocity 

Velocity data obtained from the upward-looking ADCPs resulted in erroneous 
measurements above the sea-surface, with the number of measurement cells above the 
sea-surface varying with time due to the fixed size and number of measurement cells and 
the time-varying sea-surface elevation. In order to remove these erroneous data, the sea- 
surface elevation, rj, as a function of time was computed from the co-located pressure 
data using linear wave theory (Dean and Dalrymple 1984). Due to noise in the ADCP 
velocity measurements near the sea-surface, velocity data in measurement cells greater 
than the total water depth, d> h + rj, were set to zero. 

In order to compare velocity measurements obtained in different water depths 
(due to different cross-shore instrument locations and different times in the tidal cycle at 
a single location), the velocity data are presented as a function of depth-normalized 
height above the sea-bed {dih), where d/h = 0 at the sea-bed and d/h = 1 at MSL. 
Additionally, the horizontal velocity data were rotated into the local coordinate system, 
such that the cross-shore velocity, u, was shore-normal and positive offshore, and the 
alongshore velocity, v, was shore-parallel and positive to the south. 
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D, OBSERVATIONS 

1. Mean Rip Current Vertical Structure 

The variation in the vertieal structure of the mean rip current flow as it moves 
seaward ultimately controls the amount of volumetric exchange occurring between the 
surf zone and the inner shelf. The cross-shore variation in the mean rip current velocity 
profdes measured by ADCPin and ADCPout were examined by computing hourly mean u 
and V as a function of relative-depth dih (Figures 3 and 4). The mean cross-shore velocity 
measured at ADCPin when it was inside the surf zone was depth-uniform and was greater 
in magnitude than ADCPout at the corresponding time (Figure 3). The mean cross-shore 
velocity measured at ADCPout, and at ADCPin during times when it was outside the surf 
zone, showed variation throughout the water column with maximum offshore velocities 
at the surface. The cross-shore rip current flows, both inside and outside the surf zone, 
were strongest during large wave conditions and low tides (e.g., yeardays 124 and 125), 
consistent with previous observations (e.g., Shepard and Inman 1950; MacMahan et al 
2005). The mean alongshore velocity measured at both ADCPin and ADCPout was 
relatively depth-uniform (Figure 4), and was smaller in magnitude than the corresponding 
cross-shore velocity. The difference in vertical structure of u and v outside the surf zone, 
suggests the alongshore component of the flow is not entirely due to the rip current, but 
may indicate the presence of inner shelf flows. Additionally, the fact that the alongshore 
current measured at ADCPout is greater than ADCPin when both were outside the surf 
zone (i.e., yearday 131) suggest the presence of inner shelf flows driven by pressure 
gradients, other than the tide. 

2, Very Low Frequency Motions in the Rip Current Flow 

VLF motions in rip current flows were shown to be the dominant mechanism 
responsible for exchange between the surf zone and the inner shelf on rip-channeled 
beaches (Reniers et al. 2009, 2010). The presence and intensity of VLFs throughout the 
water column is examined in the ADCP measurements inside and outside the surf zone, 
in order to evaluate depth-variable rip current transport from the surf zone to the inner 
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shelf. The VLF eontribution to the total horizontal instantaneous veloeity (U = \u^ + v^ ) 
was examined with the rms VLF total veloeity magnitude, Urms.viF, defined by 

If 0.004 Hz i i i i 

(z) = {Zj) + S^^ {zj)df, (1) 

where Suui’z-, f) and Swi^, f) are the energy density speetra of the eross-shore and 
alongshore veloeities. Spectra were calculated from 3-hour records, with dt = \ s, to 
capture the VLF signature in the flow with statistical confidence, using a window length 
of 1024 s with a Hanning window and 50% overlap, resulting in 56 degrees of freedom 
and a frequency resolution of 0.00098. Throughout the experiment, the VLF energy 
appeared to increase with increasing wave energy, both inside and outside the surf zone 
(Figure 5). During energetic conditions, Urms.viF at ADCPin was constant throughout the 
water column, while Urms.viF at ADCPout was more surface dominated and lesser in 
magnitude (Figure 5). The reduction in the VLF rip current intensity as it moved offshore 
is consistent with previous observations (MacMahan et al. 2004b). 

For the rest of the analysis, the effect of VLF motions on the strength and vertical 
structure of the rip current flow is examined by evaluating velocities measured only 
during low tides, when rip current flows are strongest (e.g., MacMahan et al. 2006), and 
during weak cross-shore wind conditions (|t™”'^x| < 0.05 N m'^), to ensure only wave- 
driven flows are examined (Figure 2). A total of thirty instances met these criteria during 
the experiment. Velocity measured during 3-hour windows, centered on low tides, 
denoted 3HRLT from here on out, were evaluated in order to ensure approximate 
stationarity in the flow, and to accurately capture the 0(10 minute) VLF component of 
the flow. 

In order to describe the strongest rip current flows, which were depth-uniform at 
ADCPin when it was inside the surf zone, and surface-dominated when it was outside the 
surf zone and at ADCPout, the VLFs were examined at a given reference level near the 
surface. The reference levels used in this work corresponded to dih = 0.7 at ADCPm and 
dih = 0.8 at ADCPout, which were the uppermost measurement cells within the water 
column that contained velocity measurements for each entire 3HRLT, and were located 
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below the lowest wave trough for each given record to eliminate the influence of the 
surface waves. 

The mean Urms.viF over all 3HRLTs at the near-surface reference level was -“O.IO 
m s'* at ADCPin and decreased to -“O.Od m s'* at ADCPout- In comparison, the mean total 
velocity Umean ovcr all 3HRLTs at the same near-surface reference level was '^0.27 m s'* 
at ADCPin and -0.08 m s'* at ADCPout- Therefore, the VLF motions are energetic, both 
inside and outside the surf zone, and the VLF component of the surface rip current flow 
accounts for approximately 37% and 50% of the total energy at ADCPin and ADCPout- 

a. VLF Pulses 

In previous rip current field observations, the cross-shore and alongshore velocity 
time series measured in a rip channel inside the surf zone had significant VLF signals 
(MacMahan et al. 2004b). The presence of VLF pulsations in the measured rip current 
flow was examined by low-pass filtering if < 0.004 Hz) the cross-shore and alongshore 
velocities {uvlf, vvlf) throughout the water column measured by ADCPin and ADCPout 
for the 3HRLTs. Examples of the uvlf pulses for high, intermediate, and low VLF energy 
conditions are shown in Figure 6. The pulsating nature of the VLF rip current flow is 
apparent, where the uvlf intensity varies in time between strong seaward velocities (warm 
colors in Figure 6) and relatively small or zero cross-shore velocities (cool colors in 
Figure 6). At ADCPin the pulses can be seen extending throughout the entire water 
column, while at ADCPout the pulses are predominantly at the surface. The magnitude 
and number of the VLF pulses decreased with decreasing VLF energy. 

The maxima in uvlf measured in the near-surface reference levels are further 
investigated to evaluate the individual pulses in the rip current flow (Figure 7). For each 
3HRLT the total number of VLF pulses, Npuhe, was determined as the total number peaks 
in Uvlf, which is assumed to represent the total number of independent VLF observations 
for that 3HRLT. In order to evaluate the maximum effect of the VLF pulses on the rip 
current transport, only significant VLF pulses were considered, which are defined as the 
successive peaks in uvlf above a given threshold. The threshold was determined by 
considering uvlf measured in the near-surface reference level over all 3HRLTs, as shown 
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by the frequency distribution in Figure 7. The frequency distribution shows the uvlf 
signal is composed of a large percentage of small pulses and a small percentage of large 
pulses. The effect of these strong pulses is of interest here, therefore the threshold used to 
determine significant VLF pulses was two times the standard deviation above the mean 
[jl + lo]. This equated to 0.68 m s"' and 0.19 m s’' for ADCPin and ADCPout, which are 

approximately seven and five times greater than the associated mean Urms.viF- A 
significant pulse was described by the duration, Tsig_puise, defined as the time from when 
Uvlf exceeded the given threshold until it decreased to a value below the threshold, and the 
amplitude, Asig^puhe, defined as the maximum uvlf during the pulse duration (Figure 8). 

Examples of the VLF rip current pulsations measured in the reference level during 
high, intermediate, and low VLF energy conditions are shown in Figure 8. The VLF 
pulses were found to be aperiodic and varied in amplitude throughout each 3HRLT. It 
was also observed that the VLF pulses did not occur simultaneously at ADCPin and 
ADCPout, and the number of significant VLF pulses measured at ADCPin did not 
necessarily equal the number measured at ADCPout- During high VLF energy conditions 
there were numerous significant VLF pulses measured at both ADCPin and ADCPout, 
with many of the pulses lasting for several minutes. As the VLF energy decreased, the 
magnitude and number of significant VLF pulses decreased. For low VLF energy 
conditions, the magnitude of the VLF pulses at ADCPin and ADCPout were equal at times, 
and there were less significant VLF pulses. For example, on yearday 128 there was only 
one, relatively long significant VLF pulse measured at ADCPout- These examples 
demonstrate the randomness of frequency and magnitude of the significant VLF pulses 
that occur during all VLF energy levels. 

The magnitude, duration, and frequency of the VLF pulses in the cross-shore 
velocity measured inside and outside the surf zone were evaluated to examine the 
influence of the VLF pulsations on the mean rip current flow. The VLF rip current flow 
for each 3HRLT was described by the number of significant VLF pulses {Nsig,pulse), the 
mean amplitude of the VLF pulses {Amean.sig,pulse), and the fraction of the time that the 
VLF pulses were occurring, defined as the sum of the individual pulse durations divided 
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by the 3-hour record length ). The relationship of the VLF pulse 

characteristics to the offshore wave height is shown in Figure 9, which highlights the 
increased VLF pulse activity for the more energetic wave conditions. The magnitude of 
Amean.sig,pulse at ADCPin was always greater than that at ADCPout for each 3HRLT, 
consistent with the observed decrease in Ums.vLF outside the surf zone. In general, more 
VLF pulses occurred during larger wave conditions, at both ADCPin and ADCPout- 
However, the number of VLF pulses measured at each ADCP was typically not equal. In 
fact, during several 3-hour records no significant pulses were measured at ADCPin while 
they were at ADCPout- This may be attributed to the higher threshold set at ADCPin, 
which incorporates measurements from when ADCPin was located both inside and 
outside the surf zone, which filters out the smaller VLF pulses when ADCPin was outside 
the surf zone. Examining the duration of the significant VLF pulses, the fraction of each 
3-hour record that the pulses were occurring at ADCPin and ADCPout was generally 
equal, with pulses occurring about 10% of the time on average. 

b. VLF Vertical Structure 

At each ADCP, the vertical profiles of uvlf corresponding to the each individual 
significant VLF pulse that occurred during all 3HRLTs were combined and bin-averaged 
based on the individual pulse amplitude Asig,puise measured in the surface reference bin. 
This resulted in the VLF rip current pulse vertical profiles being sorted into similar 
magnitudes and profile shapes (Figure 10). At ADCPin, the uvlf profiles are 
approximately depth-uniform, and increased in magnitude with increasing Asig,puise- At 
ADCPout, the Uvlf profiles are surface-dominated, with an increase in magnitude of the 
surface flow with increasing Agig, pulse- This results in a greater amount of shear in uvlf 
throughout the water column for larger Agig^puise, as the near-bottom uvlf at ADCPout are 
statistically similar for all Agig^puise- 

The effect of the VLF rip current pulses on the magnitude and vertical structure of 
the mean rip current flow was also examined by sorting the 3HRLT profiles by various 
VLF pulse characteristics (Figure 11). Inside the surf zone, the VLF pulses added to the 
mean background rip current flow, while outside the surf zone there is no mean 
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background flow so the strength of the VLF pulses eontrols the strength of the rip eurrent 
flow. In all oases, the magnitude of the mean rip eurrent oross-shore velooity profile was 
only a fraction of the measured mfzf profiles. The maximum amplitude of the VLF pulses 
in the rip eurrent surfaoe velooity Amax.sig,pulse that ooourred during eaoh 3HRLT was found 
to distinguish the 3-hour mean rip eurrent strength and vertioal profile (Figure 11, oolumn 
1). Although the strength of the individual VLF pulses that ooourred during a 3-hour 
reoord varied, the strength of the maximum VLF pulse for that reoord, even if it only 
ooourred onoe, seemed to represent the strength of the oorresponding 3-hour mean rip 
eurrent flow. In general, the mean oross-shore velooity rip eurrent profiles below the 
wave trough were primarily depth-uniform at ADCPin, with stronger flows for times 
when larger Amax.sig,pulse oocurred. For smaller Amax.sig.puise, the rip eurrent profde weakened 
and beoame more surfaoe dominated, and ADCPin was likely outside of the surf zone 
during those times. At ADCPout the rip eurrent was always surfaoe dominated, however 
the flows were stronger and there was more shear in the velooity profile for larger 
Amax.sig,pulse- Tho mean Eulerian velooity profdes presented in Figure 11 show the velooity 
measurements between the wave trough and orest, where the velooity goes to zero in this 
region at ADCPin and there are some onshore velooities in this region at ADCPout- Similar 
trends were observed when the mean velooity profiles were sorted by the mean VLF 
pulse amplitude Amean.sig.puise of eaoh 3-hour record (Figure 11, oolumn 2). These results 
suggest that although a number of smaller VLF pulses ooour during eaoh reoord, the 
maximum VLF pulse amplitude ultimately eontrols the rip eurrent flow magnitude and 
profde. 

The rip eurrent velooity profiles were also sorted based on the total number of 
signifioant pulses {Nsig,pulse) and total signifioant pulse duration i'LT^ig^putjT^hour) 3- 

hour reoord (Figure 11, oolumns 3 and 4). However, for many of the reoords only one 
signifioant VLF pulse ooourred, with the total pulse duration being a very small 
percentage of the total reoord, therefore using these characteristics to sort the mean 
velooity profiles was not helpful. 
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c. VLF Rotation over Depth 

As a rip current pulse moves offshore and ehanges in magnitude, the direction of 
the pulse likely also changes. The temporal shifts between the VLF horizontal veloeity 
time series measured as a function of depth at each ADCP is determined by eomputing 
the cross-eorrelation coefficient of u measured at each depth with u measured in the 
surface reference level, and similarly for v, and evaluating R{t = Q) (Figure 12). The 

statistical significance of the computed eorrelation values was based on Npube for eaeh 
3HRLT, where only statistieally signffieant values at the 95% signifieance level are 
shown in Figure 12. During energetic conditions (e.g., yeardays 123 to 126), the u (v) was 

strongly correlated (i?(T = O) > 0.75) throughout the entire water eolumn at ADCPm, 
while the correlation decreased with depth at ADCPout, as expeeted. 

In order to evaluate the rotation in the VLF pulses over depth, the mean horizontal 
VLF veloeity vector at multiple depths for each 3HRLT were calculated (Figure 13). At 
ADCPin, minimal rotation was observed throughout the water eolumn (not shown), as 
expected, due to the depth-uniformity that has been observed there. Therefore, the rip 
eurrent pulses were correlated over the entire depth at ADCPin, with no twisting. At 
ADCPout, the surface UmeanjiF vectors were greater in magnitude and directed more 
straight offshore, whereas the bottom Umean,vLF vectors were smaller in magnitude and 
were offset in direction (Figure 13). The angle of rotation between the surfaee and bottom 
Umean.vLF vootors Varied, with the average difference in angle being ~50°. Therefore, as 
the rip current moved seaward, the flow near the bottom was less correlated with the 
surface flow, and there was twisting of the rip current pulse over the depth. Complex 
correlation analysis was also performed, whieh gives a measure of the correlation 
between the total horizontal veloeity (given as a eomplex number) and the angle of 
displacement between them over the depth (Kundu 1976). Similar results were found, 
however they are not shown here, as the eross-eorrelation and rotation of the veetors over 
depth provides more visual apparent results. 
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E, 


DISCUSSION 


1. Offshore Movement of VLF Rip Current Pulses 

The observations of the magnitude and vertical structure of the mean (Figures 3 
and 4) and rms VLF (Figure 5) velocities measured at ADCPin and ADCPout indicate a 
distinct difference in the rip current flow dynamics as it flows from inside to outside the 
surf zone. The relationship between the VLF pulses measured at ADCPin and ADCPout 
was examined by performing cross-correlation analysis of uvlf and vvlf individually in 
the near-surface reference levels for the duration of the experiment using 3-hour record 
segments. The maximum correlation values, (t) , between ADCPin and ADCPout are 

shown in Figure 14 for the cross-shore and alongshore velocity components. The 
significance level of the correlation values was determined based on Npube measured at 
ADCPin and ADCPout, where the minimum Npube measured for each 3HRLT was used to 
determine the corresponding degrees of freedom {dof= Npube - 1) and correlation value 
for the 95% significance level, which are shown in Figure 14. The cross-shore and 
alongshore velocities were not significantly correlated between ADCPin and ADCPout for 
any of the 3HRLTs, indicating the VLF pulses inside and outside the surf zone were 
independent. 

Previous work showed that the VLF velocity was related to the mean rip current 

flow, as well as the offshore wave height and tidal conditions (MacMahan et al. 2004b). 

This was evaluated here by performing linear regression analysis between Urms.viF and 

various flow and wave parameters averaged over each 3HRLT, for all thirty 3HRLTs 

(significance level of = 0.35 at the 95% significance level) (Figure 15). There was no 

significant correlation between Urms.viF measured at ADCPout and the offshore wave 

conditions of Hrms, Tmoi, and 9mn (Figure 15a, b, c), while Urms.viF measured at ADCPin 

was only correlated with Hrms (Figure 15a). Therefore, the VLF pulses inside the surf 

zone are dependent upon the incoming wave energy, due to the forcing of the VLF 

motions by the spatially varying breaking wave groups. However, the VLF pulses outside 

the surf zone are not dependent upon the offshore wave conditions, supporting the notion 

that they originate inside the surf zone but travel offshore randomly. The correlation 

between Urms.viF and the corresponding Umem.viF (Figure 15d) shows significant 
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correlation at ADCPm (r^ = 0.73), but no correlation at ADCPout- There was no significant 
correlation between Urms.vLF and the local water depth (i.e., indicative of tidal elevation ) 
at either ADCP. The highest degree of significant correlation for both ADCPs was 
Urms.vLF and the local HrmJh (r^ = 0.74 and 0.44 at ADCPm and ADCPout). Therefore, the 
strength of the VLF pulses is dependent upon both the wave energy and water depth. 
These results further suggest that the VLF pulses at ADCPout are independent of those at 
ADCPm. 

The uvLF and vylf measured throughout the water column at ADCPout was further 
examined to evaluate the VLF pulses as they move offshore. On yearday 125 (Figure 6), 
significant uvlf pulses near the surface are apparent, while vvlf was approximately zero 
(not shown), indicating the rip current pulsation was directed predominantly offshore. For 
some other 3HRLTs, the magnitude of the uvlf and vvlf pulses was near equal at the 
same time (not shown), indicating the rip current pulsation was moving offshore at an 
angle. Both of these cases suggest the rip current flow behaved like a jet, but with the rip 
current pulses changing direction as they exit the surf zone. However, on yearday 137 
(Figure 16) there was a variation in time in the strength of the uvlf and vvlf pulses, where 
at the beginning of the record on yearday 137 uvlf is dominant and vvlf is basically zero, 
then briefly uvlf and vvlf are approximately equal, and within minutes, vvlf is dominant 
and Uvlf is basically zero, with this pattern repeating once more. This successive change 
in magnitude and direction of the rip current VLF pulse suggests the fixed point ADCP is 
observing the passage of a detached VLF eddy. The range in magnitude and direction of 
the measured VLF pulses at ADCPout is consistent with field measurements that showed a 
large range in Lagrangian drifter trajectories exiting the surf zone via a rip current on the 
same beach, with drifters moving offshore at different angles and travelling varying 
distances (MacMahan et al. 2010a; Brown et al. 2014, submitted manuscript). The 
numerical modeling results of Renters et al. (2010) showed similar virtual drifter 
trajectories due to the presence of wave-group induced VLF motions temporarily 
extending the rip current flow beyond the surf zone, resulting in the corresponding VLF 
eddies becoming subsequently detached. 
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2 . 


Surface Transport and Vertical Mixing by Rip Currents Outside Surf 
Zone 


The movement of the rip eurrent pulse as it moves at the surface out of the surf 
zone and offshore is further investigated. It is hypothesized that the rip current moves 
offshore as a jet at the surface, and its movement over the surfzone break (defined as the 
abrupt change in depth at the surfzone boundary) as it exits the surf zone results in the 
surface dominated flow offshore. For 3HRLTs with strong VLF pulses at ADCPm, the 
corresponding 3-hour mean cross-shore velocity profdes at ADCPm and ADCPout, as a 
function of vertical elevation, z, are examined, where z = 0 corresponds to MSL and is 
positive upward (Figure 17a). The cross-shore velocity profiles are shown relative to the 
mean cross-shore bathymetry profile in the rip channel. At ADCPin, the rip current 
velocity is depth uniform, extending to the sea-bed. At ADCPout, the surface-dominated 
rip current appears to be moving offshore as a jet. From a maximum near the surface, the 
velocity decreases downward with a large amount of shear to a certain depth, and then 
appears to decrease less rapidly or remain constant. The depth corresponding to this 
change in shear in the velocity appears to be the depth of the surfzone break. The jet-like 
behavior of the seaward flowing rip current measured here supports the observations of 
Smith and Largier (1995) of well-defined jets of water extending offshore from the surf 
zone as rip currents. 

The vertical mixing of material carried offshore by a rip current jet is further 
investigated with observations of dye releases that were performed at the same field site 
at Monterey, California, in June 2010 (Figure 17b-c). Fluorescent Rhodamine dye was 
released at a point inside of the surf zone and the distribution over time and space was 
measured. In particular, the vertical structure of the dye plume outside of the surf zone 
was examined by performing vertical casts from a GPS-equipped vessel with a co-located 
dye sensor and CTD throughout the deployment and at various locations. The casts 
extended to the sea-bed, and the measurements from the downcast and upcast were 
averaged to obtain mean profiles of dye concentration and temperature with depth (Figure 
17b-d). The temporal and spatial variation of the dye profiles outside the surf zone was 
examined to evaluate the material transport of the rip current jet offshore. Note that the 
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dye profiles measured outside the surf zone were from randomly ehosen locations 
between the surfzone boundary and several surfzone widths offshore, and at variable 
alongshore locations offshore of rip channels and shoals. Therefore, each profile may not 
correspond to the same rip current pulse from the same rip channel, so the evolution of 
one rip current pulses is not being described, rather the evolution of all material that is 
transported offshore due to rip currents. 

In general, the measured profiles revealed that the offshore dye plumes were 
relatively deep, encompassing greater than half, or all, of the water column. Near the 
surfzone boundary, and within the first 30 minutes after the dye was released (Figure 
17b), the vertical profiles of the dye were surface dominated, within the depth range of 
the depth at the edge of the surf zone. Later in time (Figure 17c), the vertical profiles 
taken closer to shore (200 < x < 275 m) became more vertically mixed within the depth 
range of the depth at the edge of the surf zone, while the vertical profiles farther offshore 
were still nearer to the surface and had not been mixed down yet. After an hour (Figure 
17d), the dye was well offshore of the surf zone and encompassed approximately half of 
the water column. At these cross-shore locations well offshore of the surf zone (x > 275 
m), the vertical profiles of dye were also observed to be coincident with the thermocline 
(Figure 17d). More mixing occurs in the upper water column offshore, resulting in 
warmer water and dye vertically mixed downward. At depths below the thermocline, the 
dye is no longer mixed down and the water temperature decreases. These results suggest 
that as the vertically-uniform rip current exits the surf zone, it moves offshore as a 
surface jet with its vertical extent remaining the same as at the edge of the surf zone. As 
the water depth increases outside the surf zone, the rip current transport remains near the 
surface, resulting in a surface dominated rip current flow offshore. The material 
transported offshore due to rip currents is then mixed down vertically to a depth 
coincident with the thermocline. The surface material transported by the rip current 
outside the surf zone, which encompasses a large portion of the water column, spreads 
laterally and can have a significant effect on the overall 3-D material transport. 
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F. 


SUMMARY AND CONCLUSIONS 


Field measurements of the vertical structure of a rip current, inside and outside the 
surf zone, were obtained on a rip-channeled beach at Monterey, California. The 
observations presented here provide new insight into the temporal and depth variation of 
a rip current as it moves offshore, and contribute to the previous findings of horizontal rip 
current circulation dynamics, also obtained at the same field site (MacMahan et al. 
2004a,b; MacMahan et al. 2005; Brown et al. 2009; MacMahan et al. 2010a). In these 
measurements, wave-group induced VLF motions were found to be a significant 
component of the rip current flow, accounting for 37% to 50% of the total energy inside 
and outside the surf zone. The counter-rotating SZEs associated with these VLFs, which 
are coupled to the underlying bathymetry and oscillate in the cross-shore direction, 
occasionally converge in the rip channels and these VLF pulses cause the rip current to 
extend outside the surf zone (MacMahan et al. 2004b; Renters et al. 2009, 2010). The 
VLF pulses in the near-surface rip current flow, both inside and outside the surf zone, 
were found to be aperiodic, with the number and amplitude of significant VLF pulses 
increasing with increasing wave energy. Evaluating the VLL signal in the surface rip 
current flow during 3-hour time windows centered around low tide, there are times when 
zero or only one significant VLL pulse was measured, but during larger wave conditions 
up to 20 significant VLL pulses occurred, where a significant pulse is defined as being 
greater than ji + la , which equated to 0.7 and 0.2 m s’' inside and outside the surfzone. 
On average, the significant VLL pulses occurred for approximately 10% of the time. The 
magnitude and vertical structure of the mean rip current flow was controlled by the 
strength of the surface VLL pulses, where inside the surf zone the VLL pulses were found 
to add to the strength of the background rip current flow, and outside the surf zone, where 
there was no background flow, the stronger VLL pulses resulted in stronger rip current 
flows and more shear in the vertical. 

As the rip current moved from inside the surf zone to offshore, the significant 
VLL rip current pulses at ADCPin and ADCPout did not occur simultaneously, and there 
was no significant correlation in the cross-shore and alongshore velocity measured at the 
two locations. There was no significant correlation between the VLL energy measured at 
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each ADCP and the offshore wave conditions, except for with the local HrmJh (r^ = 0.74 
and 0.44 at ADCPin and ADCPout), indicating the VLF pulses were dependent upon both 
the wave energy and water depth, and that the VLF pulses occurring at each ADCP was 
random and the VLF pulses at ADCPout are independent of those at ADCPin. 

The VLF pulses outside of the surf zone behave like a surface jet moving directly 
offshore or at an angle at times, with uvlf either dominating or being equal to vvlf at the 
surface. Other times the VLF pulse has detached and moved as an eddy passing by the 
stationary sensor, with uvlf and vvlf at the surface successively being dominant. The lack 
of correlation between inside and outside the surf zone (Figure 14) also suggests that the 
VLF pulse detaches from the surf zone and moves offshore. The vertical profiles of dye 
that were measured outside the surf zone indicate the rip current pulses encompass a 
significant portion of the water column, and therefore the VLF pulses in the surface rip 
current flow result in a significant amount of cross-shore exchange. Offshore, the dye 
profiles were coincident with the thermocline, which thereby limited the depth of the 
vertical mixing of material in the surface rip current flow. The vertical profiles measured 
by the ADCPs indicated that inside the surf zone the VLF pulses were depth-uniform 
with minimal rotation throughout the water column, while outside the surf zone they 
became surface-dominated as they moved offshore (Figures 5 and 6), and there is more 
rotation throughout the water column, with the near-bottom flow not correlated with the 
strong near-surface flow (Figures 12 and 13). These results suggest that as the depth- 
uniform rip current exits the surf zone and moves over the surfzone break, it behaves like 
a surface jet with its vertical extent remaining the same depth as the depth at the offshore 
edge of the surf zone, and remains near the surface controlled by the depth of the 
thermocline as it continues to be transported and mixed offshore. 
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III. DRIFTER MEASUREMENTS OF CROSS-SHORE 
EXCHANGE ON A RIP-CHANNELED BEACH 


A. INTRODUCTION 

The cross-shore exchange of material across the surfzone boundary, in both the 
seaward and shoreward directions, plays an important role in processes ranging from how 
land-based pollution, which enters the ocean in the surf zone, is carried offshore to the 
inner shelf, to how harmful algae blooms and offshore operational discharges, such as oil, 
which originate in the inner shelf, are transported into the surf zone. Understanding the 
transport and retention of hazardous pollutants and biological matter in the nearshore 
region is important for reasons ranging from human health factors to sustaining 
ecosystems. 

Studies of cross-shore circulation over the inner shelf, focusing on wave-driven 
flows (i.e., during weak wind forcing), demonstrate the importance of the onshore 
Lagrangian Stokes drift, ust{z), associated with shoreward propagating surface gravity 
waves (Stokes 1847), and the compensating wave-driven offshore Eulerian flow, mA^), 
referred to as undertow. The presence of wave-driven undertow inside the surf zone is 
widely known and has been measured extensively in the field (Haines and Sallenger 
1994; Garcez Faria et al. 2000; Reniers et al. 2004b). Wave-driven undertow has also 
been found to extend well offshore of the surf zone (Lentz et al. 2008; Kirincich et al. 
2009; Ohlmann et al. 2012), and to be the dominant component of the depth-averaged 
Eulerian cross-shore flow during weak wind conditions. The changes in the distribution 
of ust{z) and u^z) throughout the water column as the water depth decreases from the 
offshore waters of the inner shelf to the shallower waters near the surfzone boundary was 
shown to control the net cross-shore transport between the inner shelf and the surf zone 
(Lentz et al. 2008; Ohlmann et al. 2012). An in-depth description of the theory of Stokes 
drift and undertow, including typical profdes inside and outside the surf zone, is given in 
the next Section. 

On the inner shelf of Martha’s Vineyard, Massachusetts, and Duck, North 

Carolina, during weak wind conditions and along-shelf uniform flows, Lentz et al. (2008) 
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found that the onshore Stokes drift opposed the offshore Eulerian undertow flow 
throughout the water eolumn, resulting in zero net Lagrangian transport at all depths. 
Expanding upon these findings, Ohlmann et al. (2012) released drifters during weak wind 
eonditions on the inner shelf off the eoast of two beaehes in Southern California, whieh 
are relatively alongshore uniform and typieally do not support rip eurrents, and observed 
a deeeleration of shoreward-moving surfaee drifters as they approaehed the surfzone 
boundary. This was attributed to ehanges in the relative shapes of the onshore Stokes drift 
and offshore wave-driven undertow profiles with deereasing water depth near the 
surfzone boundary, resulting in an imbalanee of ust{z) and ue{z) with a net offshore 
Eagrangian flow near the water surface. The reason for the initial onshore motion of the 
drifters could not be discerned by Ohlmann et al. (2012). The results of Ohlmann et al. 
(2012) suggest that surface waves drive offshore Eulerian flows in the upper water 
column that may prevent the onshore transport of material into the surf zone from the 
inner shelf. 

These findings of cross-shore transport over the inner shelf on relatively 
alongshore homogeneous beaches are used as a starting point to investigate cross-shore 
exchange on rip-channeled beaches, which are known to have different dynamics, 
including the fate of material transported offshore from the surf zone to the inner shelf, 
which has never been thoroughly investigated. Rip currents are known to be capable of 
transporting material offshore, and have even been confused with wave-driven undertow 
flows (Shepard 1936), so it is conceivable that the interaction of the onshore Stokes drift 
associated with surface gravity waves and the offshore rip current flows will influence 
the net cross-shore transport on a rip-channeled beach. The conventional theory of rip 
current transport and an overview of previous rip current studies are given in the next 
Section. 

In this work, an extensive set of Eagrangian and Eulerian field measurements of 
surf zone and inner shelf rip current flows is analyzed to evaluate cross-shore exchange 
on a rip-channeled beach at Monterey, California. The observations are unique in that 
they not only describe the surf zone circulations, but they also explain flows outside of 
the surf zone, which have not been extensively measured before on a rip-channeled 
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beach. The results of these observations are presented, and the contribution to previous 
work evaluating rip current behavior is discussed. The findings of this study support a 
new theory of how material is exchanged and conserved within the nearshore region on a 
rip-channeled beach. 


B, THEORY 

1. Wave Forcing 

As surface gravity waves propagate onshore from the inner shelf to the surf zone, 
the time-averaged particle velocity in a wave results in a net onshore drift in the direction 
of wave propagation (Stokes 1847), referred to as Stokes drift. In a Lagrangian (particle 
following) reference frame, Stokes drift is depth varying (Figure 18a), where 


Us, (z) 


Hlcok cosh\_2k{z + h)~\ 
16 sinh^ kh 


cose„. 


( 2 ) 


where z is vertical elevation defined as positive upward relative to mean sea level (MSL), 
Hs is significant wave height, co is wave frequency, k is the dominant wave number, h is 
local water depth, and 6w is the mean wave direction (Stokes 1847). In a Eulerian (fixed 
in space) reference frame, the onshore volume Stokes transport occurs between the wave 
trough and crest (Figure 18a), because at a fixed point there is either no water above the 
wave trough or an onshore wave-averaged Stokes drift above the wave trough. In both 
reference frames, the onshore volume Stokes transport, Qst, is the theoretical depth- 
integrated ust{z), and is given by 


Qs, = —^ cos d.. 


16c ( 3 ) 

where g is gravitational acceleration, and c = co/k is the wave phase speed (Stokes 1847). 
In wave environments, a compensating Eulerian offshore flow, ue{z), exists below the 
wave trough, referred to as undertow (Haines and Sallenger 1994; Garcez Faria et al. 
2000; Reniers et al. 2004b). The distribution of ue{z) throughout the water column is a 
function of the local water depth, h. The summation of ue{z) and ust{z), at a given 
location, results in a wave-averaged Eagrangian velocity, ul{z), distributed throughout the 
water column, and is given by 
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Ul[z) = u^[z) + Us,{z). (^ 4 ) 

Material transport between the inner shelf and the surf zone is deseribed in terms of the 


Lagrangian transport, Ql, whieh equals the depth-integrated ul{z) and is defined as 

QL~QE'^Qst' ( 5 ) 

where Qe is the offshore Eulerian volume transport, equal to the depth-integrated ue{z), 
and estimates of the depth-averaged veloeity, or are found by dividing Qg or Qst 
by the loeal water depth h. 

On the inner shelf, in a Lagrangian referenee frame, the Stokes drift is largest at 
the water surfaee and deeays exponentially with depth (Figure 18b). During weak wind 
eonditions and alongshore-uniform flows, the eompensating wave-driven undertow is 
distributed as predicted by Hasselmann with a maximum offshore flow at the surface and 
decreasing toward the bottom (Figure 18b) (Hasselmann 1970; Lentz et al. 2008; 
Fewings et al. 2008). The wave-driven Eulerian undertow profile cancels the Lagrangian 
Stokes drift, which results in a vertical balance between ue{z) and ust{z), so ul{z) is zero at 
all depths and there is a zero net Lagrangian transport Ql (Hasselmann 1970; Xu and 
Bowen 1994; Smith 2006; Monismith et al. 2007; Lentz et al. 2008; Hendrickson and 
MacMahan 2009). 

Inside the surf zone, current dynamics are predominantly driven by breaking 
surface gravity waves and are largely dependent on the surf zone bathymetry. As the 
water depth decreases, the Lagrangian Stokes drift profile becomes more vertically 
uniform, and is constant throughout the water column shoreward of wave breaking 
(Figure 18c). In a surf zone characterized by an alongshore-homogeneous morphology, 
the onshore Stokes transport is balanced by an undertow flow with a parabolic velocity 
profile (Figure 18c) with maximum offshore velocity at mid-depth or near bottom 
(Svendsen 1984; Putrevu and Svendsen 1993; Haines and Sallenger 1994; Garcez Faria et 
al. 2000; Renters et al. 2004b). Based on Equation (4), this results in a vertical imbalance 
between ugiz) and ust{z) and a net onshore Lagrangian transport near the surface and 
offshore Lagrangian transport below the wave trough. In decreasing water depths from 
offshore to the surfzone boundary, changes in the relative shapes of the onshore 
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Lagrangian Stokes drift and offshore wave-driven Eulerian undertow profiles predieted 
by Ohlmann et al. (2012), resulted in an imbalanee between ue{z) and us^z) throughout 
the water eolumn, and an offshore Lagrangian flow near the water surfaee, eonsistent 
with the deeelerating shoreward-moving drifters on the inner shelf as they approaehed the 
surfzone boundary. 

2. Rip Currents 

Rip currents are strong, seaward-directed flows that originate near the shoreline 
and are generally coupled to the underlying beach bathymetry, characterized by deeper, 
incised channels in shallower, shore-connected shoals or alongshore bars. Rip currents 
have long been believed to extend far offshore, well beyond the surf zone and wave 
breaking, and thereby be a means of continuously transporting material offshore (Shepard 
et al. 1941; Inman and Brush 1973). On a rip-channeled beach, conventional theory 
suggests that seaward flowing rip currents are the compensating return flow, Qg.rip, that 
balance the continuous transport of water shoreward and into the surf zone by Stokes drift 
(Figure 18d). Therefore, at the surfzone boundary, a balance exists where the onshore 
Stokes transport Qst integrated over the rip channel spacing, Xrip, would equal the offshore 
Eulerian rip current transport exiting through the rip channel, defined as 

= Qp b , 

^St rip ^E,rip rip^ 

where bnp is the width of the rip current. Rip currents commonly occur in some variation 
on most sandy beaches, and have been observed in nature and studied extensively for 
many years (e.g. Shepard 1936; Shepard et al. 1941; Bowen 1969; Smith and Largier 
1995; Brander and Short 2001; MacMahan et al. 2005; MacMahan et al. 2010a). 
However, it still remains largely unknown how rip currents act as a transport mechanism, 
both temporally and spatially, between the surf zone and the inner shelf Most rip current 
studies focused on surfzone measurements, and very few field observations of exchange 
across the surfzone boundary have been obtained. 

Recently, MacMahan et al. (2010a) deployed a fleet of position-tracking surface 
drifters in different open coast beach rip current systems and obtained spatially synoptic 
estimates of rip current flow patterns, which revealed semi-enclosed, large-scale vortices 
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that were eontained within the surf zone and largely retained drifters (only 17% drifter 
exits per hour). Rip eurrent surfzone retention and exchange evaluated with a three- 
dimensional (3-D) numerical model demonstrated the importance of including Stokes 
drift and very low frequency (VLF) motions to accurately describe the behavior of rip 
current exits (Reniers et al. 2009; Reniers et al. 2010), where excluding Stokes drift or 
VLF motions resulted in an over or under estimation, respectively, of drifter exits 
(Reniers et al. 2009). Smith and Largier (1995) observed rip currents that episodically 
transported water offshore with a sector-scanning Doppler sonar in La Jolla, California, 
and concluded that rip currents could cause significant exchange between the surf zone 
and offshore, capable of flushing the surf zone in roughly three hours. 

These studies did not investigate the fate of the material (drifters) that was 
transported outside the surf zone by the rip currents. Talbot and Bate (1987) obtained 
qualitative aerial observations of dye released in rip currents and examined “exchange” 
rips that discharged dye offshore of the breaker line, which either traveled a short 
distance past the breaker line and was recycled back into the surf zone within seconds, or 
moved up to 2 km beyond the breaker line and then moved alongshore and sometimes re¬ 
entered the surf zone. Rip current flows outside the surf zone were investigated by 
Reniers et al. (2010) using a numerical model and unveiling Lagrangian Coherent 
Structures (LCSs) within the pulsating rip current surface velocity field. LCSs represent 
transport barriers hidden in the flow that control the motion of passively advected 
particles in unsteady flows. Numerical drifters released in the model demonstrated 
narrow streaks of floating material observed outside the surf zone on rip-channeled 
beaches. The model results showed that the LCSs were elongated in the cross-shore at the 
locations of rip currents, allowing the transport of drifters offshore until they reached the 
offshore extent of the LCS. Then they could move only alongshore and were shown to 
gradually re-enter the surf zone. Occasionally, the LCSs became detached from the 
nearshore rip current circulation cell at the outer surf zone due to VLF eddy motions. 
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C. FIELD EXPERIMENT 

1. Field Site 

The exchange of material between the surf zone and the inner shelf was examined 
during an extensive field experiment on a rip-channeled beach in southern Monterey Bay, 
California, in May 2009. Bathymetry, offshore waves, wind, tidal elevation and currents 
were measured throughout the field experiment. At the field site in Sand City, CA, the 
beach is composed of a relatively steep foreshore (1/10 slope) with straight and parallel 
contours, flattening out inshore (1/100 slope), with quasi-periodic, 0(125m), incised rip 
channels and continuing with a 1/20 offshore slope and straight and parallel contours 
seaward of the breaker zone (Figure 19). Feeder channels were present near the shoreline 
and converged with the incised rip channels. Regular bathymetric surveys were 
conducted throughout the experiment using a kinematic global positioning system (GPS) 
device mounted to a sonar-equipped personal watercraft (MacMahan 2001) to measure 
the inner shelf and outer surf zone, and a walking person to measure the shallow surf 
zone and dry beach. The local coordinate system used in this work is defined as positive 
offshore and to the south. 

2, Wave, Wind, Tide and Current Measurements 

Offshore waves were measured in 13 m water depth with an acoustic Doppler 
current profiler (ADCP) (Figure 20). There were gaps in the offshore ADCP wave data 
due to instrument fouling, therefore during those times wave parameters measured by 
National Data Buoy Center (NDBC) buoy 46240, located in 18.5 m water depth and just 
south of the field site, were shoaled to 13 m using linear wave theory (Dean and 
Dalrymple 1984) and hourly-averaged (Figure 20, gray lines). Wind speed and direction 
were recorded atop a 10 m tower onshore of the field site and rotated to the local 
coordinate system, and wind stress was calculated using the method by Large and Pond 
(1981) (Figure 20). Tidal elevation was obtained from the nearby National Ocean Service 
(NOS) tidal station 9413450 (Figure 20). 

Five ADCPs were deployed to obtain in-situ Eulerian current and pressure 


measurements for 17 days in an alongshore array at the approximate offshore edge of the 
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surf zone (Figure 19). ADCPs 1 to 4 had a sample period of = 2 s with a measurement 
bin size of 0.11 m, and ADCP 5 sampled at = 3 s with a measurement bin size of 0.25 
m. For the last 4 days of the experiment (yeardays 134 to 138), one ADCP was moved 
offshore (labeled ADCP 6), in line with a rip ehannel in the surf zone, to measure inner 
shelf flows, and three ADCPs remained in the alongshore array (Figure 19). During this 
time, ADCPs 1, 3 and 6 had a sample period of = 1 s with a measurement bin size of 
0.10 m, and ADCP5 sampled at = 3 s with a measurement bin size of 0.25 m. The 
hourly-mean, depth-averaged eross-shore and alongshore veloeities measured by the 
ADCP in the primary rip channel (ADCP 3: x = 140 m, y = 100 m) are shown in Figure 
20, and the depth-averaged cross-shore velocity averaged over the duration of each drifter 
deployment is given in Table 1. 

In order to estimate the surfzone-width. Lx, the mean wave conditions in 13 m 
water depth (measured by the offshore ADCP when data were available or by NDBC 
buoy 46240) were shoaled to the shoreline using a mean cross-shore beach profile, 
adjusted for the tidal elevation, with a depth-limited wave-breaking criterion of Hsih > 0.4 
(Thornton and Guza 1982). The cross-shore location of the surfzone boundary, Xsz, is 
defined as the cross-shore location of the maximum significant wave height Hs, and the 
cross-shore location of the shoreline, Xsh, is defined as the cross-shore location where h = 
0 m, with the surfzone-width Lx = Xsz - Xsh- Throughout this work distances are given in 
terms of the number of surfzone-widths from the shoreline, where one Lx represents the 
edge of the surf zone. 

3, Drifter Deployments 

A fleet of Lagrangian surfzone drifters equipped with GPS devices was deployed 
to measure surface current patterns outside the surf zone. See MacMahan et al. (2009) for 
a complete description of the drifter design and performance. The drifter positions, 
sampled every 2 s, were converted to the local coordinate system and filtered to remove 
erroneous points. Velocity estimates were computed using a forward-differencing 
scheme. Nine drifter deployments were conducted during varying wave and tidal 
conditions, each lasting approximately 3 hours with between 28 to 45 drifters being 
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released. In an attempt to eapture the spatial variation of rip-eurrent flows that exit the 
surf zone, drifters were released using three approaehes: 1) by swimmers in a eluster at 
the offshore edge of a rip ehannel, just outside of the surf zone, 2) by a boat in an 
alongshore line just outside of the surf zone, spanning four to five rip-ehannel/shoal 
systems, and 3) by a boat in a eross-shore line spanning from the edge of the surf zone 
and extending up to two surfzone widths offshore, in line with a rip ehannel (Figure 19). 
Conditions during the drifter deployments eonsisted of offshore waves with signifieant 
wave heights of 0.43 to 1.31 m, peak wave periods Tp of 6.7 to 12.6 s, total wind stress 
magnitude Ir^l < 0.04 N/m^, and Lx of 78 to 123 m (Table 1). All drifter deployments 
oeeurred during low tides (tidal elevation < MSL), to ensure rip eurrent flows would be 
indueed by the rip-ehanneled morphology, and during weak wind eonditions, defined as 
|t^|< 0.05 N/m^ by Ohlmann et al. (2012), allowing for the wave-driven transport to be 
examined separately. 

D, LAGRANGIAN OBSERVATIONS 

1. Qualitative Drifter Patterns 

In general, the drifters released outside the surf zone initially moved seaward due 
to rip eurrents and eventually returned shoreward in an areing pattern, at times re-entering 
the surf zone over shoals, with no drifters being permanently removed from the nearshore 
region. The qualitative drifter observations support the numerieal modeling results of 
Reniers et al. (2010), whieh showed LCSs in the rip eurrent flows had offshore limits 
prohibiting farther eross-shore transport and resulting in the shoreward return of drifters. 
Results from eaeh individual drifter deployment are shown in Figures 21 and 22, where 
the traek of every drifter that was released during the deployment is plotted, and the eolor 
of the drifter traek represents the speed of the drifter. The drifters typieally moved one to 
two Lx seaward of the surf zone before returning shoreward, with a maximum eross-shore 
extent of about six Lx. Two distinguishable drifter patterns were observed: 1) the drifters 
moved seaward and returned sharply baek shoreward a short alongshore distanee from 
where they exited, resulting in eross-shore exehange that was loeally eontained (Figure 
21), and 2) the drifters moved seaward and travelled farther in the alongshore direetion as 
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they gradually moved shoreward, resulting in cross-shore and alongshore exchange 
(Figure 22). During two of the drifter deployments (yeardays 135 and 136), the majority 
of the drifters moved into and remained inside the surf zone and no distinguishable 
offshore pattern was observed (not shown), therefore the results from these days will not 
be further discussed qualitatively. 

a. Locally Contained Cross-Shore Exchange 

Locally contained cross-shore exchange pattern occurred on yeardays 122, 129, 
and 132, when the rip currents were stronger during low tides and the alongshore currents 
outside the surf zone were weak (indicated by small mean alongshore drifter velocities 
outside the surf zone, as shown in Table 1). On yearday 122 (Figure 21), drifters released 
at the edge of the surf zone, offshore of a rip-channel, tended to move seaward and to the 
south, moving back shoreward and re-entering the surf zone over shoals less than 350 m 
down the beach. A cluster of drifters released about 2.5 Lx offshore in line with the rip 
channel were also pushed slightly seaward and south by the rip current pulse, then 
returned shoreward. The maximum extent of all drifters, whether released at the edge of 
the surf zone or farther offshore, was approximately the same and was 3.6 Lx. At different 
times during the deployment, drifters were observed to move shoreward and seaward at 
approximately the same location offshore of a rip channel (y = 250 m, x = 200 m), 
depending on if the rip current was pulsing or not at that time. 

On yearday 129 (Figure 21), the first drifter release occurred just after low tide, 
with a total of 27 drifters being released in an alongshore line approximately 1.5 Lx 
offshore, spanning the width of two rip channels and two shoals. In general, drifters that 
were released to the south of the primary rip channel moved directly shoreward and 
entered the surf zone over a shoal. Drifters that were released offshore of the primary rip 
channel and to the north initially moved seaward, then turned either north or south and 
moved back shoreward. The drifters that moved north travelled about 200 m farther 
seaward, reaching a maximum distance of 3.8 Lx offshore, before arcing sharply and 
moving back shoreward. These drifters returned shoreward within 200 m or less in the 
alongshore direction from where they initially moved seaward. Two hours after the initial 
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drifter release, another alongshore line of drifters was released at the approximate same 
cross-shore location spanning two shoals and a rip channel. All of these drifters moved 
shoreward, most likely due to the rising tide and the weakening of the rip current. 

Similar results were seen on yearday 132 (Figure 21), where drifters released 
offshore of shoals moved directly shoreward, entering the surfzone over the shoal. 
Drifters released offshore of the primary rip channel moved seaward and to the north, and 
returned shoreward within 300 m in the alongshore of their release location. Interestingly, 
all of these drifters were pushed back offshore when they reached y = -100 to -200 m 
(offshore of a rip channel), due to a neighboring pulsing rip current. The drifters again 
moved seaward and north, then returned shoreward within 200 m or less in the 
alongshore. The maximum offshore extent of the drifters on this day was 3.7 Lx. 

b. Cross-Shore and Alongshore Exchange 

Cross-shore and alongshore exchange pattern occurred on yeardays 125, 126, 130, 
and 137, when the rip currents were weak during high tides and the alongshore current 
outside the surf zone was stronger (Table 1). On yearday 125 (Figure 22), the first drifter 
release consisted of 29 drifters being released in an alongshore line, approximately 20 m 
apart, spanning an alongshore distance of about 600 m and five shoals and four rip 
channels. The surf zone was relatively wide {Lx = 122 m) this day, so the alongshore line 
was actually inside of the surf zone. Drifters released to the south of the primary rip 
channel (y > 100 m) meandered alongshore to the north inside the surf zone, following 
the underlying bathymetry contours. A number of these drifters ended up on the beach, 
and others were pushed seaward due to the rip current aty = 100 m. The drifters released 
near the primary rip channel and to the north, moved seaward and to the north. Once the 
drifters were outside of the surf zone they continued to move seaward and to the north, 
and then gradually moved back shoreward, with some even being pushed back offshore 
due to rip currents to the north. These drifters were picked up before they could re-enter 
the surf zone, which was 500 to 700 m alongshore to the north from the primary rip 
current where they exited the surf zone. The cross-shore extent of the drifters on this day 
was within 2 Lx. Additional drifter releases on this day consisted of clusters of drifters 
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being released in the primary rip channel near the offshore edge of the surf zone, and 
these all moved in a similar pattern seaward and to the north. A similar drifter pattern was 
observed on yearday 126 (Figure 22). Again, an alongshore line of drifters was released, 
this time just outside of the surf zone, and the drifters either moved into the surf zone or 
alongshore and were eventually pushed seaward and to the north. The cross-shore extent 
was about 2.8 Lx, and again the drifters were picked-up greater than 600 m to the north 
from where they were pushed seaward as they were beginning to gradually move 
shoreward. 

On yearday 130 (Figure 22), in general, drifters that were released outside of the 
surf zone and offshore of the primary rip channel or to the north, were pushed seaward as 
they moved alongshore to the south, and then began moving shoreward once they were 
about 300 m to the south of the rip channel (y = 400 m). Drifters that were released 
outside of the surf zone and to the south of the primary rip channel moved predominantly 
alongshore with little to no seaward movement and gradually shoreward. The cross-shore 
extent of the rip current was demonstrated by a release of drifters in a cross-shore array 
offshore of the primary rip channel, approximately 1.5 to 2.5 Lx offshore, where all of the 
drifters moved south and seaward and seemed to converge at a cross-shore extent of 
about 3.3 Lx, where they then began moving shoreward. Drifters that were released at the 
surfzone boundary and inshore were immediately pulled into the surf zone over a shoal. 
As the tide rose during the deployment, the rip current weakened and no drifters moved 
seaward, and the drifters inside of the surf zone exhibited a meandering current pattern. 
Throughout the deployment, all drifters that moved seaward and alongshore outside of 
the surf zone were picked-up 300 to 500 m south of the primary rip channel, before they 
had a chance to re-enter the surf zone, indicating a significant amount of alongshore 
exchange was occurring. 

The drifter deployment on yearday 137 (Figure 22) occurred during low tide with 
relatively small wave conditions, resulting in a narrow surf zone {Lx = 81 m). Drifters 
were released in a cross-shore array in line with the primary rip channel between 1 to 2 
Lx. The rip current pulse pushed the drifters seaward and to the south, with a maximum 
cross-shore extent of about 2.5 Lx. All of the drifters gradually returned shoreward, and 
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approached the surf zone boundary greater than 500 m from the primary rip channel 
where they were released, with a majority of the drifters then being pushed baek offshore 
and to the south exhibiting a meandering alongshore eurrent pattern outside of the surf 
zone within 2 Lx. Shear in the alongshore current outside of the surf zone ean be seen in 
the drifter results, with the alongshore current magnitude decreasing near the surfzone 
boundary. 


2. Drifter One-Particle Diffusivity Statistics on the Inner Shelf 


The movements of the drifters outside the surf zone were evaluated using one- 
particle statistics, which describe the average drifter evolution from a common release 
point. The trajectories of the drifters are represented as the collection of relative position 
displacements, r{f) = x{t'o + f) - x{t'o), for a relative time step t', computed as the 
difference between the drifter position, x = (x, y), at an earlier time t'o from the position a 
time t' later. Each recorded drifter position is considered a possible starting position, 
x{t'o), therefore a drifter time series with N total position measurements yields {N - 1) 
values for r{t' = 2 s), where the drifter positions were sampled at dt = 2 s. Assuming a 
spatially homogeneous and temporally stationary flow field, all drifters released during a 
given deployment can be viewed as originating from a common release point, and the 
random motions represented by the drifter trajectories can be described in the same 
manner as the bulk molecular diffusion (Taylor 1922). The one-particle absolute 
diffusivity, Ky, is related to the change in drifter displacement variance, a^y, with time 




where ij = x,y and y is the second moment of displacements 


(V) 


^5 " JJ r'r'pU/^d'yr'dr' 


( 8 ) 

and P{rx', ry'; t') is the probability distribution function (pdf) of the anomalous relative 
particle displacements r'={rx', ry). The anomalous relative position displacements were 
calculated for each drifter as r^t) = r{t') - Rit), where /f(t') is the mean spatial 
displacement of all drifters for the time step t' and all arbitrary t'o (Spydell et al. 2007; 
Brown et al. 2009). The pdf of the anomalous relative position displacements P{rx', ry'-, t') 
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represents the ensemble average evolution of the drifters from a eommon initial point 
souree. 

In order to evaluate the fate of material outside the surf zone, only drifter 
positions reeorded outside the surf zone (x > Lx) for durations of at least one minute were 
evaluated. The probability distributions of the anomalous displaeements for t' = 2 s, P{rx', 
ry\ t' =2 s), were relatively eentered about zero in the eross-shore direetion, indieating an 
equal amount of drifters moved seaward and shoreward (Figure 23). Examining just the 
eross-shore anomalous relative displaeements, the distributions were similar for all drifter 
deployments regardless of drifter pattern observed (Figure 23), and were slightly 
negatively skewed (skewness of -0.24 to -0.05), and had kurtosis values ranging from 2.8 
to 4.9, suggesting that the distributions were basieally Gaussian (kurtosis = 3). Therefore, 
there was a relatively equal ehanee that a drifter outside of the surf zone would move 
seaward or shoreward, whieh supports the qualitative observations that the drifters did 
not leave the nearshore system (Figure 21 and 22). 

The eross-shore absolute diffusivities Kxx (Equation 7) outside the surf zone 
showed similar behavior for all eurrent patterns, with Kxx rising to a peak then deereasing 
to an asymptotie value (Figure 24). The time to reaeh the peak Kxx was generally 15 to 19 
minutes, with the exeeption of yearday 129 whieh oeeurred at 36 minutes. The values of 
Kxx were greater during observations of loeally-eontained eross-shore exehange, with peak 
values of of 1.11 to 2.41 m^ s"' and asymptotie values of k^xx of 0.36 to 1.74 m^ s'\ 
During observations of alongshore and eross-shore exehange, the differenee in the peak 
and asymptotie values of Kxx were relatively small, with peak values of Kxx of 0.28 to 1.00 
m^ s'^ and asymptotie values of k'^xx of 0.27 to 0.73 m^ s"'. The alongshore absolute 
diffusivities Kyy outside the surf zone inereased approximately monotonieally to an 
asymptotie level for all drifter releases (Figure 24). 

During observations of loeally-eontained eross-shore exehange, k'^xx was greater 

than K^^yy (1.25 to 1.83) due to the predominant movement and spreading of the drifters in 

the eross-shore direetion while remaining within a relatively narrow alongshore region. 

However, during observations of eross-shore and alongshore exehange, values of K^yy 

(0.36 to 5.01 m^ s'^) were generally greater than k'^xx due to the drifters spreading more in 
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the alongshore direetion and remaining within a smaller eross-shore extent. For these 
reasons, k^xx was largest when there was more eross-shore exehange, and K^yy was largest 
when there was more alongshore exchange (Figure 24). 

The asymptotic values of diffusivity obtained outside the surf zone {k^xx = 0.27 to 
1.74 m^ s'\ K^yy = 0.36 to 5.01 m^ s"') were smaller than those measured in previous 
experiments inside the surf zone {k^xx = 0.10 to 2.50 m^ s"', K^yy = 2.00 to 12.60 m^ s'^) 
(Spydell et al. 2007; Brown et al. 2009), regardless of the measured current pattern (i.e. 
rip current circulation versus alongshore current). Additionally, Kxx and Kyy rise more 
rapidly inside the surf zone, and the time to approach k^xx and K^yy is longer outside the 
surf zone, similar to the limited inner shelf observations of drifter diffusivity obtained by 
Spydell et al. (2007) and Brown et al. (2009). For all drifter deployments, the absolute 
diffusivities lacked oscillations, indicating a lack of rotational structures outside the surf 
zone (Brown et al. 2009), and that different eddy processes with different spatial and time 
scales are responsible for mixing material on the inner shelf. 

3, Cross-Shore Variation in Drifter Cross-Shore Velocity Magnitude 

The behavior of the drifters moving seaward {udrift.off) and shoreward {udrift.on) was 
evaluated by examining the drifter cross-shore velocity magnitude a 

function of cross-shore location Lx. The study area was divided into 25 m bins in the 
cross-shore direction, and the seaward {udrift > 0 m/s) and shoreward {udrift < 0 m/s) 
moving drifters in each cross-shore bin were averaged separately over the entire 
alongshore distance (Figure 25). Outside the surf zone, the magnitude of the seaward and 

shoreward moving drifters were relatively equal, |Mrfri/i,on| ~ |“rfnyr,o#| > indicating the total 

alongshore-averaged cross-shore drifter velocity was approximately zero, as a function of 
cross-shore distance. At the edge of the surf zone, the shoreward drifter velocities were 
greater than the seaward component, possibly due to the surfing of the shoreward moving 
drifters in the breaking waves. Inside the surf zone, the shoreward drifter velocities over 
the shoals were larger than the seaward drifter velocities in the rip channel, due to 
continuity and differences in water depth (MacMahan et al. 2004a). The were 
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similar for all drifter deployments, that is the drifter veloeities measured for the loeally 
eontained eross-shore exchange drifter pattern were not necessarily greater or less than 
those measured for the cross-shore and alongshore exchange drifter pattern. 

The drifter cross-shore velocity magnitude tended to decrease over the 

distance of one surfzone-width offshore of the surf zone (i.e., from x = ILx to x = 2Lx), 
and was relatively constant offshore of x > 2Lx, for all drifter deployments 

(Figure 25). The average change in offshore of the surf zone was 

= ~ 0.12 iTi/s, mcludlng scaward and 

shoreward drifter velocities. The drifters exiting the surf zone due to rip current pulses 
and moving seaward had large velocities at the edge of the surf zone, which subsequently 
decreased as they continued to move seaward, resulting in a deceleration of the drifters as 
they exited the surf zone. The drifters outside the surf zone moving shoreward 
demonstrated an increase in velocity magnitude as they approached the edge of the surf 
zone. Therefore the shoreward moving drifters were actually accelerating as they re¬ 
entered the surf zone over shoals. These results differ from those of Ohlmann et al. 
(2012), who observed a tendency for drifters released on the inner shelf to decelerate as 
they moved shoreward, due to a net offshore Lagrangian surface velocity, on beaches that 
do not typically support rip currents. 

The measurements outside the surf zone were compared to theoretical 

estimates of the Stokes drift, ust- The cross-shore variation of ust{z) was computed using 
wave parameters measured in 13 m water depth and shoaled to the shoreline using the 
mean cross-shore beach profde adjusted for tides. The ust{z) was averaged over the upper 
0.5 m of the water column corresponding to the depth of the surface drifters, indicated as 
Ma(0.5 m). In all cases, was greater than the Ma(0.5 m) estimates as a function 

of cross-shore location (Figure 25). Therefore, the return of the drifters shoreward was 
not solely due to the net onshore Stokes drift due to waves. 


44 



E, EULERIAN OBSERVATIONS 

1. 3-D Variations of the Rip Current Flow 

The Eulerian measurements made by the array of ADCPs show differenees in the 
magnitude and vertical structure of the cross-shore flows between the ADCPs on the 
shoals (ADCPs 1 and 5) and the ADCPs in the rip channel (ADCPs 2, 3, 4, and 6). To 
examine these differences, the cross-shore velocities measured during low tides, when rip 
current flows are known to be strongest (MacMahan et al. 2006), were evaluated. For 
each ADCP, the vertical current profiles were interpolated to a depth-relative reference 
frame zlh, where zlh = 0 is at the sea-bed and zlh = 1 is at MSL, and a 3-hour mean 
centered about low tide was computed. These mean current profdes were then sorted by 
local Hsih and bin-averaged, as shown in Figure 26. 

In general, at a given ADCP location, the magnitude of the cross-shore velocity 
was greater for larger HJh, representative of larger waves and/or lower tides. These 
Eulerian measurements support the Fagrangian drifter observations of decreasing 

\udrift,ojfio\ with increasing distance offshore, where smaller HJh are indicative of 

measurements farther outside of wave breaking. Additionally, for a given HJh range, the 
magnitude of the seaward cross-shore velocities were greatest in the center of the rip 
channel near the edge of the surf zone at ADCP 3, and were smallest over the shoals at 
ADCPs 1 and 5 (Figure 26). 

The measured Eulerian cross-shore velocity profiles revealed more shear over the 
vertical, particularly near the surface, within the rip channel, and depth-uniform currents 
below the wave trough over the shoals. In the center of the rip channel near the surfzone 
boundary at ADCP 3, the measured mean rip current velocity profile was surface 
dominated and had more vertical structure for smaller HJh, and was more depth-uniform 
for larger HJh representative of times when it was near the onset of wave-breaking and 
just outside the surf zone (Figure 26). Well offshore of the surf zone, at ADCP 6, the 
measured seaward rip current flow was also surface dominated, demonstrating the cross¬ 
shore extent of the rip current flows that exit the surf zone. These results support the 
laboratory measurements of Haas and Svendsen (2002), and reveal that the vertical 
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profile of the rip eurrent veloeity is strongest and depth-uniform as it exits the surf zone, 
then beeomes surfaee dominated and decays in magnitude with distance offshore. The 
cross-shore velocity profiles measured by ADCPs 1 and 5 over the shoals were vertically 
uniform below the wave trough for all values of Hsih, and were even onshore at times for 
Hsih > 0.4 when they were inside the surf zone. 


2. Volumetric Transport Estimates 


The volumetric transport on the rip-channeled beach is evaluated using the 
Eulerian current measurements, and is compared to conventional theory. If the 
conventional theory of rip current transport (Figure 18e) is valid, a volumetric balance 
across the surfzone boundary on a rip-channeled beach would exist where the continuous 
deposition of water into the surf zone via the Stokes transport over the shoals {Qst Kip) 
would be balanced by the steady compensating offshore Eulerian rip current transport 
through the rip-channels {QE.rip brip). Therefore, the depth-averaged offshore Eulerian rip 
current flow would equal 


^s,Kip 

■ =- - 

E,rip 7 

^rip 


( 9 ) 


according to Equation (6), where the depth-averaged Stokes drift associated with the 
surface gravity waves is 
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Observed hourly-mean cross-shore velocity profiles measured by each ADCP 
were used to estimate the depth-averaged Eulerian flows below the wave trough, . The 
velocity was assumed uniform between the lowest measurement (generally 0.5 m above 
the sea-bed) and the sea-bed (Eentz et al. 2008), and was estimated using trapezoidal 


integration and averaging over the depth of the wave trough. The measured was 
normalized by {\6c/gh), represented as ^, and were bin-averaged as a function of HJh, 
and are shown in Figure 27. Theoretical estimates of were computed using Equation 
(9), where the rip-channel spacing =125 m and width trip = 40 m were based on the 
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measured bathymetry, and also normalized by {\6clgh), represented as are 

shown in Figure 27 (dashed blaek line). Additionally, the theoretieal normalized 

based on Equation (10), as a funetion of HJh, is shown in Figure 27 (solid blaek line), 
whieh represents the balanee between u^z) and Ma(z) found on the inner shelf by Lentz et 
al. (2008). The loeation of the onset of wave breaking is defined as HJh > 0.4 and is 
shown in Figure 27 to deseribe the observations as being inside or outside of the surf 
zone. 

In general, field observations of ^ inereased with inereasing HJh, assoeiated 

with inereasing wave heights and/or lower tides. For small waves and/or higher tides 
{HJh < 0.2), measured was between the two theoretieal volumetrie balanees 

(Equations (9) and (10)), and was generally equal to (Figure 27, solid black line). 

This is consistent with observations on the inner shelf by Lentz et al. (2008), and 
indicates zero net Lagrangian cross-shore transport well outside of the surf zone. A 
transition in exists approaching the surfzone boundary (0.2 < HJh < 0.4), where for 

a given range of HJh, the magnitude of ^ was dependent upon the alongshore location 
of the ADCP (i.e., within the rip channel). The measured by ADCP 1 (Figure 27, 

blue) and ADCP 5 (Figure 27, yellow) remained consistent with Equation (10) and Lentz 
et al. (2008). This is attributed to ADCPs 1 and 5 being located at the borderline of the rip 
current width, yet not over the shoals, and therefore the measured velocities were 
typically small and not representative of strong rip current nor onshore flows. However, 
£ measured within the rip channel by ADCPs 2, 3, and 4 (Figure 27, green, magenta, 

and cyan) tended to exceed the theoretical Stokes drift transport (Equation (10)) and 
the conventional rip current volumetric transport (Equation (9)). The fact that the 

measured ^ in the rip current was greater than ^ near the edge of the surf zone 

indicates an imbalance in the volume transport across the surfzone boundary on the rip- 
channeled beach, with more offshore transport exiting the surf zone through the rip 
channels than the onshore Stokes transport entering the surf zone over the shoals. 
Additionally, there is a greater imbalance (larger offshore transport) at the surfzone 
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boundary that decreases with distance offshore (decreasing HJh), consistent with the 
drifter observations. The consistency in the measurements by ADCPs 2, 3, and 4 
demonstrates that each ADCP was measuring seaward rip current flows, and the 
alongshore spacing of the ADCPs verifies the narrow rip current width, trip = 40 m, used 
in Equation (9). Inside of the surf zone {HJh > 0.4), measured on the shoals was less 

than the theoretical ^ (Equation (9)), but unfortunately no measurements of ^ in 
the rip currents were obtained shoreward of wave breaking. 

F. DISCUSSION 

1, Mass Balance Controlling Cross-Shore Exchange on a Rip-Channeled 
Beach 

The extensive data set presented here, incorporating measurements of Eulerian rip 
current flows and transport inside and outside the surf zone, as well as Eagrangian drifter 
observations on the inner shelf, combined with previous surf zone drifter observations 
(MacMahan et al. 2010a), is contrary to the conventional paradigm of rip current 
transport. Based on the conventional theory (Equation (6)), if onshore Stokes transport 
was continuously depositing water into the surf zone over the shoals {Qst 2„p), a steady 
compensating offshore transport through the rip channels would exist {QE.rip brip). 
However, the measured ^ in the rip channel was greater than suggested by 

theory (Eigure 27), indicating a greater amount of offshore transport by the rip current 
across the surfzone boundary than the filling of the surfzone by Stokes transport over the 

shoals. Also, the drifter observations show that ^ wa(0.5 m), indicating that 

Stokes transport associated with short waves does not seem to transport mass shoreward 
in the nearshore and across the surfzone boundary on the rip-channeled beach. The excess 
£ measured above does not indicate a net offshore transport, however, because 

the drifter observations on the inner shelf show that no drifters completely left the 

nearshore system. The seaward and shoreward drifter velocity magnitudes were relatively 

equal outside the surf zone (Eigure 25), implying there was no net Eagrangian transport, 

consistent with Lentz et al. (2008). These results suggest that a recirculation exists across 

the surfzone boundary. Additionally, based on the high retention of drifters within the 
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surf zone observed on this same beaeh by MacMahan et al. (2010a), the offshore flow via 
the rip current ^ could not be constant. This is consistent with the findings of Reniers 

et al. (2009, 2010), who found rip current exchange is driven by VLF rip current 
pulsations forced by wave groups for a rip-channeled beach with 0(10 minute) 
timescales. Therefore, as material is transported offshore by rip current pulsations, some 
factor besides Stokes transport must be responsible for pulling material shoreward on the 
inner shelf and into the surf zone. This also explains why the drifters did not decelerate as 
they approached the surfzone boundary, as observed by Ohlmann et al. (2012) on the 
inner shelf where there were no rip channels or surfzone rip current circulations present. 

It is hypothesized that a mass balance exists on the rip-channeled beach, where 
the surfzone rip current circulations both eject material offshore in the rip channels as 
VLF pulses, as well as pull material back into the surfzone over shoals (Figure 18e). This 
is investigated by further examining the Lagrangian and Eulerian measurements during 
the time of the drifter deployment on yearday 132, when locally-contained cross-shore 
exchange was observed. The mean cross-shore velocity observations of the drifters from 
the entire deployment on this day were re-evaluated by spatially binning and averaging 
the observations into 10 m by 10 m bins, with the mean cross-shore movement of the 
drifters at each bin represented as a percentage of time that the drifters were moving 
seaward, i.e. always (100% of the time, red) or never (0% of the time, blue), as shown in 
Figure 28. It is evident that the drifters moved out of the surf zone and seaward through 
the rip channel, and returned shoreward and re-entered the surf zone over the shoals. 
Figure 28 also shows that for this particular drifter deployment, out of the ADCPs that 
were within the rip channel, only ADCPs 2 and 3 were located where drifters were 
exiting the surf zone, with ADCPs 4 and 5 being located where drifters were moving both 
seaward and shoreward, and ADCP 1 being located where a slightly higher percentage of 
drifters were moving shoreward. This demonstrates the narrow width of the rip current, 
and highlights the asymmetry in the re-circulation of the flow. Therefore, the values of 
£ shown in Figure 27 measured by the ADCPs in the rip channel (ADCPs 2, 3, and 4) 

are assumed to describe the seaward rip current transport, and support the observations of 
the seaward moving drifters. However, ADCPs 1 and 5 were on the edge of the rip 
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current flow and did not necessarily measured the strong onshore flows occurring over 
the shoals, which is why the values of ^ measured by ADCPs 1 and 5 do not match 

the observations of the shoreward moving drifters. 

The vertical structure of the cross-shore flows measured by the ADCPs during the 
drifter deployment on yearday 132 (excluding ADCP 2, which was inoperative this day) 
is shown in Figure 29 (top row). The profdes show offshore flows, which were surface 
dominated, measured inside the rip channel (ADCPs 3 and 4), and onshore flows, which 
were vertically uniform below the wave trough, measured by ADCPs 1 and 5, where the 
measured onshore flows were a fraction of the measured offshore rip current flows. The 
presence of VLF pulses in the measured rip current surface flow was examined by low- 
pass filtering if < 0.004 Hz) the near-surface cross-shore velocity measured by the 
ADCPs (Figure 29, bottom row). Oscillations were present in the low-passed fdtered 
offshore surface rip current velocity (ADCP3 and 4) at the VLF timescale with relatively 
large amplitudes (0.5 to 1.0 m s'*) that represent the VLF rip current pulsations 
responsible for the rip current exits. Additionally, similar VLF oscillations of smaller 
magnitude (0.1 to 0.2 m s'*) were present in the onshore surface velocities over the 
shoals. These results demonstrate the VLF pulsating nature of the seaward rip current 
flows as well as the shoreward flows re-entering the surf zone. 

A new paradigm of rip current behavior is presented (Figure 18e). It is proposed 
that the onshore Stokes transport associated with short waves does not transport mass into 
the surf zone on a rip-channeled beach, but rather deposits momentum at the wave group 
timescale that drives the cross-shore exchange and forces episodic, VLF rip current 
pulsations. To balance the mass transported offshore by the rip current exits, the surfzone 
rip current circulations over the shoals pull material back into the surf zone at the same 
VLF time scale. Ultimately little or no material is completely removed from the 
nearshore system. 

G. SUMMARY AND CONCLUSIONS 

Cross-shore exchange on a rip-channeled beach at Monterey, California, was 
examined using an extensive set of Lagrangian and Eulerian field measurements of surf 
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zone and inner shelf rip eurrent flows. The results revealed differenees in eross-shore 
transport from that measured on the inner shelf of alongshore homogeneous beaehes, and 
is eontrary to the eonventional theory of mass transport due to rip eurrents. Position- 
traeking drifters released on the inner shelf were observed to move seaward due to rip 
eurrent pulses, then return shoreward in an areing pattern, re-entering the surf zone over 
shoals. Two drifter patterns were observed: 1) loeally-eontained eross-shore exehange, 
where the drifters moved seaward and returned sharply baek shoreward a short 
alongshore distanee from where they exited (Figure 21), and 2) eross-shore and 
alongshore exehange, where the drifters moved seaward and traveled farther in the 
alongshore direetion as they gradually moved shoreward (Figure 22). The eross-shore 
extent of the drifters was typieally one to two surfzone-widths beyond the surfzone 
boundary. The qualitative drifter observations were similar to the numerieal model results 
of Reniers et al. (2010), whieh showed that at the loeations of rip eurrents LCSs 
eonstrained the offshore extent of drifter transport in the eross-shore where they then 
moved alongshore and re-entered the surf zone. Additional LCS model simulations with 
drifters seeded on the inner shelf would help to further explain the fate of the detaehed 
LCSs and eross-shore exehange of material on the inner shelf 

On the inner shelf of the rip-ehanneled beaeh, the rate of spreading and mixing 
{k^xx = 0.27 to 1.74 m^ s'\ K^yy = 0.36 to 5.01 m^ s’') was slower than that measured 
inside the surf zone during previous studies (Brown et al. 2009), suggesting there is less 
shear in the flow field outside the surf zone. The probability distributions of the drifter 
eross-shore anomalous relative displaeements outside the surf zone were basically 
Gaussian, indicating there was a relatively equal chance that a drifter anywhere on the 
inner shelf would move seaward or shoreward, consistent with the qualitative 
observations that no drifters left the nearshore system. The seaward and shoreward 
moving cross-shore drifter velocities were approximately equal in magnitude as a 
function of cross-shore location, resulting in a total Lagrangian velocity of approximately 
zero The seaward moving drifters were observed to decelerate after exiting the surf zone 
through rip channels, while the shoreward moving drifters were observed to accelerate as 
they re-entered the surf zone over shoals. Near the surfzone boundary, the depth-averaged 
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Eulerian flows measured in the rip ehannel were largest near the edge of the surf 

zone, eonsistent with the drifter observations, and were greater than the offshore rip 
eurrent transport suggested by eonventional theory to balanee the filling of the 

surfzone by Stokes transport over the shoals. The magnitude of the shoreward moving 
drifters outside the surf zone was greater than the near-surfaee theoretieal Stokes drift, 
therefore Stokes drift was not responsible for moving the drifters shoreward and baek 
into the surf zone, and Stokes transport does not provide a eontinuous, wave mass flux 
into the surf zone aeross the surfzone boundary that would induee a eontinuous, 
eompensating mass flux out of the surf zone via rip currents on a rip-channeled beach, as 
conventionally thought. The results support a new theory of how material is exchanged 
and conserved within the nearshore region on a rip-channeled beach: 1) short waves 
deposit momentum into the surf zone at the wave group timescale, then 2) spatial 
variations in the wave groups force VLF motions and rip current circulations inside the 
surf zone that episodically expel material offshore through the rip channels, and 3) the 
surfzone rip current circulations act to pull the material back onshore into the surf zone 
over shoals to satisfy a mass balance, with ultimately little or no material being 
completely removed from the nearshore system on the rip-channeled beach. 
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IV. MIXING AND TRANSPORT ON A STEEP BEACH 


A. INTRODUCTION 

The mixing and transport of material in the nearshore, ineluding suspended 
sediment, pollutants, biologieal matter, and nutrients, influences beach composition, 
water quality, and ecosystem stability. There have been few field studies on mass 
transport and cross-shore exchange across the surfzone boundary, and most have 
occurred on wide, dissipative beaches in Southern California (Inman et al. 1971; Boehm 
2003; Grant et al. 2005; Clark et al. 2010). These beaches support communities with 
large human populations, and an understanding of the mixing and transport of pollutants 
on dissipative beaches has both health and economic relevance (Grant et al. 2005; Given 
et al. 2006; amongst others). In contrast, steep beaches typically do not support a lot of 
recreational activities, owing to their hazardous shore break. However, steep reflective 
beaches are important because they are common and are typically situated between rocky 
outcrops where there are often dense and diverse communities of organisms. Therefore, 
steep beaches can be studied to understand exchange processes that are important for 
those communities, such as biological recruitment. In general, a better understanding of 
how material is dispersed on beaches with varying topography and wave and current 
conditions is needed. 


Mixing and transport in the surf zone results from several physical processes with 
distinctive temporal and spatial scales, which can be separated into two dominant 
mechanisms that act to disperse material horizontally: 1) advection, which describes the 
bodily transport of material due to wave-driven surfzone currents, and 2) diffusion, which 
describes the spreading of material due to turbulent motions. These processes can be 
described by the two-dimensional (2-D), depth-averaged, advection-diffusion equation 
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where C is the depth-averaged tracer concentration, h is the total water depth, u and v are 
the cross-shore (x) and alongshore (y) components of the time- and depth-averaged 
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velocity, Kx and Ky are the eddy diffusion coefficients, and 5* is a source/sink term (Fischer 
et al. 1979). Solving Equation (11) for an instantaneous point source of dye mass Mq 
input at X = 0, y = 0, t = 0 (S' = Mo6{x,y,t)), with a no-flux boundary condition applied at 
the shoreline (x = 0) using the method of images (Fischer et al. 1979), and assuming 
constant diffusion coefficients, results in the solution (Rutherford 1994) 


C{x,y,t) = ^- - 1 — 

h 47ItJK^K 


2 exp 


[x-Ut) 


(y-Vt) 


where C{x,y,t) is in ppb, h = 0.5 m is the mean surfzone water depth, and U and V 
are the surfzone averaged cross-shore and alongshore currents. 


Previous studies used fluorescent dye as a tracer to examine surfzone mixing and 
transport on wide, dissipative beaches, focusing on either alongshore dispersion (Harris et 
al. 1963; Inman et al. 1971; Boehm 2003; Grant et al. 2005; Clarke et al. 2007) or cross¬ 
shore dispersion (Clark et al. 2010), but often lacked high resolution, 2-D measurements 
and utilized simple models to describe the 2-D evolution of the surfzone mixing and 
transport processes. Lagrangian surface drifters have also been used to study dispersion 
on alongshore-uniform (Spydell et al. 2007) and rip-channeled (Brown et al. 2009) 
dissipative beaches, however drifters do not perform well in the energetic, narrow surf 
zone of a steep beach. Visual observations of dye released in the surf zone reveal rapid 
vertical mixing (Clark et al. 2010), and initial 2-D horizontal (cross-shore and 
alongshore) dispersion until the surf zone was saturated (approximately uniform cross¬ 
shore concentration) (Harris et al. 1963; Inman et al. 1971; Clarke et al. 2007; Clark et al. 
2010). The dye was dispersed by advection and diffusion; however, once the width of the 
dye patch was of the order of the width of the surf zone (100 m to over 1 km alongshore 
from the dye release location), mixing by diffusion was not important and advection 
dispersed the dye, becoming one-dimensional (1-D) in the alongshore direction (Harris et 
al. 1963; Inman et al. 1971; Grant et al. 2005; Clarke et al. 2007; Clark et al. 2010). 


A measure of the amount of cross-shore exchange occurring on a beach can be 
expressed as the ratio of alongshore transport to cross-shore transport (Grant et al. 2005). 
In the absence of an ambient current, local wave conditions control the wave-driven 
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currents and influence the exchange ratio VjU. The case where V/U ~0 indicates 
complete cross-shore exchange, which occurs when waves break shore-normally and 
V jU increases with increased alongshore advection and decreased cross-shore exchange, 
which occurs when waves approach the shore obliquely. On a wide (50 m surf zone 
width), dissipative (1/50 beach slope), relatively alongshore uniform beach, that spanned 
0(km) in the alongshore, at Huntington Beach, California, Grant et al. (2005) found the 
cross-shore flux was 50 to 300 times smaller than the alongshore flux. The resulting ratio 
of V/U >50 indicates that alongshore advection dominates, and therefore suggests that 
most of the material remains within the surf zone and travels alongshore between the 
shoreline and active wave breaking with only a small amount of cross-shore exchange 
occurring continuously across the surfzone boundary. Field measurements of alongshore 
mass transport on dissipative beaches indicate an exponential decay of tracer as a 
function of distance along the shoreline from the release point (Inman et al. 1971; Boehm 
2003; Grant et al. 2005), consistent with a conceptual model where alongshore advection 
dominates both alongshore mixing and cross-shore mixing (Boehm 2003; Grant et al. 
2005). 

Previous studies estimated surfzone diffusivities by fitting tracer measurements to 
a Fickian diffusion model (Inman et al. 1971; Clarke et al. 2007; Clark et al. 2010). Using 
discrete tracer samples at the shoreline, and assuming constant alongshore currents and 
depth, and prescribing no-flux boundaries at the shoreline and wave-breaker line, Inman 
et al. (1971) and Clarke et al. (2007) found that in general Kx » Ky, however these 
estimates were not statistically reliable. Clark et al. (2010) estimated surfzone cross-shore 
diffusivities on a wide, dissipative beach at Huntington Beach, California, using 
ensemble-averaged cross-shore profiles of tracer concentration and a Fickian diffusion 
model with a no-flux boundary at the shoreline. They observed shoreline-attached cross¬ 
shore tracer concentration profiles, with Kx = 0.5 to 2.0 m^ s'\ and decreased tracer 
concentrations downstream of the source and seaward of the surf zone, consistent with 
lower diffusivity outside the surf zone (Harris et al. 1963; Inman et al. 1971; Clark et al. 
2010). Surfzone diffusivities on dissipative beaches have also been estimated using 
surface drifters, with asymptotic (long time) Kx = 0.7 to 1.5 m^ s'* and Ky = 2.0 to 4.5 m^ s' 
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' found for an alongshore uniform, dissipative beaeh (Spydell et al. 2007), and asymptotic 
Kx = 0.9 to 2.2 m^ s"' and Ky = 2.8 to 3.9 m^ s"' found on a rip-channeled beach (Brown et 
al. 2009). 

In this work, field observations of the mixing and transport processes that occur 
on a sandy steep, reflective beach are investigated and compared to the findings on 
dissipative beaches. Dye was released in the narrow surf zone and the temporal and 
spatial evolution was measured, resulting in simultaneous alongshore and cross-shore dye 
concentration profiles, and estimates of alongshore and cross-shore transport of dye and 
diffusion coefficients were obtained. 

B, FIELD EXPERIMENT 

1. Field Site, and Wave and Current Measurements 

Field observations were obtained during an experiment in June and July 2011 on a 
sandy, steep, reflective beach at Carmel River State Beach (CRSB), California. 
Bathymetry was surveyed multiple times throughout the experiment using a kinematic 
global positioning system (GPS) mounted to an electric kayak equipped with an 
echosounder to measure the offshore profile and to a backpack carried by a walking 
person to measure the dry beach, and the beach step was measured by a person swimming 
cross-shore transects with a survey rod affixed with a GPS atop at numerous alongshore 
locations. The survey data were combined and interpolated to a spatial grid, resulting in a 
single bathymetry. The alongshore-averaged, cross-shore beach profile was composed of 
a 1/7.6 subaerial beach slope, a steep 1/3 subaqueous beach step, and a 1/19 subaqueous 
offshore profde (Figure 30). Due to the steepness of the beach, depth-limited wave 
breaking occurred very close to the shore, resulting in a narrow surf zone that was 
intermittent in width at the frequency of the sea-swell waves. The plunging/surging 
breakers were effectively the run-up associated with the swash zone that alternately wets 
and dries the beach. Therefore, the surf zone on this beach is considered to be more of a 
large swash zone rather than a typical surf zone. 

Offshore waves and currents (Figure 31) were measured with a cross-shore array 
of four bottom-mounted, upward-looking acoustic Doppler current profders (ADCPs) and 
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pressure sensors (Figure 30), sampled at 1 Hz, throughout the experiment. Additionally, a 
vertical array of six Electromagnetic current meters (EMCMs) spaced at 0.2 m intervals 
was deployed in 0.4 m water depth relative to mean sea level (MSE) and sampled at 16 
Hz for 5 days (yeardays 169 to 173) to observe the vertical structure of the cross-shore 
and alongshore currents associated with the steep beach surf zone. Hourly, depth- 
averaged currents (uavg, Vavg) were computed from the ADCPs and EMCM array 
measurements to evaluate the mean currents inside and outside the surf zone (Eigure 
31c,d). The tidal elevation with respect to MSE during the experiment was approximately 
± 1 m (Eigure 31e), so the EMCMs were located inside and outside the region of active 
wave breaking as a function of tidal elevation and wave energy. A local wave breaking 
criterion of Hnns/h > 0.42 (Thornton and Guza 1982) was used to determine when the 
EMCMs were inside the surf zone, which was typically during low tides. Throughout the 
experiment, waves approached the field site obliquely from the northwest, resulting in a 
predominant alongshore current to the south (Eigure 3Id). In this work, a local coordinate 
system is used where positive is offshore and to the south (Eigure 30). 

2, Dye Experiments 

Dye was released on five different days over the course of a 10-day window 
during varying wave and alongshore-current conditions (Eigure 3Id), resulting in 15 
successful dye releases. An initial volume of 200 mL of fluorescent Rhodamine dye (20% 
weight per volume) was released as a slug in the surf zone (Eigure 32a) at mid-depth in 
the shallow water at the shoreline. Rhodamine dye was chosen since its decay rate is 
minimal and it is relatively stable in ambient light (Smart and Eaidlaw 1977; Grant et al. 
2005). The dye sensors used to measure dye concentration were 1 Hz internal logging 
WET Tabs fluorometers with a 0 to 230 ppb range. Clark et al. (2009) found that the 
instantaneous 1 Hz concentration errors were '-I ppb in a bubbly, sediment suspended 
surf zone with medium-grained quartz sand and minimal fine sediments. Similar 
measurement errors are assumed at CRSB, which had coarse sand with minimal fines. 

The temporal and spatial distribution of the dye was evaluated with an alongshore 
array of eight dye sensors (Eigure 30) mounted on vertical poles jetted into the sea-bed on 
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the upper beaehfaee (Figure 32e-e) with an approximate alongshore spaeing of 25 m. The 
dye sensors were mounted vertically with their optical sensors faced downward 
approximately 10 cm off of the sea-bed, such that they were submerged for the longest 
duration in the surfzone waves. Additionally, during some of the dye releases, swimmers 
each equipped with a dye sensor and a Global Positioning System (GPS) device 
performed repeated cross-shore transects at specified alongshore locations, which is 
described in more detail in Section 4.3. Owing to the fixed location of the stationary dye 
sensor array, the individual dye studies were performed at approximately the same tidal 
elevation, about 0.5 m above MSL, such that the dye sensors were located in the active 
surf zone. The dye was generally released between poles toward the northern end or in 
the middle of the array owing to the predominant southerly surfzone alongshore current. 

The location of the alongshore array in the active surf/swash zone resulted in the 
dye sensors being alternately submerged and not submerged within the sea-swell water 
level range, which varied due to tidal elevation and minimal infragravity motions (Figure 
32c-e), resulting in dropouts in measurements of dye concentration. Including these 
dropouts would bias the measurements toward low concentrations; therefore to account 
for this in the measured dye concentration time series, an envelope was fitted to the 
maxima of concentration peaks associated with the sea-swell waves (Figure 33) to 
describe the actual concentration of dye in the surf/swash zone as a function of time. No 
procedure allowed for 100% automation, so the envelope was manually selected after 
some initial automation, which also removed erroneous values. 

C. EULERIAN CURRENT MEASUREMENTS 

1. Steep Beach Surf/Swash Zone Characteristics 

The surf zone on the steep beach at CRSB has been described as being a large 
swash zone. When the EMCM array was inside of the surf/swash zone, the sensors 
became alternately submerged and not submerged, depending on their location in the 
water column, with the frequency of the incident waves. A one-minute time series of 
cross-shore velocity measured by each sensor in the EMCM array when it was inside of 
the surf zone is shown in Eigure 34 (top). At times just prior to and during wave breaking 


58 



(e.g., = 0 s), there were strong onshore (negative) veloeities measured by all sensors in 
the vertieal array. As the flow switehed direetion, the sensors nearest to the water surfaee 
(> 0.63 m above the sea-bed) immediately beeame dry and did not measure baekwash 
(positive) veloeities, while the sensor 0.43 m above the sea-bed (red line in Figure 34) 
measured part of the baekwash then beeame dry, resulting in a veloeity measurement of 
zero, and only the bottom-most sensors (0.03 and 0.23 m above the sea-bed, blue and 
green lines in Figure 34) measured offshore veloeities throughout the entire baekwash 
(Figure 34, top). Although the EMCM array was in a mean water depth of 0.7 m during 
this example time series, the measured veloeities still exhibit swash zone veloeity 
eharaeteristies. When the EMCM array was outside of the surf/swash zone, all of the 
sensors in the vertieal array were submerged and measured veloeities at all times (Eigure 
34, bottom). 

2. Undertow Inside and Outside Surf Zone 

The vertieal strueture of the eross-shore flow inside and outside of the surf zone 
was evaluated using half-hour averaged EMCM and ADCPl veloeity measurements 
(Eigure 35). When the EMCM array was inside the surf zone, an undertow profile was 
observed, with the eharaeteristie parabolie shape and maximum offshore veloeity at mid¬ 
depth, whieh is typieally observed within the surf zone on alongshore-homogeneous 
dissipative beaehes (Gareez Earia et al. 2000; Reniers et al. 2004a). The undertow 
veloeity ehanged from a small offshore veloeity throughout the water eolumn when the 
EMCM array was outside of the surf zone (Eigure 35b, solid eirele, square, and triangle 
symbols) to a strong offshore flow at mid-depth as the EMCM array transitioned to being 
inside of the surf zone (Eigure 35b, solid star and diamond symbols). This demonstrates 
that the undertow within the surf zone on this steep, refleetive beaeh was not destroyed 
by standing waves, whieh are eommon on steep, refleetive beaehes (Wright 1982). The 
mean measured alongshore veloeity was also greatest when the EMCM array was inside 
the surf zone (Eigure 35e, solid star and diamond symbols). 

Outside of the surf zone at ADCPl, loeated in 3 m water depth approximately 35 
m from the shoreline, an undertow profile was also observed, whieh had maximum 
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offshore flow near the surfaee and decreased with depth (Figure 35b, open symbols), 
consistent with undertow profiles measured on the inner shelf by Lentz et al. (2008). The 
magnitude of the undertow profile at ADCPl was similar to that measured with the 
EMCM array when it was outside the surf zone for corresponding times; the magnitude 
of the alongshore velocity at ADCPl was less than that measured by the EMCM array for 
corresponding times (Eigure 35c). These results indicate that outside of the surf zone 
there was a decrease in velocity with increasing water depth. 

D, DYE MEASUREMENTS 

1. Qualitative Dye Observations 

The dye was visually observed to mix rapidly throughout the water column and in 
the cross-shore immediately after the release, apparently saturating the surf zone within a 
few incident waves (Eigure 32b), and began seeping offshore of the surf zone shortly 
thereafter. The predominant pattern observed during the experiment consisted of the dye 
moving alongshore to the south, while also slowly being transported offshore. Eor times 
when the alongshore current was weak, the offshore movement of the dye was more 
visually noticeable because it stayed in approximately the same alongshore location, 
however it was found to move offshore at the same slow rate as when the alongshore 
current was strong. A summary of the visual dye patterns is given in Table 2. 

The movement of the dye was also examined with the dye concentration time 
series, C{t), measured at each alongshore sensor (Eigure 36). Unfortunately, there were 
times when the dye concentration was greater than the maximum limit of a dye sensor 
and the data were clipped, or some sensors were lost or removed when dye was still in the 
study area, so a complete dye mass balance could not be performed and observations 
from these dye releases are not discussed. Repeated cross-shore transects performed by 
swimmers during some of the dye releases also provide a measure of the offshore 
movement of the dye plume. The dye sensors in the alongshore array and/or the 
swimmers that were used to compute mass transport estimates during each dye release 
are given in Table 2. 
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The dye releases on yeardays 162, 169 and 172 oeeurred during relatively small 
offshore wave eonditions {Hrms < 0.5 m at ADCPl), whieh foreed relatively weak 
alongshore eurrents {Vavg < 0.1 m s'^ at ADCPl) so the dye appeared to move primarily 
offshore. During the dye releases on these yeardays, two distinet peaks in C{t) were 
measured by the dye sensor nearest to the dye release, indieating a separation in the dye. 
The C{t) measured by the alongshore array of dye sensors during dye release 162a is 
shown in Figure 36 (top) as an example of the weak alongshore eurrent ease. The dye 
initially reaehed the dye sensor nearest to the dye release loeation and moved offshore, 
seen by a narrow peak and sharp deerease in C{t). Then the remaining dye reaehed the 
sensor and gradually moved offshore, seen by a seeond, broader peak in C{t) measured by 
the same dye sensor. The C{t) measured by the subsequent dye sensors to the south in the 
alongshore array show broader profiles with smaller peak eoneentrations, indieating that 
the dye eontinued to spread alongshore and travel downstream, whieh ean also be seen in 
the time evolution of C(y) (Figure 36), whieh was eonsistent with visual observations. 
Due to the weak alongshore eurrent and redueed number of stationary dye sensors that 
measured eoncentration time series, quantitative estimates using the alongshore array of 
stationary sensors were limited. However, cross-shore profiles of dye concentration 
measured by the swimmers on yeardays 169 and 172 were used to quantify the cross¬ 
shore movement of the dye (Table 2). 

The dye releases on yeardays 165 and 166 occurred during relatively larger 
offshore wave conditions {Hms > 0.5 m at ADCPl), which forced relatively strong 
alongshore currents {vavg > 0.1 m s'^ at ADCPl) so the dye was released at the northern 
end of the array and appeared to move quickly alongshore to the south inside the surf 
zone with slight cross-shore movement. The C{t) measured by the alongshore array of 
dye sensors during dye release 166c is shown in Figure 36 (bottom) as an example of the 
strong alongshore current case. C{t) measured at the sensor nearest to the dye release 
showed a sharp rise followed by a rapid decay resulting in a relatively narrow 
concentration peak. As the dye reached the subsequent downstream alongshore sensors, 
the rise in concentration was less steep and it decayed more slowly, resulting in a broader 
concentration peak. The measured peak concentration moved quickly alongshore, and 
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also decreased with alongshore distance from the dye release, as seen in the time 
evolution of C(y) (Figure 36, bottom). There were no traces of dye measured by the 
sensors after 35 minutes during these releases. Swimmer data could not be used to 
quantify the cross-shore movement of dye during these releases because the dye moved 
too quickly past the locations of the swimmer transects and was not adequately captured. 


2. Alongshore Mixing and Transport Inside Surf Zone 


a. Surfzone A veraged Alongshore Current 


In-situ surfzone velocity measurements were not obtained during most of the dye 
releases. Therefore, a Eulerian measure of the alongshore current inside the surf zone was 
estimated as the alongshore velocity that the dye was advected past each dye sensor, 
using the dye concentration time series measured with the alongshore array of stationary 
dye sensors, which are assumed to represent the surfzone-averaged and depth-averaged 
concentration. The surfzone averaged alongshore current, V(y), at each alongshore 
instrument location, y, was determined by interpolating the alongshore velocity between 
sensors to the entire length of the alongshore array as a function of distance from the dye 
release location at y = 0. The alongshore velocity was computed by a forward- 
differencing scheme, defined as 


y= 

^Loi{y) 


(13) 


where Ay is the linear distance between sensors and Atmoiiy) is the difference in time 
associated with the first moment of the dye concentration time series (e.g. Figure 36), 
defined as 


^tC{y,t)dt 

l.m{y) = 7 - (14) 

jc{y,t)dt 

0 

where C(y,t) is the dye concentration as a function of alongshore distance, y, and time, t, 
from the dye release (y = 0, t = 0) and T is the duration of the dye measurements. Due to 
separations in the dye patch as it moved alongshore, it was found using tmoi, which 
represents the mean time associated with the dye concentration profile, provided a better 
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estimate of F(y) rather than using the time assoeiated with the leading edge of the dye 
plume, the time assoeiated with peak dye eoneentration, or the time assoeiated with the 
trailing 10% of the peak dye eoneentration. The effeet of alongshore diffusion on the 
estimates of tmoj are assumed negligible in the V(y) ealeulations. Estimates of V(y) tended 
to deerease slightly with distanee alongshore (Figure 37a), with estimates of the 
alongshore-averaged surfzone averaged alongshore eurrent, V , ranging from 0.1 to 0.3 
m/s for all dye releases (Table 2). The deerease in V(y) with distanee alongshore to the 
south eould be due to the dynamies assoeiated with the eurvature of CRSB and is 
investigated more in the Diseussion. 

The estimates of V from the dye inside the surf zone were eompared with the 
depth-averaged Eulerian veloeities measured by the EMCM array and ADCPl (Figure 
3Id, red). Estimates of V from the dye were eonsistent with the alongshore veloeities 
measured by the EMCM array when it was inside of the surf zone, (Figure 3 Id, green) 
whieh typieally oeeurred during low tides. The dye estimates of V inside the surf zone 
were faster than the depth-averaged alongshore veloeities measured outside the surf zone 
by ADCPl (Figure 3 Id, blue) and by the EMCM array when it was outside the surf zone 
(Figure 3Id, green), whieh typieally oeeurred during high tides. 

b. Alongshore Mass Transport 

Comparisons of the mass of dye measured at eaeh subsequent alongshore sensor 
during a dye release deseribe the alongshore adveetion and dispersion of the dye, and any 
loss of dye between alongshore sensors is attributed to eross-shore exehange. Assuming 
the dye is well-mixed (throughout the water eolumn and in the eross-shore) and is 
moving with V , the total mass of dye adveeted past eaeh alongshore dye sensor, 
Mmeasiy), measured in grams, where 1 ppb = 10'^ g m'^, was estimated by 

T 

M..M = ]KVC(y,l)dl ( 15 ) 

0 

where Ac is the eross-seetional area of the surf zone, whieh was eonsidered to be a wedge, 
and was eomputed as Ac = Q.5hhX„, where /?* = 1 m and is the water depth at the offshore 
edge of the surf zone, and Xw = 10 m and is the approximate surfzone width. The initial 
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amount of dye released in the surf zone, Mmeasiy = 0), was Mq = 46 g. The mass transport 
measured as a function of alongshore distance from the dye release was computed using 
estimates of V and is shown in Figure 37b. The Mmeasiy) decreased by a value of 1/e 
within 10 to 50 m from the dye release location, with a more gradual decrease in dye 
mass measured with increasing distance alongshore thereafter. 

The statistical significance of the alongshore decay of Mmeasiy) was examined by 
evaluating the mean dye concentration, C, measured at each dye sensor in the 

alongshore array for each dye release. The dependence of MmeosCy) on f C{y,t)dt , can 

J 0 


j N 

concentration measurements made in time, such that T = Ndt. Given J C{t)dt = '^C^dt 

n-l 

then [ C{t)dt = TC and M^^(y) = A^TC ■ 10^^, and the statistics associated with 

C can be applied to Mmeasiy)- Confidence intervals were estimated for C , and the decay 
with distance alongshore was found to be statistically significant at the 95% confidence 
level (not shown), which is used to demonstrate the statistical significance in the decay of 
dye mass Mmeasiy) with distance alongshore. 


^ C„ , and N is the number of discrete 


_ _ 1 

be written in terms of C , where C = — / C 

Nti 


c. Alongshore Diffusivity 


The rate of spreading, referred to as diffusivity (/c/), of the dye in the direction i 
X, y is estimated by 

1 dof 


K-,(t) = 


2 dt 


(16) 


where af is the second moment of the concentration profile. 




_ c{x,,t)d. 


(17) 


and fii is the dye center of mass given by the first moment of each dye concentration 
profile. 
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( 18 ) 


jX;C(Xp?)jX; 

jc[xj,t)dXj 

where C{xi,t) is the dye concentration profile in the z-direction (Clark et al. 2010). The 
alongshore spreading of the dye inside the surf zone from a narrow peak to a broad 
profile over a short alongshore distance (< 100 m) as a function of time (Figure 36) is 
described by the alongshore diffusivity, Ky. Equations (16), (17) and (18) are evaluated 
for i = y, using Ciy,t) measured by the alongshore array and integrating over the length of 
the alongshore array. The most dye releases occurred on yearday 166, with the fewest 
instances of data being clipped, so these results were ensemble averaged and used to 

2 

characterize the alongshore diffusivity at CRSB. Linear regression of versus tmoi 

(Figure 38) is used to fit a line to the data from all dye releases on yearday 166 (r^ = 0.44, 
significant in the 95% confidence interval), where half of the slope of the line represents 
Ky (Equation (18)). Values of Ky varied for the individual deployments. The mean Ky 
inside of the surf zone was found to be 0.5 m^ s’*. The alongshore diffusion inside the 
surf zone at CRSB is attributed to shear in the alongshore surfzone current. The 
distribution of the alongshore current from zero at the shoreline to a maximum velocity at 
a cross-shore location within the narrow surf zone (< 10 m) results in a considerable 
amount of shear in the alongshore current, which is a known mechanism for generating 
alongshore spreading and diffusion (Clarke et al. 2007). 

3. Cross-Shore Mixing and Transport across Surfzone Boundary 

The cross-shore movement of the dye was estimated using cross-shore dye 
concentration profiles measured by swimmers on yeardays 169b and 172a,b. People 
equipped with a dye sensor and GPS swam repeated, shore-normal, cross-shore transects 
from the shoreline to beyond the offshore edge of the dye plume and back in 
approximately the same alongshore location for the duration of the dye deployment 
(Figure 39, upper panels), resulting in a measure of the dye concentration as a function of 
cross-shore distance and time, C{x,t). Each swimmer performed 11 to 20 cross-shore 
transects (out and back) at a given alongshore location during a given dye release, 
resulting in 22 to 40 measurements of C(x,f) per swimmer per dye release. As the dye 
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spread offshore, the swimmers had to travel farther offshore, swimming between 5 to 45 
m offshore in a transect, taking on average 2 minutes to complete a transect. To 
summarize the offshore movement of dye, the large number of swimmer concentration 
profiles was averaged over 5-minute time windows to obtain mean concentration profiles 
as a function of cross-shore distance (Figure 39, lower panels). The cross-shore location 
of the dye plume, iJ.x{t), was found at each alongshore swimmer transect location for each 
5-minute mean C{x,t) measured by the swimmer using Equation (18) with i = x, where 
the limits of integration were from the shoreline to beyond the offshore extent of the dye 
plume (x = 80 m). 

The cross-shore dye concentration profiles at each alongshore location show the 
dye was initially shoreline-attached and saturated the surf zone in less than 5 minutes. 
The measured dye concentration decreased with time, and jix moved offshore as the dye 
concentration profile spread in the cross-shore. The Hx moved outside the surf zone after 
10 minutes, but remained within 20 m (2 surfzone widths) of the edge of the surf zone, 
indicating a decrease in the cross-shore dispersion of the dye outside the surf zone. Due 
to the presence of the shoreline, jix is expected to move offshore as the dye disperses in 
the cross-shore, however this could be due to seaward advection and/or plume widening 
by cross-shore diffusion (Clark et al. 2010; Spydell and Feddersen 2012). Unlike the 
shoreline-attached cross-shore concentration profiles and offshore jix movement with 
increasing alongshore distance from the dye source observed by Clark et al. (2010), 
which was attributed to cross-shore plume widening, the cross-shore concentration 
profdes observed at CRSB showed the peak concentration moving offshore as well as jix, 
indicating the dye was being advected offshore as well as being spread by cross-shore 
diffusion. 


a. Cross-Shore Diffusivity Inside Surf Zone 

Unfortunately, the swimmer transects were performed a considerable distance (> 
50 m) downstream of the dye release, and due to the steepness of the beach and the 
narrowness of the surf zone, the dye plume had already spread across the surf zone and 
moved offshore. Thus, the cross-shore spreading of the dye within the surf zone was not 
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accurately captured with the swimmer transeets. Due to the visual observations of rapid 
eross-shore mixing in the shallow surf/swash zone during the dye releases, the eross- 
shore diffusivity inside of the surf zone is attributed to breaking wave-indueed turbulenee 
generated within the surf zone by the intense breaking of plunging/surging waves near 
the shoreline on this steep beaeh (Feddersen 2007). The eross-shore turbulent eddy 
diffusion eoeffieient, Kx, inside of the surf zone is estimated by the depth-averaged, 
breaking wave-indueed turbulent eddy viseosity, Vt, using the Reynolds analogy whieh 
states that turbulent eddies transport momentum as well as mass, therefore momentum 
and tracer mass are assumed to mix identieally {kx ~ Vt) (Rutherford 1994; Brown et al. 
2009; Clark et al. 2011). Battjes (1975) found that v, assoeiated with eross-shore mixing 
is related to the surfzone-averaged wave energy dissipation, Sb, by 



where 77? is the signifieant wave height, whieh represents the turbulent mixing length, p is 
the density of seawater, and c is a ealibration factor (c = 0.7 (Brown et al. 2009)). The 
dissipation of ineident wave energy by breaking waves within the surf zone, Sb, is 
obtained using the wave energy balanee, assuming straight and parallel eontours, and 
integrating over the width of the surf zone. 



where E = pgH^ is the wave energy and Cgx is the wave group veloeity in the eross- 

shore direction determined using linear wave theory (Battjes 1975). Outside of the surf 
zone the ineident wave energy flux is eonstant, so measurements made by ADCPl in 3 m 
water depth are used to deseribe the wave energy flux at x = Xw, and the wave energy flux 
at the shoreline, x = 0, is assumed equal to zero. To aeeount for the refleetive nature of 
the steep beaeh, the ineident wave energy flux is multiplied by a faetor of (1 - 7?^), where 
7?^ is the refleetion eoeffieient. Henee, the surfzone-averaged wave energy dissipation is 
found by 
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( 21 ) 


(Ec\ (l-i?") 
p _ V g^/ADCPlV _ i 

where = 0.4 is the reflection coefficient for CRSB computed using the energy flux 
method by Sheremet et al. (2002). Estimates of Sb from Equation (21) were used to solve 
for Vt using Equation (19), where Kx ~ Vt. 

Recently, Eeddersen (2012) derived a surf zone turbulent dissipation rate scaling 
assuming weak mean currents and a vertically uniform length-scale, resulting in a non- 
dimensional surf zone eddy diffusivity coefficient, K. The scaling is also based on a 
balance between vertical turbulent diffusion and energy dissipation, analogous to Battjes 
(1975), however a vertically-uniform length-scale equal to a fraction, y, of the mean 
water depth, h, is used (/ = yh). The non-dimensional eddy diffusivity K as a function of 
height about the sea-bed, z, is 

^ (^) 4/3 /I 1/3 ( ^ 

= ( 22 ) 

'■UtJ 

where dFj dx = d ydx = e^ from Equation (21), A = e/ exp (a), where 

a = (3/2)*'^^C^/7 and = 0.57 (Eeddersen 2012). This parameterization allows the surf 

zone turbulence to be scaled to a fraction of the water depth that is representative of the 

small-scale turbulent eddies, whereas the calibration factor c used in Equation (19) was 

based on a best-fit to drifter diffusivity observations, which are due to horizontal motions 

at scales greater than the water depth. Due to the surf zone at CRSB effectively being a 

large swash zone, the mean water depth h used in Equation (22) was set equal to the wave 

height Hs measured at ADCPl, and the turbulent length-scale was assumed equal to the 

entire water column {y = 1). Values of the depth-averaged K{z) are assumed equal to the 

cross-shore turbulent eddy diffusivity Kx inside of the surf zone. Nearly identical results 

were obtained using Equation (19) and Equation (22) to estimate Kx, with Kx ranging from 

0.2 to 0.5 m^ s"' for the times of the dye releases. Estimates of Kx inside of the surf zone 

are the same order of magnitude as Ky found inside of the surf zone, and are similar to 

asymptotic surfzone Kx measured on dissipative beaches using drifters (Spydell et al. 

2007; Brown et al. 2009) and dye (Clark et al. 2010). 
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b. Surfzone Averaged Cross-Shore Current 

The cross-shore concentration profiles measured by the swimmers indicate that as 
the dye was spreading in the cross-shore due to turbulent diffusion, it was also being 
advected offshore with a mean cross-shore current. In order to compute the surfzone 
averaged cross-shore current at each alongshore instrument location U{y) for each dye 
release, an initial guess of f/ = 0.2 m s"' was used in Equation (12), using values 
representative of the cross-shore averaged surfzone (x = xJ2 = 5 m, h = 0.5 m, Ac = 
O.ShbXw = 5 m^), to estimate of C(y,t), which was used to compute the mass transport, 
Mcompiy), using Equation (15), and the error between Mmeasiy) and Mcompiy) was 
computed. Values of U were adjusted incrementally until the error between Mmeasiy) and 
Mcompiy) at the given alongshore instrument location was less than 0.25 g. This method 
was applied using the computed V from the given dye release, and estimated surfzone 
diffusivities of Kx = 0.5 m^ s"' and Ky = 0.5 m^ s'^ This resulted in a measure of the cross¬ 
shore advection as a function of alongshore distance y from the dye source, Uiy), which 
was largest closest to the dye release point and decreased by approximately half within 50 
m (Eigure 37c) due to the initial rapid cross-shore mixing, then effectively leveled off to 
an asymptotic value U = 0.05 to 0.09 m s"' (Table 2). The net values of U are similar to 
Uavg measured by the EMCM array when it was inside the surf zone, and are greater than 
the depth-averaged Uavg measured outside the surf zone by ADCPl (Eigure 31). The 
alongshore variation in Uiy) is examined in the Discussion. 

c. Cross-Shore Diffusivity Outside Surf Zone 

The spreading of the dye outside the surf zone is expected to be slower than the 
initial rapid cross-shore mixing inside the surf zone due to the lack of turbulence (Harris 
et al. 1963; Inman et al. 1971; Clark et al. 2010). The dye outside the surf zone is 
assumed to be well mixed over the vertical and behave independently from the processes 
inside the surf zone. In order to examine the spreading of the dye outside the surf zone 
using the swimmer transects during dye releases 169b, 172a, and 172b, the total area 
under C(x,0 was calculated iCtot), as well as the area under Cixoutd) corresponding to the 
portion of the dye concentration profde outside the surf zone, Xout > Xw (Cout). A ratio of 
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CoJCtot > 85% was used to conservatively determine when swimmer transects were 
predominantly outside the surf zone. The cross-shore diffusivity, k^, seaward of the surf 
zone was calculated using Equations (16), (17) and (18) for i = x, and integrating C{x,t) 
from the shoreline to the seaward limit of the transects for the seaward transects and vice 
versa for the shoreward transects. Estimates of cross-shore diffusivity outside the surf 
zone were found using linear regression of <7^ versus tmoi (Eigure 40), where half of the 
slope of the line represents Kx (Equation (18)). The correlation coefficients, r, were 
computed and those that were significant at the 95% confidence level were used. The 
mean Kx outside of the surf zone was found to be 0.01 m^ s"', indicating that seaward of 
the surf zone the dye is spreading slowly. Estimates of Kx outside of the surf zone are 20 
to 50 times smaller than that inside of the surf zone. 


The relative importance of the cross-shore diffusivity versus the cross-shore 
advection outside the surf zone is evaluated by examining the cross-shore advective and 
diffusive terms in Equation (11). The cross-shore advection at the edge of the surf zone is 
estimated as uhC = 0.7 ppb m^ s'^ using the net U = 0.07 m s’', ht = I m, and a 
maximum C = 10 ppb measured outside the surf zone. The cross-shore diffusivity outside 


• • 1 d{hC) dC dh , 2 i • i 

the surf zone is estimated as —-— = K^h^— + K^C— =0.02 ppb m s , using the 

ox ox ox 

same values for hb and maximum C outside the surf zone, and dCjdx ~ 1 ppb m'^ (where 
C decays from a maximum of 10 ppb to 0 ppb over a distance of 10 m), and an offshore 
slope 3/i/3x=l/10. This estimate indicates that outside the surf zone the cross-shore 
transport by advection is thirty-five times greater than the transport by diffusion. These 
results are consistent with the findings of Clark et al. (2010) for a dissipative beach, and 
suggest that the processes outside of the surf zone are disconnected from the surfzone 
dispersion processes, and that once the dye exits the surf zone it moves offshore as a 
constant patch. 


70 



E, 


DISCUSSION 


1. Stages of Mixing and Transport Inside Surf Zone on a Steep Beach 

Based on the qualitative and quantitative observations obtained with the 
alongshore array of dye sensors, two separate stages of mixing and transport are apparent 
inside the surf zone on the steep beaeh at CRSB (Figure 41). The results show differenees 
in the deeay of dye mass and cross-shore transport between the region within 
approximately 10 to 50 m of the dye release location and the region farther downstream 
alongshore (Figure 37). These differing stages of mixing and transport within the surf 
zone are analogous to the stages of mixing that occur in rivers, which occurs in three 
stages (Fischer et al. 1979). A tracer discharged into a river channel initially mixes in 
three dimensions until it becomes vertically well mixed. The tracer then continues to mix 
across and along the channel primarily by turbulence as it is advected downstream by the 
mean current. The region between the point of discharge and the downstream location 
where the tracer is fully mixed vertically and across the channel is considered the near¬ 
field. Downstream of this region is considered the far-field where longitudinal shear flow 
dispersion is the dominant mixing process. The stages of mixing in a surf zone can be 
compared to those in a river, where the surf zone is considered a channel with the 
boundaries being the shoreline and the offshore edge of the surf zone. 

The near-field stage of mixing and transport in the surf zone on the steep beach at 
CRSB is defined as the region between the point of the dye release and the location 
downstream alongshore where the dye was completely mixed vertically and in the cross¬ 
shore within the narrow surf zone. The extent of the near-field region is determined using 
results from dye release 172a when measurements of dye concentration are available both 
in the alongshore (by the alongshore array of stationary sensors) and in the cross-shore 
(by swimmers equipped with sensors). The dye appears completely mixed across the surf 
zone after 5 minutes after the release (Figure 39). At this time the dye is assumed to be 
completely mixed vertically as well because the surf zone is approximately 10 times 
wider than it is deep. The time of 5 minutes associated with the near-field stage of mixing 
corresponds to an alongshore distance of approximately 25 m from the dye release 
location, using the alongshore surfzone velocity of 0.08 m s‘^ measured during dye 
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release 172a. This alongshore loeation eorresponds with the loeation where Mmeasiy) 
deereased by a value of 1/e, whieh is used to define the near-field region on this steep 
beaeh for the different dye releases. Therefore, the near-field region is eonsidered to be 
eonservatively 50 m downstream of the point of the dye release (0 < y < 50 m), and the 
far-field region is eonsidered to be beyond 50 m downstream from the point of the dye 
release (y > 50 m). 

a. Near-Field Mixing and Transport 

In the near-field region, the dye was rapidly mixed throughout the water eolumn 
immediately after the dye release, and experieneed 2-D horizontal mixing and transport 
until the surf zone was saturated for y < 50 m (Figure 41). The dye was mixed in the 
eross-shore {kx) by turbulenee due to breaking waves, and spread in the alongshore {xy) 
due to shear in the alongshore current. The dye was advected alongshore with the 
surfzone averaged alongshore velocity V , while also being transported offshore by the 
surfzone averaged cross-shore velocity U(y), which was found to be relatively large 
(Figure 37c). Mmeasiy) was observed to decrease by a value of 1/e within y < 50 m. The 
strong decay in Mmeasiy) suggests that cross-shore diffusion, as well as cross-shore 
transport, is an important processes controlling the initial 2-D mixing and transport on 
this steep beach. Therefore, the surfzone mass transport in the near-field region is due to 
alongshore and cross-shore advection, as well as longitudinal dispersion and cross-shore 
diffusion, as summarized in Figure 41. 

b. Far-Field Mixing and Transport 

In the far-field region, the dye in the surf zone was completely mixed, and the 
transport was predominantly due to an advection balance (Figure 41). Mmeasiy) decayed in 
the alongshore more gradually, as the dye was observed to be advected alongshore with 
V , and the cross-shore transport leveled off to a relatively small constant value U 
(Figure 37). Additionally, the dye continued to be dispersed alongshore in the far-field 
(Figure 36). The decrease in cross-shore transport in the far-field region is due to the dye 
being completely well mixed, and therefore represents the steady-state cross-shore 
exchange that occurs between the surf zone and inner shelf on the steep beach. These 
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findings of mass transport being due to an advection balanee in the far-field on the steep 
beach at CRSB is consistent with the findings of Grant et al. (2005) for the advectively 
dominated transport conditions observed at the dissipative beaches in southern California. 

The exchange ratio VjU for CRSB was found to range from 1.14 to 3.50 (Table 
2), and is much smaller than that found by Grant et al. (2005) (y/C > 50 ). This indicates 
there was a greater amount of cross-shore exchange in the far-field on the steep beach at 
CRSB than measured at dissipative, alongshore current dominated beaches. The net 
offshore transport of dye at CRSB is consistent with the undertow velocity measurements 
obtained inside and outside of the surf zone during the field experiment. 

2, Offshore Pulses of Dye 

It was visually observed that the dye moved offshore in episodic pulses (not 
captured in quantitative dye measurements). To examine the possible transport associated 
with different frequency ranges, the cross-shore excursion length is estimated from the 
frequency-band-passed cross-shore velocities at ADCPl (Figure 42). Even though the 
velocities are largest for sea-swell (SS) waves (0.04 Hz < / < 0.2 Hz) and infragravity 
(IG) waves (0.004 Hz </< 0.04 Hz) (Figure 42, top), their duration is relatively short and 
their excursion is small (< 2 m) compared with very low frequency (VLF) motions 
(0.0005 Hz </< 0.004 Hz) (Figure 42, bottom). The excursions of the VFFs can be three 
times larger than SS and IG waves, and owing to their longer duration, 0(15 minutes), 
therefore the VFFs can be the mechanism for the episodic cross-shore exchange. This 
suggests that the dye pulses on a timescale of about 5 to 15 minutes due to VFF motions 
as it slowly moves offshore with the mean flow. Clark et al. (2010) similarly found that 
inside of the surf zone low-frequency horizontal rotational velocities with a length scale 
of the width of the surf zone (surfzone eddies) played an important role in cross-shore 
surfzone tracer mixing. 

It is hypothesized that the episodic offshore pulses of dye observed at CRSB are 
due to VFF surfzone eddies (SZEs). VFF SZEs are generated due to alongshore 
variations in wave breaking, which can occur as a result of several mechanisms: 1) 
alongshore variations in the underlying bathymetry (Bowen 1969; MacMahan et al. 
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2004), 2) short-crested wave breaking (Peregrine 1998; Clark et al. 2012), and 3) the 
interaetion of temporally and spatially varying wave groups (Johnson and Pattiaratehi 
2004b; Reniers et al. 2004b; Long and Ozkan-Haller 2009). 

Spatial variability in the underlying surfzone bathymetry ean result in alongshore 
variations in the breaking waves that foree offshore flows, typically in the form of rip 
eurrents (Bowen 1969; MaeMahan et al. 2004b). However, this is not believed to explain 
the offshore pulses observed at CRSB because the bathymetry eontours are straight and 
parallel and no alongshore variability was measured. The development of vortiees in the 
surf zone due to short-crested wave breaking and alongshore gradients in wave height has 
been investigated theoretically (Peregrine 1998) and observed in the field (Clark et al. 
2012). Adjacent regions of breaking and non-breaking along individual wave erests 
results in non-uniform forcing that generates rotational eddies at the end of the breaking 
wave erest, whieh ean have positive or negative vortieity, and is a meehanism for 
dispersion near the shoreline (Clark et al. 2012). Alongshore gradients in the breaking 
short-erested waves was observed at CRSB, however Clark et al. (2012) found that the 
vortieity generated by short-crested wave breaking decayed over 20 to 60 s, whieh is a 
shorter time-scale than the offshore pulses of dye that was observed. 

VLF SZEs ean also develop as a funetion of wave-group foreing assoeiated with 
direetionally-spread short waves, whieh foree intermittent eirculation eells that persist on 
0(15 minutes) and have the ability to transport material outside the surf zone (MaeMahan 
et al. 2010b). On an alongshore uniform beaeh, numerical model results show that 
alongshore variations in wave breaking can result from an ineident wave speetrum that 
ineludes a range of wave frequeneies and direetions, where the interaetion of the 
individual wave components ereates temporally and spatially varying wave heights (wave 
groups) (Reniers et al. 2004b; Long and Ozkan-Haller 2009). Long and Ozkan-Haller 
(2009) found that vortiees were generated for oblique wave groups on an alongshore 
uniform beach, and the resulting vortiees were advected by the mean surfzone eurrent. 
Field observations under direetionally-spread waves found vortex pairs manifest as 
transient rip eurrents (Johnson and Pattiaratehi 2004b). These model and field results 
show vortieity motions at time seales 0(5 min), similar to the time seale of the pulses 
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observed at CRSB. The formation of these VLF SZEs is not steady state, and it is 
hypothesized that the eonstruetive and deeonstructive interaetion of these VLF SZEs are 
the mechani s m for the episodic cross-shore exchange observed at CRSB. 

3. Effect of Local Bathymetry on Alongshore Surfzone Transport and 
Cross-shore Exchange 

The observed decrease in V with distance alongshore to the south (Figure 37a) 
could be explained by the decrease in wave angle relative to the shoreline with distance 
alongshore to the south due to the curvature of the bay. Waves coming from the north and 
driving the alongshore current to the south actually tend to be more shore-normal at the 
southern end of the experiment area and therefore force weaker alongshore currents. In 
numerical model simulations investigating biological transport during the CRSB 
experiment under the same conditions discussed here, Fujimura et al. (2013) found that 
the alongshore current is to the south and follows the contours of the embayed beach until 
approximately y = 100 m, where there is a decrease in the alongshore velocity and the 
current moves predominantly offshore. The southern end of the alongshore array of dye 
sensors is located where this change in current from alongshore to offshore occurs, which 
explains the decrease in V(y) with distance alongshore. 

4, Effect of Local Wave Conditions on Cross-shore Exchange 

It has been noted that mixing and mass transport on natural beaches is likely 
sensitive to local wave conditions (Boehm et al. 2003; Grant et al. 2005). On dissipative 
beaches, waves approaching the shoreline at an angle drive alongshore currents, and as 
the incident wave angle approaches shore-normal rip currents are formed. During 
obliquely incident, alongshore-current wave conditions, mass transport is advectively 
dominated (Grant et al. 2005). Additionally, cross-shore wave-driven mixing across a 
surf zone with alongshore current shear enhances the alongshore spreading of dye (Clarke 
et al. 2007). During shore-normal wave conditions, rip currents form resulting in 
increased cross-shore advection, and turbulent diffusion due to breaking waves becomes 
more important. Therefore, on dissipative beaches, the exchange ratio VjU is largest for 
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obliquely incident waves during alongshore-current conditions, and decreases as the 
incident wave angle decreases and a greater amount of cross-shore exchange occurs. 

Similarly, on steep beaches, obliquely incident waves drive alongshore currents; 
however, in the near-field, there is a greater amount of cross-shore exchange, likely due 
to wave breaking, and therefore, U is relatively large. In the far-field, after the dye is 
mixed across the surf zone, U decreases and levels off to a constant, and the mass 
transport is advectively dominated and is similar to that observed by Grant et al. (2005) 
on a dissipative beach under alongshore current conditions. In the far-field, the exchange 
ratio VjU is generally less than that on a dissipative beach for oblique wave angles, 
indicating more cross-shore transport occurs on steep beaches under alongshore current 
conditions. 

F. CONCLUSIONS 

Nearshore processes of waves and currents and their impact on mass transport and 
cross-shore exchange were observed for a steep, reflective beach at Carmel River State 
Beach, California. Owing to the steep slope (1/7) of the beach, plunging/surging waves 
broke close to the shoreline resulting in an intense shore break and a large swash zone 
rather than a typical surf zone generally associated with dissipative beaches. Fluorescent 
Rhodamine dye was released in the surf zone and an alongshore array of stationary dye 
sensors and repeated cross-shore transects performed by swimmers equipped with dye 
sensors at multiple alongshore locations were used to measure the temporal and spatial 
dye concentration evolution and to examine the alongshore and cross-shore mixing and 
transport associated with the hydrodynamics on the steep, reflective beach. 

Visual dye observations indicate that the dye quickly mixed vertically and 
dispersed in the cross-shore, completely saturating the narrow surf zone after a few 
incident waves. Dye concentration time series measured by an alongshore array of 
stationary dye sensors inside the surf zone revealed sharp-rising, narrow spikes of dye 
concentration that decayed rapidly near the dye release location, which broadened and 
decayed more slowly and decreased in peak concentration with distance alongshore, 
downstream from the dye release location. Mixing and transport processes occur in two 
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states on the steep beaeh, whieh are summarized in a coneeptual diagram in Figure 41. In 
the near-field, within approximately 50 m alongshore from the dye release location, dye 
was mixed in the cross-shore by turbulent diffusion due to breaking waves {kx = 0.2 to 
0.5 m^ s"') that resulted in a relatively large amount of cross-shore exchange {U{y < 50 m) 
= 0.1 to 0.4 m s"'), and was also advected downstream with the mean alongshore surfzone 
current {V = 0.1 to 0.3 m s"') and spread in the alongshore due to shear in the alongshore 
current {xy = 0.5 m^ s'^). In the far-field, once the dye was completely mixed vertically 
and in the cross-shore, the cross-shore transport leveled off to a constant {U = 0.05 to 
0.09 m s'^). Dye continued to be transported alongshore with the mean alongshore 
surfzone current and spread in the alongshore, while also being slowly transported 
offshore. Even though the offshore velocity is small, it represents a continuous loss and a 
substantial net offshore transport. Outside the surf zone, the cross-shore diffusion 
coefficient Xx was 0.01 m^ s"' indicating the dye was spreading more slowly than it was 
being transported offshore, and that it was essentially moving offshore as a constant 
patch. The offshore movement of the dye plume also pulsed on the order of 5 to 15 
minutes, attributed to VLF motions. 

The mixing and transport processes observed on this steep beach resemble those 
observed on dissipative beaches; however, some differences do exist. The mass transport 
measured on several dissipative beaches in Southern California was found to be 
dominated by advection, and the dye mass was observed to decay exponentially with 
distance alongshore (Inman et al. 1971; Boehm 2003; Grant et al. 2005). On the steep 
beach, in the near-field there was a large amount of cross-shore exchange initially and 
advection and diffusion processes were important. In the far-field, the mass transport was 
dominated by advection and was similar to that on dissipative beaches under alongshore 
current conditions; however, the far-field solution on the steep beach was reached in a 
much shorter alongshore distance than on dissipative beaches. The cross-shore exchange 
ratio V jU in the far-field for CRSB ranged from 1 to 4, and was an order of magnitude 
less than that found by Grant et al. (2005) for a dissipative beach with similar wave and 
alongshore current conditions. Therefore, a greater amount of cross-shore exchange 
across the surfzone boundary was observed on the steep, reflective beach. 
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V. SUMMARY 


Cross-shore exchange on two contrasting natural beaches has been described: 1) a 
dissipative, alongshore-inhomogeneous, rip-channeled beach, and 2) a reflective, 
alongshore-homogeneous, steep beach. The morphology, beach slope, wave conditions, 
and resulting current patterns differ significantly on these beaches, and therefore so do 
the mixing and transport processes and the subsequent material exchange and retention. 

Material transport was investigated on a rip-channeled beach at Monterey, 
California. The surfzone morphology consisted of transverse bars, which were relatively 
dissipative with a slope of about 1/100, and quasi-periodic incised rip channels, with 
straight and parallel contours outside the surf zone and an offshore slope of about 1/20. 
Rip currents were observed to occur during low tides and normally incident wave 
conditions, with increased wave energy resulting in stronger rip current flows. During 
these conditions, the surf zone width was approximately 100 m, on average. A 
combination of Eulerian and Lagrangian instruments were used to measure the temporal 
and spatial evolution of the rip current flow inside and outside the surf zone. These 
observations, in combination with results from previous studies, give a comprehensive 
description of the 3-D transport and exchange between the surf zone and the inner shelf 
due to rip currents, which is summarized here. 

The rip current circulations inside the surf zone form when incident waves 
approach the beach shore-normally, and the alongshore perturbations in the bottom 
topography induce alongshore variations in the depth-limited wave breaking, resulting in 
alongshore gradients in pressure and momentum driving water from the shoals to the rip 
channels (Bowen 1969, Haller et al. 2002). The directionally spread incident waves 
deposit momentum into the surf zone at the wave group time scale (0(10 minutes)), and 
the spatial variations in the breaking wave groups force VLF motions inside the surf zone 
(Renters et al. 2007, Long and Ozkan-Haller 2009). These energetic motions for 
frequencies less than 0.004 Hz are outside the gravity restoring region in frequency- 
wavenumber space, indicating they are horizontal vortices. These large-scale vortices, or 
surf zone eddies (SZEs) are coupled to the underlying bathymetry and are predominantly 
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contained within the surf zone (MaeMahan et al. 2004b). However, as eounter-rotating 
SZEs oscillate in the eross-shore direction inside the surf zone, they oeeasionally 
eonverge at the loeations of rip ehannels, and the presenee of the VLF rip eurrent pulses 
eause the rip eurrent to extend outside the surf zone (Reniers et al. 2009, 2010). 

The VLF eomponent of the rip eurrent flow aecounted for a signifieant portion of 
the total energy inside and outside the surf zone (37% and 50%, respectively), with a 
reduetion in the VLF rip eurrent intensity as it moved offshore, eonsistent with previous 
observations (MaeMahan et al. 2004b). The VLF pulses in the near-surfaee rip eurrent 
flow responsible for the rip eurrent exits were found to be infrequent and aperiodie, both 
inside and outside the surf zone. Evaluating 3-hour time windows eentered around low 
tide, when rip eurrents are known to be strongest, signifieant VLF pulses (greater than 
11 +2a) exeeeded 0.7 and 0.2 m s’' and oecurred one to twenty times inside and outside 
the surfzone, with the number and amplitude of significant VLF pulses inereasing with 
inereasing wave energy. On average, the significant VLF pulses oecurred for 
approximately 10% of the time. This is eonsistent with observations of relatively high 
amounts of surfaee surfzone retention (Castelle et al. 2010, MaeMahan et al. 2010b). The 
strength of the surfaee VLF rip eurrent pulses was found to eontrol the magnitude and 
vertieal strueture of the mean rip eurrent flow. Inside the surf zone the VLF pulses added 
to the strength of the baekground rip eurrent flow, and outside the surf zone, where there 
was no baekground flow, the stronger VLF pulses resulted in stronger rip eurrent flows 
and more shear in the vertieal. 

As the rip eurrent moved seaward, the VLF pulses inside and outside the surf 
zone did not oeeur simultaneously and were not eorrelated. The mean VLF eomponent of 
the veloeity measured inside the surf zone was significantly correlated with the offshore 
wave height, indieating the VLF pulses inside the surf zone were dependent upon the 
ineoming wave energy, due to the foreing of the VLF motions by the spatially varying 
breaking wave groups. However, outside the surf zone the mean VLF eomponent of the 
veloeity was not eorrelated with the offshore wave eonditions, supporting the notion that 
the VLF pulses originate inside the surf zone but travel offshore randomly. Additionally, 
the mean VLF eomponent of the veloeity was significantly correlated with the loeal wave 
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height to water depth ratio both inside and outside the surf zone. These results suggest 
that the VLF pulses outside the surf zone are independent of those inside the surf zone. 
Outside the surf zone, the VLF rip eurrent pulses moved offshore like a surfaee jet, either 
direetly offshore or at an angle. However, at times there was evidenee that the rip eurrent 
VLF pulse was moving offshore as a detaehed eddy, eonsistent with numerieal model 
results by Renters et al. (2010). 

The rip current flow was depth-uniform inside the surf zone, with minimal 
rotation throughout the water column. As the rip current exited the surf zone it moved 
over the surfzone break as a surface jet with its vertical extent remaining the same as at 
the edge of the surf zone. As the water depth increased outside the surf zone, the rip 
current transport remained near the surface, resulting in a surface dominated rip current 
flow offshore. There was rotation over depth in the rip current flow outside the surf zone, 
with an -'50° difference between surface and bottom, on average. The material transport 
offshore at the surface by the rip current flow was then mixed vertically downward to a 
depth corresponding to the thermocline, often times encompassing greater than half of the 
water column. 

Material that exited the surf zone at the surface due to VLF rip current pulses was 
observed to move seaward and ultimately return shoreward in an arcing pattern, re¬ 
entering the surf zone over shoals, with ultimately no material being completely removed 
from the nearshore system. The surface material either moved seaward and returned 
sharply back shoreward a short alongshore distance from where it exited, resulting in 
locally-contained cross-shore exchange, or moved seaward and travelled farther in the 
alongshore direction as it gradually moved shoreward, resulting cross-shore and 
alongshore exchange. The cross-shore extent was typically one to two surfzone-widths Lx 
beyond the surfzone boundary. The rate of spreading and mixing {k^xx = 0.27 to 1.74 m^ 
s'^ K^yy = 0.36 to 5.01 m^ s'') outside the surf zone was slower than that measured inside 
the surf zone during previous studies (Brown et al. 2009), suggesting there is less shear in 
the flow field outside the surf zone. 

The probability distributions of particle movements at the surface outside the surf 

zone were basically Gaussian, indicating there was a relatively equal chance that a 
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particle anywhere on the inner shelf would move seaward or shoreward, consistent with 
the qualitative observations that no drifters left the nearshore system. The cross-shore 
velocities of seaward and shoreward moving particles at the surfaee had a maximum at 
the edge of the surf zone, and were approximately equal in magnitude as a function of 
cross-shore location, resulting in a total Lagrangian velocity of approximately zero. 
Surface material moving seaward decelerated after exiting the surf zone through rip 
ehannels, while surface material moving shoreward accelerated approaching and re¬ 
entering the surf zone over shoals. The magnitude of the shoreward moving material 
outside the surf zone was greater than the near-surface theoretieal Stokes drift; therefore, 
Stokes drift was not solely responsible for moving material shoreward and baek into the 
surf zone. 

These results support a new theory of how material is exchanged and conserved 
within the nearshore region on a rip-ehanneled beaeh. 

1) Initially, short waves deposit momentum into the surf zone at the wave group 
timescale. 

2) Then spatial variations in the wave groups force VLF motions and rip current 
circulations inside the surf zone, where the rip current cells are topographically fixed in 
the alongshore direction, but oscillate in the cross-shore direction within the surf zone. 

3) The episodie convergence of these SZEs eauses the rip eurrent to pulse and 
expel material offshore through the rip channels. 

4) Material is transported seaward like a jet at the surfaee, with the depth-uniform 
rip current flow inside the surf zone becoming surface dominated as it moves over the 
surfzone break, remaining between the surface and the depth at the edge of the surf zone. 

5) Outside the surf zone, material at the surface spreads laterally slowly, and is 
also mixed vertically downward to a depth eontrolled by the thermoeline, often 
encompassing greater than half the water column. 
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6) Ultimately, the surfzone rip eurrent eireulations act to pull the material back 
onshore into the surf zone over shoals to satisfy a mass balance, with little or no material 
being completely removed from the nearshore system on the rip-channeled beach. 

Material transport was also investigated on a steep beach at Carmel River State 
Beach, California. The bathymetry was alongshore uniform, and had a cross-shore profile 
consisting of a linear beach face with a slope of about 1/7, connected to a steep 1/3 
subaqueous beach step, and continuing offshore with a slope of about 1/19. Due to the 
steepness of the beach, plunging/surging waves broke very close to the shoreline, 
resulting in a narrow surf zone that was approximately 10 m wide. The intense shore 
break was effectively the run-up that alternately wets and dries the beach; therefore, the 
surf zone on this steep beach is considered to be more of a large swash zone than a 
typical surf zone generally associated with dissipative beaches. Lagrangian measurements 
were obtained to quantify the movement and exchange of material. 

Mixing and transport processes occur in two states on the steep beach. 1) In the 
near-field, within approximately 50 m alongshore from the source, material quickly 
mixed vertically and dispersed in the cross-shore, completely saturating the narrow surf 
zone after a few incident waves. The material was mixed in the cross-shore by turbulent 
diffusion due to breaking waves {kx = 0.2 to 0.5 m^ s'*) that resulted in a relatively large 
amount of cross-shore exchange {U{y < 50 m) = 0.1 to 0.4 m s'*), and was also advected 
downstream with the mean alongshore surfzone current {V = 0.1 to 0.3m s'*) and spread 
in the alongshore due to shear in the alongshore current (% = 0.5 m^ s'*). 2) In the far- 
field the material was completely mixed vertically and in the cross-shore, and advection 
processes were dominant. The material continued to be transported alongshore with the 
mean alongshore surfzone current and spread in the alongshore, while also being slowly 
transported offshore, with the cross-shore transport leveling off to a constant {U = 0.05 
to 0.09 m s'*). This resulted in a continuous loss and a substantial net offshore transport. 
Outside the surf zone, the cross-shore diffusion coefficient Kx was 0.01 m^ s'* indicating 
material was spreading more slowly than it was being transported offshore, and that it 
was essentially moving offshore as a constant patch. 
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The natural beaches studied in this work represent two vastly different beach 
states, which commonly occur in nature. In reality, a given beach can change beach states 
over time based on the incident waves and resulting currents. Additionally, changes in the 
incident wave characteristics on a given beach can result in greatly different current 
patterns and exchange processes. Therefore, measuring mixing and transport for all types 
of beach states will provide information on how material might be transported as the state 
of a beach or the current patterns on a beach progresses. For example, rip currents are 
known to form when waves approach shore-normally on a beach with transverse bars and 
incised rip channels, however when the waves approach the same beach obliquely, an 
alongshore sinuous current pattern develops where the flow follows the underlying 
bathymetry contours (MacMahan et al. 2010b). The mixing and transport associated with 
these two current patterns are likely to be different. Forthcoming analysis examining field 
measurements of exchange processes during varying wave and current conditions on the 
rip-channeled beach will provide further insight into the potentially variable mixing and 
transport processes. 

Ultimately, future work investigating cross-shore exchange on these beaches will 
require the use of a numerical model, supplemented by the field observations, to quantify 
temporal and volumetric exchange estimates. While field measurements are crucial for 
understanding current dynamics, the number of observations points is ultimately limited 
spatially, therefore a numerical model allows for the transport and mixing processes to be 
evaluated everywhere in space (3-D), as well as in time. However, the model forcing 
parameters and subsequent hydrodynamic output must be validated by the field 
measurements to ensure the transport and mixing processes are accurately implemented 
(i.e., accurate wave-group forcing on the rip-channeled beach is critical for producing 
VLFs in the surfzone, which were shown to control the magnitude and amount of rip 
current transport). Ultimately, numerical models can be used to estimate surfzone 
flushing times and volumetric transport rates, based on the quantified mixing and 
transport processes, on all types of natural beaches. 
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TABLES 


Table 1. Summary of conditions observed during the drifter 
deployments at Monterey, California, grouped by observed drifter 
pattern. Tidal elevation was obtained from NOS tidal station 9413450; 

wave parameters and surfzone characteristics were obtained from 
offshore ADCP in 13 m water depth when data was available, otherwise 
were obtained from NDBC buoy 46240; wind data was recorded atop a 
10 m tower onshore of the field site; depth-averaged cross-shore velocity 
Uavg was measured by ADCP 3, and alongshore velocity Vavg was 
computed from the drifters; cross-shore location of the surfzone 
boundary Xsz and surfzone width were determined by shoaling waves 
from the ADCP in 13m water depth to the shoreline using linear wave 

theory. 


*sensor was buried 


yearday 

drifter 

pattern/ 

exchange 

type 

tidal 

elevation 

(m) 

Hs 

(m) 

Tp 

(s) 

\^s\ 

(N/m^) 

Uavg 

(m/s) 

Vavg 

(m/s) 

Az 

(m) 

A 

(m) 

maximum 

drifter 

cross¬ 

shore 

extent 

(t/x) 

122 

cross-shore 

-0.87 

0.43 

10.6 

0.00 

0.01 

0.05 

147 

78 

3.6 

129 

cross-shore 

-0.30 

1.31 

9.5 

0.02 

0.70 

0.00 

180 

123 

3.8 

132 

cross-shore 

-0.51 

0.77 

7.6 

0.04 

0.35 

-0.04 

157 

97 

3.7 

125 

cross-shore 
& alongshore 

-0.14 

0.97 

11.5 

0.02 

0.03 

-0.14 

176 

122 

1.9 

126 

cross-shore 
& alongshore 

-0.05 

1.22 

12.6 

0.00 

0.20 

-0.15 

174 

122 

2.8 

130 

cross-shore 
& alongshore 

-0.02 

0.92 

9.2 

0.04 

* 

0.13 

162 

111 

3.3 

137 

cross-shore 
& alongshore 

-0.54 

0.55 

8.2 

0.01 

0.08 

0.13 

142 

81 

2.5 

135 

none 

-0.63 

0.63 

6.7 

0.04 

0.04 

0.02 

149 

86 

3.7 

136 

none 

-0.70 

0.86 

8.3 

0.01 

0.07 

0.02 

169 

105 

6.2 
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Table 2. Description of dye releases and observations at Carmel River 
State Beach, California, and sensors used to estimate mixing and 
transport. Surfzone averaged alongshore current V is equal to the 
alongshore-mean V(y) computed using Equation (13). Surfzone averaged 
cross-shore current U is for the far-field and is computed iteratively 

using Equation (12). 


Release 

Dye Release 
Location 
(between 
poles) 

Visual Dye 
Observations 

Sensors Used 
for Flux 
Calculations 

E (m s-') 

U (m s"') 

Exchange 

Ratio, 

y/c 

162a 

3 &4 

predominantly 

offshore 

1-8 

0.11 

0.06 

1.83 

165a 

3 &4 

predominantly 

alongshore 

1-6, 8 

0.12 

0.07 

1.71 

165b 

1 &2 

predominantly 

alongshore 

1, 3-6, 8 

0.14 

0.05 

2.80 

165c 

1 &2 

predominantly 

alongshore 

1, 3-6, 8 

0.10 

0.05 

2.00 

165d 

1 &2 

predominantly 

alongshore 

1,3-6 

0.13 

0.06 

2.17 

166a 

1 &2 

predominantly 

alongshore 

1, 3-6, 8 

0.18 

0.08 

2.25 

166b 

1 &2 

predominantly 

alongshore 

1, 3-6, 8 

0.22 

0.07 

3.14 

166c 

1 &2 

predominantly 

alongshore 

1, 3-6, 8 

0.28 

0.08 

3.50 

166d 

2&3 

predominantly 

alongshore 

1,4-6, 8 

0.28 

0.09 

3.11 

166e 

2&3 

predominantly 

alongshore 

1, 3-6, 8 

0.27 

0.08 

3.38 

166f 

2&3 

predominantly 

alongshore 

1, 3-6, 8 

0.18 

0.07 

2.57 

169b 

5&6 

predominantly 

offshore 

3 swimmers 




169c 

4&5 

offshore and 
alongshore 

4-8 

0.09 

0.05 

1.80 

172a 

4.5 & 5 

predominantly 

offshore 

4.5-8, 

3 swimmers 

0.08 

0.07 

1.14 

172b 

4.5 & 5 

predominantly 

offshore 

2 swimmers 
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FIGURES 



alongshore (m) 

Figure 1. Measured bathymetry in the local coordinate system, where the 
contours are labeled with the elevation in meters relative to MSL. 
ADCP locations are shown by the black circles. Dashed black line 
represents the approximate edge of the surf zone. 
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time (yeardays) 


Figure 2. Wave eonditions measured by the ADCP in 13 m water depth of (a) 
root-mean-square wave height (b) mean wave period Tmoi, (c) 
and mean wave direction , (d) cross-shore wind stress computed 
for winds measured onshore of the field site, and (e) tidal elevation 
measured at NOAA/NOS tidal station 9413450 relative to MSL. 
Gray shaded regions represent 3-hour records, centered around low 
tide, used in the analysis of the VLF rip current pulses. 
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Figure 3. Conditions of (a) offshore root-mean-square wave height Hms (left 
axis) and tidal elevation (right axis), shown for reference, and the 
hourly mean cross-shore velocities measured by (b) ADCPout and (c) 
ADCPin, as a function of relative-depth, dih. Colorbar in (b) 
represents the cross-shore velocity magnitude for (b) and (c). Red 
dots at the bottom in (a) and (c) represent times when ADCPin was 
located inside the surf zone {HmJh > 0.4). 
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Figure 4. Conditions of (a) offshore root-mean-square wave height Hrms (left 
axis) and tidal elevation (right axis), shown for reference, and the 
hourly mean alongshore velocities measured by (b) ADCPout and (c) 
ADCPin, as a function of relative-depth, dih. Colorbar in (b) 
represents the alongshore velocity magnitude for (b) and (c). Red 
dots at the bottom in (a) and (c) represent times when ADCPin was 
located inside the surf zone {HrmJh > 0.4). 
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Figure 5. Conditions of (a) offshore root-mean-square wave height Hms (left 
axis) and tidal elevation (right axis), shown for reference, and the 3- 
hour averaged root-mean-square VLF total horizontal velocity, 
Urms.vLF, mcasurcd by (a) ADCPout and (b) ADCPin, as a function of 
relative-depth, d/h. Colorbar represents the Urms.viF magnitude, note 
the different scales. Red dots at the bottom in (a) represent times 
when ADCPin was located inside the surf zone {HrmJh > 0.4). 
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Figure 6. Examples of low-pass filtered eross-shore veloeity, uvlf, as a 
funetion of relative-depth, d/h, measured at (left column) ADCPin 
and (right column) ADCPout during 3HRLTs with (row 1 and 2) 
high, (row 3 and 4) intermediate, and (row 5 and 6) low VLF energy 
conditions. Colorbars represent UfTF magnitude. 
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Figure 7. Frequency distribution of uvlf measured during all 3HRLTs at (top) 
ADCPout and (bottom) ADCPin. Vertical gray dashed lines represent 
the threshold used to determine significant VLF pulses, which was 
two times the standard deviation above the mean (^ + 2 (t) . 
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Figure 8. Examples of low-pass filtered eross-shore veloeity, uvlf, measured 
in the reference level at ADCPm (blue) and ADCPout (red) during 
3HRLTs with (row 1) high, (row 2) intermediate, and (row 3) low 
VLF energy conditions. Colored horizontal dashed lines represent 
the threshold used to determine significant VLF pulses, and colored 
arrows indicate time and amplitude of the significant pulses. 
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Figure 9. Conditions of (a) offshore root-mean-square wave height Hrms, (b) 
number of significant VLF pulses, Npuise, (c) mean amplitude of 
significant VLF pulses, Amean,pulse, and (d) the fraction of the time 
that the VLF pulses were occurring, defined as the sum of the 
individual pulse durations divided by the 3-hour record length, 

’ Pleasured at ADCPi„ (blue) and ADCPout (red) 

during the 3HRLTs (indicated by gray shaded regions in (a)) 
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Figure 10. Vertical profiles of the low-passed filtered cross-shore velocity, uvlf, 
as a function of relative-depth, dih, measured at (left) ADCPm and 
(right) ADCPout , bin-averaged as a function of the mean amplitude 
of significant VLF pulses, Amean.puhe, measured in the reference level 
(colored lines, defined in legend). 
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Figure 11. Vertical profiles of the mean cross-shore velocity, Umean, as a 
function of relative-depth, dih, measured at (top row) ADCPom 
and(bottom row) ADCPm, bin-averaged as a function of (column 1) 
the maximum amplitude of significant VLF pulses, Amax.puke, 
(column 2) the mean amplitude of significant VLF pulses, Amean.puke, 
(column 3) the number of significant VLF pulses, Npuise, and 
(column 4) the duration of significant VLF pulses, > 

measured in the reference level (colored lines). 
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Figure 12. The cross-correlation R{t = 0) as a function of relative-depth, d/h, 

of the (a, c) cross-shore and (b, d) alongshore velocity measured at 
each depth with the relative velocity measured in the surface 
reference level, measured at (a, b) ADCPout and (c, d) ADCPm. 
Colorbars show correlation values. 
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Figure 13. Examples of mean horizontal VLF velocity vectors measured at 
ADCPout at various relative-depths, dih (colored lines, defined in 
legend) during 3HRLTs, showing differences in angle between the 
surface and sub-surface VLF velocities. 
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Figure 14. The eross-correlation of the cross-shore (blue) and 

alongshore (red) velocity measured at ADCPout with ADCPin. Gray 
shaded regions represent SFIRLTs and horizontal black lines 
represent significant correlation values for each 3HRLT. 
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Figure 15. Linear regression of the 3-hour averaged root-mean-square VLF 
total horizontal velocity, Umis.vLF, measured in the reference level 
and the 3-hour averaged (a) offshore root-mean-square wave height 
Hrms, (b) offshore mean wave period Tmoi, (c) offshore mean wave 
direction 0^^^^, (d) mean total horizontal velocity, Umean, measured in 
the reference level, (e) local water depth, h, and (f) local root-mean- 
square wave height normalized by the local water depth, HrmJh, 
measured by ADCPin (blue) and ADCPout (red) during the 3HRLTs. 
Colored lines represent bet-fit lines, and correlation coefficients, , 
given in legends. Significant correlation coefficient is = 0.35 at the 
95% significance level for all cases. 
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Figure 16. Example of low-pass filtered (top) eross-shore veloeity, uvlf, and 
(bottom) alongshore veloeity, vvlf, as a function of relative-depth, 
dih, measured at ADCPout during a 3HRLT on yearday 137. 
Colorbars represent uvlf and vvlf magnitude. 
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Figure 17. (a) 3-hour mean eross-shore veloeity profiles eorresponding to times 

of strong VLF pulses measured at ADCPm and ADCPout showing the 
rip eurrent moving offshore as a jet. (b-d) Vertieal measurements of 
dye eoneentration (solid) and de-meaned temperature (dashed) at 
varying eross-shore loeations for time intervals after the dye release 
of (b) 0 to 30 minutes, (e) 30 to 60 minutes, and (d) 60 to 90 
minutes. Mean eross-shore bathymetry profile in a rip ehannel (solid 
blaek line) and on a shoal (dashed blaek line). 


103 



































a) Stokes transport b) inner shelf c) surf zone 

z z z 



d) conventional rip current transport e) observed rip current transport 

®E,rip 



Figure 18. (a) Two-dimensional (x, z) onshore volume Stokes transport Qst in a 

Lagrangian (particle following, red line) and a Eulerian (fixed in 
space, red shaded area) reference frame. Two-dimensional (x, z) 
wave-driven circulation in a Lagrangian reference frame typically 
observed (b) over the inner shelf, where the onshore Lagrangian 
Stokes drift ust{z) (red) and the offshore wave-driven Eulerian 
undertow ue{z) (blue) cancel resulting in a net Lagrangian velocity 
ul{z) = 0 (black dashed line) at all depths, and (c) inside the surf 
zone, where the onshore Lagrangian Stokes drift ust{z) (red) is depth- 
uniform and the offshore wave-driven Eulerian undertow W£(z) 
(blue) is parabolic resulting in a vertical imbalance and ul{z) 4- 0 
(black dashed line). Two-dimensional (x,y) conceptual diagrams of 
(d) the conventional rip current theory, where a continuous onshore 
Stokes transport over the shoals {Qst Kip) is balanced by a steady 
offshore rip current transport through the rip channels {Qe.hp brip), 
and (e) the observed rip current cross-shore exchange, where 
surfzone rip current circulations (solid line block arrows) largely 
retain material inside the surf zone, and episodic, VLF rip current 
exits transport material offshore (offshore block arrow) and surf 
zone rip current circulations pull material back into the surf zone 
(onshore block arrow) to satisfy a mass balance, with no material 
being completely removed from the nearshore system. 
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Figure 19. Field experiment site at Monterey, California, in the local coordinate 
system, where the thin black lines are the bathymetry contours 
(labeled with the elevation in meters relative to MSL) and the thick 
black line represents the mean shoreline. In-situ ADCP locations are 
shown, where black circles are the locations for yeardays 121 to 133 
and red triangles are the locations for yeardays 134 to 138. Examples 
of the three drifter release approaches are shown from yearday 130 
as crosses: 1) in a cluster at the offshore edge of a rip channel, just 
outside the surf zone (blue), 2) in an alongshore line outside the surf 
zone, spanning multiple shoals and rip channels (magenta), and 3) in 
a cross-shore line outside the surf zone, in line with a rip channel 

(green). 
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Figure 20. Conditions measured during the field experiment of (a) signifieant 
wave height Hs, (b) peak wave period Tp, (e) depth-averaged eross- 
shore (red) and alongshore (blue) velocities measured by ADCP 3 in 
the rip channel, (d) total wind stress magnitude, and (e) tidal 
elevation relative to MSL. In (a) and (b) the black lines are measured 
by the ADCP in 13 m water depth, and gray lines are measured by 
NDBC buoy 46240 shoaled to 13 m water depth. Gray shaded 
regions represent times of drifter deployments. 
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Figure 21. Drifter positions and speed, where the color of the line represents the 
speed, for drifter deployments exhibiting the locally-contained cross¬ 
shore exchange drifter pattern on yeardays 122, 129, and 132. 
Bathymetry contours are shown in the background in black, the 
shoreline is shown as the solid black line, and the approximate 
surfzone boundary is shown as the dashed gray line. Colorbar 
represents the drifter speed. Green circles show drifter release 
locations and red circles show drifter retrieval locations. Note there 
are different cross-shore and alongshore scales between plots. 
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Figure 22. Drifter positions and speed, where the color of the line represents the 
speed, for drifter deployments exhibiting the cross-shore and 
alongshore exchange drifter pattern on yeardays 125, 126, 130, and 
137. Bathymetry contours are shown in the background in black, the 
shoreline is shown as the solid black line, and the approximate 
surfzone boundary is shown as the dashed gray line. Colorbar 
represents the drifter speed. Green circles show drifter release 
locations and red circles show drifter retrieval locations. Note there 
are different cross-shore and alongshore scales between plots. 
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Figure 23. Log-log plots of the probability density function (pdf) of the 

anomalous relative drifter position displacements outside the surf 
zone for a relative time step of L = 2 s for (top row) drifter 
deployments exhibiting the locally-contained cross-shore exchange 
drifter pattern on yeardays 122, 129, and 132, and (bottom row) 
drifter deployments exhibiting the cross-shore and alongshore 
exchange drifter pattern on yeardays 125, 126, 130, and 137. 
Contours are \ogw[P{rx, ry'; t'=2 s)] = -4.5, -4, -3.5, ... 1, 1.5. 
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Figure 24. Absolute diffusivity (one-particle) statistics for the cross-shore (Kxx, 
red) and alongshore (Kyy, blue) for drifter positions recorded outside 
the surf zone for (top row) drifter deployments exhibiting the 
locally-contained cross-shore exchange drifter pattern on yeardays 
122, 129, and 132, and (bottom row) drifter deployments exhibiting 
the cross-shore and alongshore exchange drifter pattern on yeardays 
125, 126, 130, and 137. 
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Figure 25. Drifter cross-shore velocity magnitude as a function of 

cross-shore distance, measured in surfzone-widths Lx, for each 
drifter deployment, where onshore velocities shown in 

red and offshore velocities are shown in blue, and error bars 

are shown. Black line represents the theoretical estimate of the 
Stokes drift averaged over the upper 0.5 m of the water column, 
W5((0.5 m), which correspond to the depth of the surface drifters. 
Rows 1 and 2 show drifter deployments exhibiting the locally- 
contained cross-shore exchange drifter pattern on yeardays 122, 129, 
and 132, rows 3 and 4 show drifter deployments exhibiting the cross¬ 
shore and alongshore exchange drifter pattern on yeardays 125, 126, 
130, and 137, and row 5 shows drifter deployments with no 
distinguishable exchange pattern on yeardays 135 and 136. 
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Figure 26. 3-hour mean vertical profiles, centered around low tide, of cross¬ 

shore currents measured at each ADCP, sorted by HJh and bin- 
averaged for HJh ranges given in the legend, where ADCPs 2, 3 and 
4 were located at the surfzone boundary inside the rip channel, 
ADCPs 1 and 5 were located at the surfzone boundary over shoals, 
and ADCP 6 was located outside the surf zone offshore of the rip 
channel. Note the variable velocity scales. 
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Figure 27. Normalized depth-averaged cross-shore Eulerian currents ^ 

measured at each ADCP and bin-averaged by HJh, where ADCPs 2, 
3 and 4 were located at the surfzone boundary inside the rip channel, 
ADCPs 1 and 5 were located at the surfzone boundary over shoals, 
and ADCP 6 was located outside the surf zone offshore of the rip 
channel. Solid black line represents theoretical Stokes drift that 

balances ^ according to Lentz et al. (2008), and the dashed black 
line represents theoretical ^ through the rip channel that would 

balance the continuous onshore Stokes transport over the shoals 
according to the conventional rip current theory. 
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Figure 28. Drifter observations spatially binned into 10 m by 10 m bins and 
averaged for yearday 132, with the mean eross-shore movement in 
each bin represented as a percentage of time during the drifter 
deployment that the drifters were moving seaward. Red indicates 
drifters moved seaward 100% of the time, blue indicates drifters 
only moved shoreward (moved seaward 0% of the time), and lighter 
colors indicate a tendency for drifters to move in both directions. 
Bathymetry contours shown in black and ADCP locations shown as 

black circles. 
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Figure 29. Eulerian velocity measurements at each ADCP (columns) during the 
drifter deployment on yearday 132, showing (top) the depth-relative 
mean vertical profile averaged over the duration of the drifter 
deployment, and (bottom) the surface cross-shore velocity (gray) and 
low-pass filtered (f < 0.004 Hz) surface cross-shore velocity (black), 
representing VLF rip current pulsations, where ADCPs 3 and 4 were 
located at the surfzone boundary inside the rip channel. Note the 
variable velocity scales in the top plots. 
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Figure 30. Carmel River State Beaeh, California, field site shown with a Google 
Earth image and bathymetry eontours in meters overlaid in white. 
The inset shows the mean cross-shore beach profile, relative to mean 
sea level. Red dots represent the alongshore array of stationary dye 
sensors, cyan squares represent the cross-shore array of ADCPs, and 
the green square represents the location of the EMCM array. The 
local coordinate system is shown, where positive is offshore and to 

the south. 
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Figure 31. Hourly-mean (a) root-mean square wave height, (b) mean wave 
period, Tmoi, (c) depth-averaged cross-shore velocity, (d) depth- 
averaged alongshore velocity, and (e) tidal elevation measured at 
ADCPl in 3 m water depth. Red dots and error bars in (c) and (d) 
represent surfzone currents estimated with the dye. Green lines in (c) 
and (d) represent depth-averaged velocities measured by the EMCM 
array in 0.4 m water depth. Vertical dashed magenta lines represent 
times of dye releases. Velocities are in the local coordinate system 
where positive is offshore and to the south. 
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Figure 32. Photos showing examples of (a) the initial dye slug release, (b) the 
well-mixed surf/swash zone shortly after the dye release, and (c-e) a 
dye sensor attached near the sea-bed on a pole (shown by the green 
arrows) in the surf zone, highlighting the intermittency of the dye 
observations as the sensor is alternately submerged and not 

submerged. 
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Figure 33. Example of measured dye concentration as a function of time (black) 
and the dye envelope fit to the measurements (magenta). 
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Figure 34. Example of eross-shore veloeities measured by the EMCM array for 
times when the EMCM array was (top) inside the surf zone and 
(bottom) outside the surf zone. Cross-shore veloeity is in the loeal 
eoordinate system where positive is offshore. The height of eaeh 
EMCM sensor above the sea-bed is given in the legend. 
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Figure 35. (a) Water depth measured at ADCPl (dashed line) and at EMCM 

array (solid line), with times when the EMCM array was inside of 
the surf zone shown by the thick black line, (b) half-hour mean 
cross-shore velocity profiles measured by ADCPl (open symbols) 
and EMCM array (solid symbols), and (c) half-hour mean 
alongshore velocity profiles measured by ADCPl (open symbols) 
and EMCM array (solid symbols), as a function of depth-normalized 
height above the sea-bed. Symbols in the velocity profiles 
correspond to times in the water depth time series. Velocities are in 
the local coordinate system where positive is offshore and to the 

south. 
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Figure 36. Dye eoneentration as a funetion of time (left) and as a function of 
alongshore distance from the dye release (right) measured by the 
alongshore array of stationary dye sensors for (top) yearday 162a 
during weak alongshore current conditions, and (bottom) yearday 
166c during strong alongshore current conditions. Colored circles in 
left plots represent tmoi computed with Equation (14). 
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Figure 37. Estimates based on the dye concentration time series measured by 
the alongshore stationary dye sensor array of (a) the surfzone 
averaged alongshore current V using Equation (13), (b) the 
alongshore decay of dye mass M(y) using Equation (15), and (c) the 
surfzone averaged cross-shore current U using Equation (12), as a 
function of alongshore distance y from the dye release. In (b) the 
horizontal black dashed-dotted line represents where Mo decreases 
by 1/e, and the in (b) and (c) the vertical black dashed lines represent 
the division between the near-field and far-field regions of the flow. 
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Figure 38. Dye dispersion inside of the surf zone measured with the alongshore 
array of stationary sensors as a funetion of weight-mean 
concentration time tmoi- The dashed line represents the line of best-fit 
for all measurements on yearday 166, which has an = 0.44 and is 
significant in the 95% confidence interval. The alongshore 
diffusivity Ky is equal to half of the slope of the best-fit line. 
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Figure 39. Example of data eolleeted by swimmers for dye releases (left) 172a 
and (right) 172b. (top) Repeated eross-shore transeets performed by 
different swimmers, and (bottom) the eorresponding mean dye 
coneentration profiles as a function of cross-shore distance for 5 
minute time windows. Each color represents a different swimmer, 
and solid lines represent their mean cross-shore transect with ± one 
standard deviation represented by shaded areas. In bottom plots, 
solid vertical black lines represent approximate location of the 
shoreline, dashed vertical black lines represent approximate edge of 
surf zone, and asterisks indicate the first moment of the dye 
concentration profile representing the cross-shore location of the dye 
plume /ix- The local coordinate system is used where positive is 
offshore and to the south. 
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Figure 40. Dye dispersion outside the surf zone measured with the cross-shore 
swimmer transects as a function of weight-mean concentration time 
tmoh with error bars representing the 95% confidence levels for the 
dispersion estimates for each transect. Solid lines of each color 
represent the dispersion estimates from the swimmer data for that 
dye release. Dashed lines of each color represent the lines of best-fit 
for the swimmer data for that dye release. Coefficients of 
determination are given for each swimmer, with all being significant 
in the 95% confidence interval. The cross-shore diffusivity Kx outside 
the surf zone is equal to half of the slope of the best-fit line. 
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Figure 41. A conceptual diagram of the mixing and transport of dye (pink) that 
was observed on the sandy, steep beach at CRSB. Dye released as an 
instantaneous point source at the shoreline (star at t = 0, x = 0, y = 0), 
which is assumed to be rapidly vertically well-mixed, initially 
dispersed in two dimensions in the near-field {t< 5 min, y < 50 m), 
where advection and diffusion processes were important; then when 
the dye was completely mixed across the surf zone, the dye was 
transported due to advection in the far-field {t> 5 min, y > 50 m). In 
the near-field, the dye spread in the cross-shore (kx) by turbulent 
diffusion due to breaking waves and spread in the alongshore (Ky) 
due to shear in the alongshore current, while also being advected 
downstream with the surfzone averaged alongshore current (V ) and 
also being advected offshore, resulting in a relatively large amount 
of cross-shore exchange (U(y)). In the far-field, the dye was 
advected downstream with the surfzone averaged alongshore current 
(V ) and the cross-shore transport leveled off to a constant (U). 
Outside the surfzone (x > x^), the dye moved as a constant patch, 
with a small amount of cross-shore spreading. 
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Figure 42. Frequency band-passed (top) cross-shore velocities and (bottom) 
excursion lengths computed from the time-integration of cross-shore 
velocity, measured at ADCPl for the given frequency bands. 
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